Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,591,499 C→T 5.9% intergenic (‑220/‑158) btsT ← / → tsr pyruvate:H(+) symporter/methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,591,4990CT 5.9% 426.9 / 10.5 134intergenic (‑220/‑158)btsT/tsrpyruvate:H(+) symporter/methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (65/61);  new base T (2/6);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

Only 100 of 136 total aligned reads displayed.

TAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAA  >  NC_000913/4591350‑4591648
                                                                                                                                                     |                                                                                                                                                     
tAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTc                                                                                                                                                       >  2:1135767/1‑150 (MQ=255)
   aTTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCata                                                                                                                                                    <  2:385645/150‑1 (MQ=255)
   aTTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCata                                                                                                                                                    <  1:1297061/150‑1 (MQ=255)
      gACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATca                                                                                                                                                 >  1:522675/1‑150 (MQ=255)
        ccGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATcaca                                                                                                                                               >  1:992941/1‑150 (MQ=255)
           cTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACatat                                                                                                                                            <  2:499679/150‑1 (MQ=255)
            tATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATAtt                                                                                                                                           <  1:341581/150‑1 (MQ=255)
                 tGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtatcg                                                                                                                                      <  1:667229/150‑3 (MQ=255)
                  gTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATAtt                                                                                                                                     >  1:934788/1‑150 (MQ=255)
                  gTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATAtt                                                                                                                                     >  2:1117229/1‑150 (MQ=255)
                  gTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATAtt                                                                                                                                     >  2:219242/1‑150 (MQ=255)
                  gTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtatcgt                                                                                                                                     <  2:75898/150‑4 (MQ=255)
                       gATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGaa                                                                                                                                <  2:867496/150‑1 (MQ=255)
                             gCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAAttt                                                                                                                          <  1:607835/150‑1 (MQ=255)
                             gCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAAttt                                                                                                                          >  1:871279/1‑150 (MQ=255)
                             gCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAAttt                                                                                                                          >  1:1033824/1‑150 (MQ=255)
                              ccAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTa                                                                                                                         <  2:1258767/150‑1 (MQ=255)
                              ccAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTa                                                                                                                         >  2:437850/1‑150 (MQ=255)
                              ccAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTa                                                                                                                         >  2:360253/1‑150 (MQ=255)
                                 gagTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGt                                                                                                                      <  2:655319/150‑1 (MQ=255)
                                  agTTACGCGCGCGTAAAGTAATGTAAATAACTGATTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTa                                                                                                                     <  1:1225890/150‑1 (MQ=255)
                                     tACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCa                                                                                                                                                      <  2:445862/114‑1 (MQ=255)
                                     tACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCa                                                                                                                                                      >  1:445862/1‑114 (MQ=255)
                                      aCGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATcaca                                                                                                                                               <  1:64035/120‑1 (MQ=255)
                                      aCGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATcaca                                                                                                                                               >  2:64035/1‑120 (MQ=255)
                                       cgcgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGACTCATATCACATATTTATATTGTGAATAATTTATGTAaataa                                                                                                                >  1:155117/1‑150 (MQ=255)
                                        gcgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataat                                                                                                               >  2:225606/1‑150 (MQ=255)
                                         cgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATg                                                                                                              >  1:756576/1‑150 (MQ=255)
                                          gcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATgg                                                                                                             >  2:1209345/1‑150 (MQ=255)
                                             cgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCtt                                                                                                          <  2:287649/150‑1 (MQ=255)
                                              gTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCttt                                                                                                         <  2:319005/150‑1 (MQ=255)
                                                aaaGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCttttt                                                                                                       >  1:949280/1‑150 (MQ=255)
                                                 aaGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTa                                                                                                      >  2:1156709/1‑150 (MQ=255)
                                                  aGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTaa                                                                                                     <  1:459535/150‑1 (MQ=255)
                                                   gTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTaaa                                                                                                    <  2:967253/150‑1 (MQ=255)
                                                    tAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAt                                                                                                   <  1:1413106/150‑1 (MQ=255)
                                                    tAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAt                                                                                                   <  2:940874/150‑1 (MQ=255)
                                                    tAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAt                                                                                                   >  1:109968/1‑150 (MQ=255)
                                                    tAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAt                                                                                                   >  2:760015/1‑150 (MQ=255)
                                                      aGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                 <  1:360253/150‑1 (MQ=255)
                                                       ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                >  2:747461/1‑150 (MQ=255)
                                                       ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                <  1:1135767/150‑1 (MQ=255)
                                                       ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                <  2:701323/150‑1 (MQ=255)
                                                       ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                >  2:568097/1‑150 (MQ=255)
                                                         tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                 <  1:607230/147‑1 (MQ=255)
                                                         tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                 >  2:607230/1‑147 (MQ=255)
                                                         tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                              >  2:452956/1‑150 (MQ=255)
                                                         tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                              >  2:1121675/1‑150 (MQ=255)
                                                         tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATAATTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                              >  1:1155455/1‑150 (MQ=255)
                                                               aCTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGa                                                                                        >  2:782840/1‑150 (MQ=255)
                                                                cTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGaa                                                                                       <  2:222471/150‑1 (MQ=255)
                                                                   aGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATaa                                                                                    <  2:63843/150‑1 (MQ=255)
                                                                      ggTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAAt                                                                                 <  2:1147205/150‑1 (MQ=255)
                                                                         tATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTATATTCAGAGTGTGAATAAAATTCa                                                                              >  2:823065/1‑150 (MQ=255)
                                                                          aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                             >  1:1327669/1‑150 (MQ=255)
                                                                          aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                             <  2:611436/150‑1 (MQ=255)
                                                                          aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                             >  1:423050/1‑150 (MQ=255)
                                                                          aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                             >  1:807779/1‑150 (MQ=255)
                                                                           ttttAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                              >  2:729760/1‑132 (MQ=255)
                                                                           ttttAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                              <  1:729760/132‑1 (MQ=255)
                                                                                gggATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCg                                                                       <  1:231564/150‑1 (MQ=255)
                                                                                 ggATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGt                                                                      >  2:719007/1‑150 (MQ=255)
                                                                                 ggATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGt                                                                      <  1:1070696/150‑1 (MQ=255)
                                                                                  gATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTa                                                                     >  1:1177557/1‑150 (MQ=255)
                                                                                             tCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCggg                                                          <  2:100136/150‑1 (MQ=255)
                                                                                              cAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGa                                                         >  2:553677/1‑150 (MQ=255)
                                                                                               aGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGat                                                        <  1:437850/150‑1 (MQ=255)
                                                                                                    tGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCa                                                   <  1:568097/150‑1 (MQ=255)
                                                                                                    tGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCa                                                   <  1:349515/150‑1 (MQ=255)
                                                                                                      cGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATa                                                 >  1:1244327/1‑150 (MQ=255)
                                                                                                        gTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATaaa                                               >  1:1175859/1‑150 (MQ=255)
                                                                                                        gTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATaaa                                               >  1:371790/1‑150 (MQ=255)
                                                                                                         ttGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAg                                              <  2:605245/150‑1 (MQ=255)
                                                                                                                gCAAAGGCGTCGTTTAAGTGACAATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttcc                                       >  2:743995/1‑150 (MQ=255)
                                                                                                                 cAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttcct                                      >  1:470617/1‑150 (MQ=255)
                                                                                                                     ggCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttcctttcc                                  <  2:934788/150‑1 (MQ=255)
                                                                                                                      gCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCa                                 <  1:219242/150‑1 (MQ=255)
                                                                                                                           gtTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCg                            >  2:512187/1‑150 (MQ=255)
                                                                                                                            tttAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGa                           >  2:1336572/1‑150 (MQ=255)
                                                                                                                               aaGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGaaaa                        >  2:900044/1‑150 (MQ=255)
                                                                                                                                acagACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTAAGTAAATAATGGCTTTTAAAATTCTGAGTGTGAATAAAATTGACACGGCCTAATCTCCGCGGGATATTCATAAAGTTTTTCATTTACAGGCCGAAAAt                       >  2:957183/4‑150 (MQ=255)
                                                                                                                                aGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTtttcct                                      >  2:1140896/1‑135 (MQ=255)
                                                                                                                                 gTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATc                      <  2:1430998/150‑1 (MQ=255)
                                                                                                                                 gTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATc                      >  2:557128/1‑150 (MQ=255)
                                                                                                                                  tGACGATAAATATGTGACTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                     <  2:155117/150‑1 (MQ=255)
                                                                                                                                   caataTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                     <  2:850217/145‑1 (MQ=255)
                                                                                                                                    aCGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTg                   <  1:1121675/150‑1 (MQ=255)
                                                                                                                                     cGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGc                  <  2:522675/150‑1 (MQ=255)
                                                                                                                                      gATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCa                 >  1:125896/1‑150 (MQ=255)
                                                                                                                                        tAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATc               <  2:1033824/150‑1 (MQ=255)
                                                                                                                                         aaaTATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCg              <  1:225606/150‑1 (MQ=255)
                                                                                                                                          aaTATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCgg             >  2:492337/1‑150 (MQ=255)
                                                                                                                                           atatGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGt            >  1:1134968/1‑150 (MQ=255)
                                                                                                                                            tatGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTc           <  1:512187/150‑1 (MQ=255)
                                                                                                                                              tgtgATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCca         >  1:38887/1‑150 (MQ=255)
                                                                                                                                                 gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAg      >  2:49252/1‑150 (MQ=255)
                                                                                                                                                 gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAg      >  1:536569/1‑150 (MQ=255)
                                                                                                                                                 gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAg      >  2:351015/1‑150 (MQ=255)
                                                                                                                                                 gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGATTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCAGAGGTATATTCATAAAGTTTTTCCTTTCAAGGCCGAAAATCTTGCATCTGTCCACAg      >  1:1430438/1‑150 (MQ=255)
                                                                                                                                                     cATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGaaa  >  1:471510/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAA  >  NC_000913/4591350‑4591648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: