Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,562,708 A→T 7.8% intergenic (‑33/+35) yjiJ ← / ← yjiK putative transporter YjiJ/uncharacterized protein YjiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,562,7080AT 7.8% 269.0 / 8.9 87intergenic (‑33/+35)yjiJ/yjiKputative transporter YjiJ/uncharacterized protein YjiK
Reads supporting (aligned to +/- strand):  ref base A (39/41);  new base T (2/4);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GCAACATAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTATGCGAC  >  NC_000913/4562561‑4562854
                                                                                                                                                   |                                                                                                                                                  
gCAACATAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAg                                                                                                                                                  >  2:1566067/1‑150 (MQ=255)
   aCATAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTta                                                                                                                                               >  1:380044/1‑150 (MQ=255)
   aCATAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTta                                                                                                                                               >  2:1247358/1‑150 (MQ=255)
      tAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTtaata                                                                                                                                            <  2:989952/150‑1 (MQ=255)
      tAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTtaata                                                                                                                                            <  2:405790/150‑1 (MQ=255)
       aGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTtaataa                                                                                                                                           >  2:594716/1‑150 (MQ=255)
        ggCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAAc                                                                                                                                          >  1:62744/1‑150 (MQ=255)
         gCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTATTAACTTGTTAacaaca                                                                                                                                         <  1:1599622/150‑1 (MQ=255)
          cGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTATTAACTTGTTAAcaacaa                                                                                                                                        <  1:224177/150‑1 (MQ=255)
              taAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAg                                                                                                                                    >  2:1315034/1‑150 (MQ=255)
              taAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAg                                                                                                                                    >  1:685062/1‑150 (MQ=255)
                aaGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGt                                                                                                                                  >  1:1519785/1‑150 (MQ=255)
                aaGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGt                                                                                                                                  <  1:591393/150‑1 (MQ=255)
                aaGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGt                                                                                                                                  >  1:389275/1‑150 (MQ=255)
                 aGGAAGCGGCCTAAACCCATACCTAATGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGTATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTg                                                                                                                                 <  2:1048727/150‑1 (MQ=255)
                   gAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGaa                                                                                                                               <  2:317760/150‑1 (MQ=255)
                         gCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGc                                                                                                                         >  1:1112652/1‑150 (MQ=255)
                         gCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTATTAACTTGTTAACAACaaaaatgac                                                                                                                                 >  1:1622223/1‑138 (MQ=39)
                         gCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTATTAACTTGTTAACAACaaaaatgac                                                                                                                                 <  2:1622223/142‑5 (MQ=39)
                             aaaCCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGc                                                                                                                     >  1:328117/1‑150 (MQ=255)
                             aaaCCAATACCTAACGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGc                                                                                                                     >  2:526935/1‑150 (MQ=255)
                                 ccATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGt                                                                                                                 <  2:780015/150‑1 (MQ=255)
                                   aTACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCa                                                                                                                       >  1:225680/1‑142 (MQ=255)
                                   aTACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCa                                                                                                                       <  2:225680/142‑1 (MQ=255)
                                    tACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTaa                                                                                                              >  1:1315275/1‑150 (MQ=255)
                                    tACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTaa                                                                                                              <  1:781994/150‑1 (MQ=255)
                                       cTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGc                                                                                                           <  1:229510/150‑1 (MQ=255)
                                          aGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCgat                                                                                                        <  1:83722/150‑1 (MQ=255)
                                            gTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCgatga                                                                                                      >  1:1587797/1‑150 (MQ=255)
                                              cAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGAtt                                                                                                    <  2:62744/150‑1 (MQ=255)
                                               aGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTg                                                                                                   <  1:1367103/150‑1 (MQ=255)
                                                   ccAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAggg                                                                                               >  2:1453094/1‑150 (MQ=255)
                                                     aGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGt                                                                                             <  1:526935/150‑1 (MQ=255)
                                                        aTCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTaaa                                                                                          >  2:803604/1‑150 (MQ=255)
                                                          cccGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGc                                                                                        <  1:92760/150‑1 (MQ=255)
                                                            cGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGa                                                                                      >  2:1355750/1‑150 (MQ=255)
                                                               aaaGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATaa                                                                                   <  2:1560498/150‑1 (MQ=255)
                                                                 aGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATaaaa                                                                                 <  2:389275/150‑1 (MQ=255)
                                                                     cgGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGa                                                                             <  1:975387/150‑1 (MQ=255)
                                                                       gTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGAtt                                                                           >  2:1326260/1‑150 (MQ=255)
                                                                        tGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGAttt                                                                          <  2:1187609/150‑1 (MQ=255)
                                                                         gCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTg                                                                         >  1:599372/1‑150 (MQ=255)
                                                                          cTGAAAGAAAAGCGTTCTACGGGATGCGTGTTCGAATGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAATTTAATAACAAAAAGGTGAAGCATTCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTgg                                                                        <  2:277050/150‑1 (MQ=255)
                                                                          cTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTgg                                                                        <  1:1228777/150‑1 (MQ=255)
                                                                           tGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGt                                                                       <  1:1278574/150‑1 (MQ=255)
                                                                            gaaagaaaAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGtt                                                                      >  1:286154/1‑150 (MQ=255)
                                                                            gaaagaaaAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGtt                                                                      <  2:1315275/150‑1 (MQ=255)
                                                                             aaagaaaAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTc                                                                     <  2:328117/150‑1 (MQ=255)
                                                                                gaaaAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCt                                                                  <  2:685062/150‑1 (MQ=255)
                                                                                  aaaGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTaa                                                                >  1:1358379/1‑150 (MQ=255)
                                                                                   aaGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACGAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAAc                                                               <  2:1499304/150‑1 (MQ=255)
                                                                                   aaGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAATCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAAc                                                               >  1:193189/1‑150 (MQ=255)
                                                                                      cGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGat                                                            <  2:240844/150‑1 (MQ=255)
                                                                                       gTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGata                                                           >  1:738036/1‑150 (MQ=255)
                                                                                          cTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATaa                                                        >  2:1144322/1‑150 (MQ=255)
                                                                                           tACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATaaa                                                       <  1:967105/150‑1 (MQ=255)
                                                                                            aCGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAat                                                      <  1:1144862/150‑1 (MQ=255)
                                                                                             cGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAata                                                     <  1:376950/150‑1 (MQ=255)
                                                                                               ggATGCGTGGACGAAGGCATGTTTAACCTTGTTTCTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGAGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTAGATAACGATATAAATAtt                                                   >  1:279125/1‑150 (MQ=255)
                                                                                                gATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTa                                                  >  2:111471/1‑150 (MQ=255)
                                                                                                 aTGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTAc                                                 <  1:594716/150‑1 (MQ=255)
                                                                                                  tGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTAcc                                                <  1:1566067/150‑1 (MQ=255)
                                                                                                  tGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTAcc                                                <  1:194154/150‑1 (MQ=255)
                                                                                                    cGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCag                                              >  1:1637934/1‑150 (MQ=255)
                                                                                                      tGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCagag                                            <  1:812876/150‑1 (MQ=255)
                                                                                                       ggACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAgg                                           <  2:1112652/150‑1 (MQ=255)
                                                                                                        gACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGc                                          >  1:1397525/1‑150 (MQ=255)
                                                                                                                  tGTTTAACCTTGTTTGTCATTTTTCTAGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCAATCGCa                                >  1:1664665/1‑150 (MQ=255)
                                                                                                                   gTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAAC‑AAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAAc                              <  1:1247358/150‑1 (MQ=255)
                                                                                                                   gTTTAACCTTGTTTGTAATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCaa                               >  1:1425547/1‑150 (MQ=255)
                                                                                                                      tAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAAccc                            >  1:52805/1‑150 (MQ=255)
                                                                                                                                 gtCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGt                 >  2:1584191/1‑150 (MQ=255)
                                                                                                                                 gtCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACACCTTCGGCCTGt                 >  2:932512/1‑150 (MQ=255)
                                                                                                                                  tCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGtt                <  1:1144322/150‑1 (MQ=255)
                                                                                                                                   cctttttGTTATTAACTTGTTAACAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGatat                                                          <  2:409825/105‑1 (MQ=38)
                                                                                                                                   cctttttGTTATTAACTTGTTAACAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGatat                                                          >  1:409825/3‑107 (MQ=38)
                                                                                                                                   cATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGTCATCCATCCCAACCCCTTCGGCCTtttt               >  2:1556415/1‑150 (MQ=255)
                                                                                                                                     tttttGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGa             <  2:256302/150‑1 (MQ=255)
                                                                                                                                      ttttGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGAt            <  2:430263/150‑1 (MQ=255)
                                                                                                                                       tttGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATg           <  1:803604/150‑1 (MQ=255)
                                                                                                                                       tttGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATg           <  2:821803/150‑1 (MQ=255)
                                                                                                                                       tttGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATg           >  1:122743/1‑150 (MQ=255)
                                                                                                                                           ttgttAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTAt       <  2:1594171/150‑1 (MQ=255)
                                                                                                                                             gttAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTATGc     <  1:1315034/150‑1 (MQ=255)
                                                                                                                                               tAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTATGCGa   >  1:809877/1‑150 (MQ=255)
                                                                                                                                               tAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTATGCGa   >  1:47971/1‑150 (MQ=255)
                                                                                                                                                aaCAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTATGCGAc  >  2:626795/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAACATAGGCGTATAAAGGAAGCGGCCTAAACCCATACCTAAGGTCAGAACCAGCATCCCGAAAAGCGCGGTGCTGAAAGAAAAGCGTTCTACGGGATGCGTGGACGAAGGCATGTTTAACCTTGTTTGTCATTTTTGTTGTTAACAAGTTAATAACAAAAAGGTGAAGCATGCCAGCCTGTTAATGCGATGATTGAGGGGTAAAGCGATAAAAACGATTTGGTTCGCTAACGATATAAATATTACCAGAGGCATCCATCGCAACCCCTTCGGCCTGTTTGATGTTATGCGAC  >  NC_000913/4562561‑4562854

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: ‑