Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,364 G→T 9.8% intergenic (+21/+26) slyB → / ← slyA outer membrane lipoprotein/global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,3640GT 9.8% 532.8 / 14.0 169intergenic (+21/+26)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (106/41);  new base T (3/19);  total (109/60)
Fisher's exact test for biased strand distribution p-value = 2.07e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

TGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720216‑1720513
                                                                                                                                                    |                                                                                                                                                     
tGGAAATTCGTAAAGACGATGGTAATACAATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGt                                                                                                                                                      >  2:1272168/1‑150 (MQ=255)
 ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                                                                                                                                                     >  1:916023/1‑150 (MQ=255)
 ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                                                                                                                                                     >  1:2832755/1‑150 (MQ=255)
 ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGATACTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAGGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTAAGCGTGTGGTCAGGtt                                                                                                                                                     >  1:1563423/1‑150 (MQ=255)
    aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                                                                                                                                                  <  2:2172139/150‑1 (MQ=255)
     aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTg                                                                                                                                                 <  1:2576588/150‑1 (MQ=255)
       tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:619327/150‑1 (MQ=255)
         gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                             <  2:1995728/150‑1 (MQ=255)
          tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcac                                                                                                                                            >  1:3620005/1‑150 (MQ=255)
           caccatgATGGTATTACCATCATGGTGGTACAGAAACAATGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca                                                                                                                                           >  1:1527807/7‑150 (MQ=255)
           aaaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACcaca                                                                                                                                           <  1:593876/150‑1 (MQ=255)
             aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACg                                                                                                                                         <  2:3807071/150‑1 (MQ=255)
                  aTGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACTccccc                                                                                                                                    >  2:1139854/1‑150 (MQ=255)
                  aTGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGccccc                                                                                                                                    >  2:2458177/1‑150 (MQ=255)
                      tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa                                                                                                                                >  2:2890617/1‑150 (MQ=255)
                       aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAt                                                                                                                               >  2:3161585/1‑150 (MQ=255)
                            catcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          >  1:3190594/1‑150 (MQ=255)
                              tcttggtggTAAAGAAACAAGCAAACACTGGTTTCTCTCGGGGCAAAGTTTTGGTAGTGGCAAGCATTGGCTGTCTGGTGACCTTTACTCCGCGCTAATTATTTTGGGGTGTGTTAAGTTTCCTGACCACACGCCCCCTTCATTTCAccc                                                                                                                        <  2:471059/147‑1 (MQ=255)
                              tcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAccc                                                                                                                        >  1:414185/1‑150 (MQ=255)
                                    tggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGc                                                                                                                  >  1:2861049/1‑150 (MQ=255)
                                      gcaCAGACACACGCCAACCCTCGTTTCTCTCCGGGCCAACGTGGCGTCCTGGCCAGCAAGGGCCGTCAGGGGACCGTTTCTCCGCGAAAAATACGTTGGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCt                                                                                                                <  2:3830464/148‑1 (MQ=255)
                                        aCAGAACCACGCCAACAGTGTTTTCTCGCGGGGCCACGGTGTGGTAGGGGCCACAATTGGCAGTCTGGTGACGGTTTCCCGGCGCTAAATATGTTGGGGTGTGGTCTGTTAACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                              <  2:3706214/150‑1 (MQ=255)
                                          aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  1:1752493/1‑150 (MQ=255)
                                          aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  2:1728376/1‑150 (MQ=255)
                                          aGAAACAAGGAAACACTCTTTTCTCTCCGGGCCAACGTTTCGTACTGGCCTGCAATGGCTGTCTGGTGACCGTTTCTCCGCGATAAATATTTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            <  2:2903121/150‑1 (MQ=255)
                                          aGAAACAAGGAAACACTCGTTTCTCTCCGGGCAAACGTGTCTTACTGGCCAGAATTGGCTGTCAGGTGACCTTTTCTCGGCGATAAATAAGTTTGCGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            <  2:1617889/150‑1 (MQ=255)
                                           gAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAAc                                                                                                           >  2:1794/1‑150 (MQ=255)
                                            aaaCATGCCAACACTGTTTTCTCTCGGGGCCAAGGTGTCTTAGTGGCATGCATTGGCTGTCTGGTTACGTTTTCTCGGCGCTAATTATTTTTGGGTGTGGTCCTGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACt                                                                                                          <  2:28040/150‑1 (MQ=255)
                                             aaCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGTCCACACGCCCCCTTCATTTCAGCCTTTGGCCTGTAACTc                                                                                                         >  2:1744695/1‑150 (MQ=255)
                                              aCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCa                                                                                                        >  1:3893930/1‑150 (MQ=255)
                                                  ggCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATACGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGa                                                                                                    <  2:3538557/150‑1 (MQ=255)
                                                  ggCAACACTCGTTTCTCTCCGGCCCAACGTGTCGTACTGGCCAGAAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTCCGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCATTGGCCTGTAACTCAATGa                                                                                                    >  2:512278/1‑150 (MQ=255)
                                                    cAACACTCTTTTCTCTCGGGGCCAAGGTGTCGTAGTGGCCAGAATTGGCTGTCTGGTGCCCGTTTCTCCGCGATAATTACGTTTGCGTGTGGTCTGTTTACTGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGata                                                                                                  <  2:3754352/150‑1 (MQ=255)
                                                    cAACACTCGTTTCCCGCGGGCCCAAGGTGTCTTCCGGGCCAGCAAGGGTTGTCCGGGGACCTTTGCCCCGCGCAAATTACGTTGGGGTGTGGACCGGACCCTCACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGata                                                                                                  <  2:1426614/150‑1 (MQ=255)
                                                      acacTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtatt                                                                                                >  1:1826589/1‑150 (MQ=255)
                                                        acTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtattat                                                                                              >  1:3150447/1‑150 (MQ=255)
                                                           cGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCt                                                                                           >  2:700765/1‑150 (MQ=255)
                                                             tttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCGTCATTTCACCCTTTGGCCTGTAACTCAATGATATTCTGCTCa                                                                                         >  1:3021079/1‑150 (MQ=255)
                                                                ctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGt                                                                                      >  1:3453315/1‑150 (MQ=255)
                                                                 tctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtt                                                                                     >  1:3786156/1‑150 (MQ=255)
                                                                   tcCGGGCCACCGTGTCGTAGTGGCCAGCATGGGCTGTCTGGTGACCGTTTCTCCGCGATAAATAAGTTTGCGTGTGGTCCGGTTACTGCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   <  2:1603530/150‑1 (MQ=255)
                                                                   tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:2563053/1‑150 (MQ=255)
                                                                   tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:572094/1‑150 (MQ=255)
                                                                   tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  1:3234008/1‑150 (MQ=255)
                                                                   tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:3894311/1‑150 (MQ=255)
                                                                   tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:3892432/1‑150 (MQ=255)
                                                                   tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTAGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:3894166/1‑150 (MQ=255)
                                                                    ccGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTg                                                                                  <  2:2624809/150‑1 (MQ=255)
                                                                     cGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGc                                                                                 >  2:852332/1‑150 (MQ=255)
                                                                     cGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGAGGACAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAAAGATATTATGCTCAAGTTTTGc                                                                                 >  1:2332665/1‑150 (MQ=255)
                                                                       ggCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGa                                                                               >  1:1263919/1‑150 (MQ=255)
                                                                         ccAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATg                                                                             >  2:3457630/1‑150 (MQ=255)
                                                                         ccAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATg                                                                             >  2:3457547/1‑150 (MQ=255)
                                                                         ccAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATg                                                                             >  2:1046182/1‑150 (MQ=255)
                                                                         ccAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATg                                                                             >  1:2838766/1‑150 (MQ=255)
                                                                              gtgtCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGAGt                                                                        >  1:1106904/1‑150 (MQ=255)
                                                                                    tACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAg                                                                  >  1:1198981/1‑150 (MQ=255)
                                                                                      cgggCCAGCCCTGGCTGTCTGGTGACGGTTTCTCCGCGCTAAATAAGTTTGCGTTTGGCCCGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGtt                                                                <  2:2861038/148‑1 (MQ=255)
                                                                                          ccAGCATTGGCTGTCTGGTGACCGTTTCTCCGCGCTAATTATGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                            <  1:2597252/150‑1 (MQ=255)
                                                                                          caaGCATTGGCTGTCTGGTGACCGTTTCTCCGCGATAAATATGTTTGCGTGTGGTCTGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                            <  1:3554800/148‑1 (MQ=255)
                                                                                           cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           >  2:152360/1‑150 (MQ=255)
                                                                                           cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           >  2:1215407/1‑150 (MQ=255)
                                                                                               aaTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGtt                                                       >  1:3363927/1‑150 (MQ=255)
                                                                                                aTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTc                                                      >  2:1288387/1‑150 (MQ=255)
                                                                                                 tGGCTGTCTGGTGACCGTTTCTCCGCGCTACATTAGTTTGCGTGTGGTACGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTcc                                                     <  2:952899/150‑1 (MQ=255)
                                                                                                 tGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTcc                                                     >  1:2204541/1‑150 (MQ=255)
                                                                                                     aGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTg                                                 >  2:935757/1‑150 (MQ=255)
                                                                                                      gTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCGTTGGCCTGTAACTCAATGATATTATGCTCAAGGTTTGCGATGAGCGTAATCAGTTGCTCCAGGTCCTCTGc                                                >  2:427672/1‑150 (MQ=255)
                                                                                                           gcccccGTTCCCCCGCTCTCAATAAGGTTGCGTGTGGTCGTTTCACGCACCCCACGCCCCCATCCATTCACCCTTTGGGCTGTAACACCATGATATTATGCACCAGTTTTGCGTTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           <  1:1346156/146‑1 (MQ=255)
                                                                                                           gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           >  2:3340375/1‑150 (MQ=255)
                                                                                                             gccgttttCTCCGCGCTAATTATGTTTGGGTGTGGTCTGTTTTGTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCGCATGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                         <  2:1777911/145‑1 (MQ=255)
                                                                                                             gACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                         <  2:2255726/150‑1 (MQ=255)
                                                                                                             gACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                         >  1:1028838/1‑150 (MQ=255)
                                                                                                                ctttcccccGCGATAAATATGTTTGGGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                      <  1:1943307/143‑1 (MQ=255)
                                                                                                                cGTTTCTCCGCTCAAAATATTTTTGCTGGTGGTCCGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                      <  2:3347485/150‑1 (MQ=255)
                                                                                                                    tctcCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTa                                  >  2:3174983/1‑150 (MQ=255)
                                                                                                                      tcCTCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGATCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTATGCGGAGATGACATGTAat                                >  2:2393838/1‑150 (MQ=255)
                                                                                                                      tcCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCAGGTAat                                >  2:3580485/1‑150 (MQ=255)
                                                                                                                        cgcgCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCCGTTCCTCTGCGGAGATGCCAGGTAatat                              >  2:186883/1‑150 (MQ=255)
                                                                                                                         gcgcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAtt                             >  2:2671383/1‑150 (MQ=255)
                                                                                                                                  atttttGGGTGTGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgg                    <  2:3762354/147‑1 (MQ=255)
                                                                                                                                  aaGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgg                    <  1:2279668/150‑1 (MQ=255)
                                                                                                                                  aaGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgg                    <  2:3232780/150‑1 (MQ=255)
                                                                                                                                     tttGCGGTTGGTCCTGTTACGGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgggct                 <  1:3229952/150‑3 (MQ=255)
                                                                                                                                        ggttttgGTCTGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg              <  1:1076938/145‑1 (MQ=255)
                                                                                                                                        gggtgtGGTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg              <  2:1263919/148‑1 (MQ=255)
                                                                                                                                          gtgtGGTCAGGTTACTGACCACAGGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGGGATGAGCGTAATCTGTTGCTCCAGTTCCTCTGCGGGGATGCGATGTATTATTTCCGCGCGGGTTTTGtt            >  2:805660/1‑150 (MQ=255)
                                                                                                                                            ggggTCCGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaa          <  2:2327388/148‑1 (MQ=255)
                                                                                                                                             tgGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  1:347868/150‑1 (MQ=255)
                                                                                                                                             tgGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTCACTCAATGATATTATGCTCAAGTTTTTCGATGAGCGTAATAAGTTGCTCAAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         >  1:3699511/1‑150 (MQ=255)
                                                                                                                                               gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       >  2:1149833/1‑150 (MQ=255)
                                                                                                                                                tccggTTACGGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc      <  2:3111462/147‑1 (MQ=255)
                                                                                                                                                tCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTCATATTTCCGCGCGGGTTTTGTTAATAAc      >  1:2797812/1‑150 (MQ=255)
                                                                                                                                                tCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc      >  1:2088935/1‑150 (MQ=255)
                                                                                                                                                tCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc      >  2:3300516/1‑150 (MQ=255)
                                                                                                                                                 cttgTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa     <  2:3033354/147‑1 (MQ=255)
                                                                                                                                                 cAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa     <  2:866633/150‑1 (MQ=255)
                                                                                                                                                  aGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAg    <  1:887247/150‑1 (MQ=255)
                                                                                                                                                    gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  2:2721138/1‑150 (MQ=255)
                                                                                                                                                    gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  2:2842069/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720216‑1720513

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑