Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,591,500 A→G 6.1% intergenic (‑221/‑157) 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,5000AG 6.1% 400.3 / 10.7 131intergenic (‑221/‑157)btsT/tsrpyruvate:H(+) symporter/methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base A (65/58);  new base G (2/6);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

ATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAG  >  NC_000913/4591353‑4591649
                                                                                                                                                   |                                                                                                                                                     
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)
                  gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATtgtg                                                                                                                                   <  1:987912/150‑1 (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:360253/1‑150 (MQ=255)
                           ccAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTa                                                                                                                          >  2:437850/1‑150 (MQ=255)
                              gagTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGt                                                                                                                       <  2:655319/150‑1 (MQ=255)
                               agTTACGCGCGCGTAAAGTAATGTAAATAACTGATTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTa                                                                                                                      <  1:1225890/150‑1 (MQ=255)
                                  tACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCa                                                                                                                                                       >  1:445862/1‑114 (MQ=255)
                                  tACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCa                                                                                                                                                       <  2:445862/114‑1 (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)
                                              aaGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTa                                                                                                       <  2:1132547/150‑1 (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                                                                                                 <  1:1135767/150‑1 (MQ=255)
                                                    ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                 >  2:568097/1‑150 (MQ=255)
                                                    ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                 <  2:701323/150‑1 (MQ=255)
                                                    ggTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                                 >  2:747461/1‑150 (MQ=255)
                                                      tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                  >  2:607230/1‑147 (MQ=255)
                                                      tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTc                                                                                                  <  1:607230/147‑1 (MQ=255)
                                                      tAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                               >  2:452956/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:423050/1‑150 (MQ=255)
                                                                       aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                              >  1:1327669/1‑150 (MQ=255)
                                                                       aTTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCAc                                                                              <  2:611436/150‑1 (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:371790/1‑150 (MQ=255)
                                                                                                     gTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATaaa                                                >  1:1175859/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)
                                                                                                                    cgtcgtTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAg                                 <  2:109968/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:557128/1‑150 (MQ=255)
                                                                                                                              gTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATc                       <  2:1430998/150‑1 (MQ=255)
                                                                                                                               tGACGATAAATATGTGACTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                      <  2:155117/150‑1 (MQ=255)
                                                                                                                                caataTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                      <  2:850217/145‑1 (MQ=255)
                                                                                                                                caataTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt                      >  1:850217/5‑149 (MQ=255)
                                                                                                                                 aCGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTg                    <  1:1121675/150‑1 (MQ=255)
                                                                                                                                  cGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGc                   <  2:522675/150‑1 (MQ=255)
                                                                                                                                   gATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCa                  >  1:125896/1‑150 (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:351015/1‑150 (MQ=255)
                                                                                                                                              gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAg       >  1:536569/1‑150 (MQ=255)
                                                                                                                                              gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAg       >  2:49252/1‑150 (MQ=255)
                                                                                                                                              gATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGATTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCAGAGGTATATTCATAAAGTTTTTCCTTTCAAGGCCGAAAATCTTGCATCTGTCCACAg       >  1:1430438/1‑150 (MQ=255)
                                                                                                                                                  cATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGaaa   >  1:471510/1‑150 (MQ=255)
                                                                                                                                                   atatCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAag  >  1:806448/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAG  >  NC_000913/4591353‑4591649

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: