Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC24.8% 278.5 / 112.8 129intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (17/80);  new base C (1/31);  total (18/111)
Fisher's exact test for biased strand distribution p-value = 4.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 153 total aligned reads displayed.

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCG  >  NC_000913/2947317‑2947549
                                                                                                                                                      |                                                                                    
ccAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                      <  4:1251938/151‑1 (MQ=255)
 cAAATCCGTGAGCTGGTTATAATTAAGGTTGAATGAATACGCCAGCGGAGAAGACTCCGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaa                                                                                     <  7:916425/151‑1 (MQ=14)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                    <  5:1723071/151‑1 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                    <  7:377746/151‑1 (MQ=255)
    aTCGGTGATCGTGGTAAAATTAAGGTTGCATGAAAACTAGTGCGGAGTATAGTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                  <  2:1226543/151‑3 (MQ=17)
      gggtgcgcGGGTTAAAAACAAGGTTGCATGAAAACGCGAGCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                <  6:1828395/144‑7 (MQ=18)
      ccGTGATGGGGGTAAAAAATAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                <  2:1330614/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGTGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                               <  7:1632899/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                               <  5:299828/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                               <  5:1505541/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCT‑GTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAAttgat                                                                              >  8:556481/1‑148 (MQ=255)
        gTGATCGTGGTAAAACAACGGGTGCAACAAATAGCGAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgat                                                                             <  6:1139927/151‑1 (MQ=17)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                              <  8:270593/151‑1 (MQ=255)
         ggATCGGGTTAAAATACCGGGTGCATGAAAACGCGAGCGGAGTATACTGCGCATCCCCGGACGCGGGGTGGAGCTGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacc                                                                             <  3:1531544/150‑10 (MQ=17)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                             <  7:1417587/150‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                             <  6:531709/150‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTtga                                                                              >  4:529832/1‑148 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCCGCCGGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAACTCCGGCCCCCGCACCCAA‑TAAAATTtggtggc                                                                          >  8:1484113/1‑149 (MQ=21)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                           <  8:167793/151‑1 (MQ=255)
             ggtgttAAAAAAAAGGTTGCATGAAAACTCTAGCGGAGTATAGTGCGCATACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑t                                                                       <  7:73008/146‑1 (MQ=17)
             gggggTAAAATAAAGCCTGAATGAAATCTAGAGCGGAGTATACTGCCCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑t                                                                       <  2:77562/150‑1 (MQ=14)
             gggggTAAAACAAAGGGGGCAGCAAAACCCGATCGGAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGaa                                                                          <  1:766834/148‑1 (MQ=14)
             gggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                         <  1:1283602/150‑1 (MQ=255)
             cGGGGTATAAAAAACGGTGCACGAAATCGCGAGCGGAGTATTCTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGaa                                                                          <  6:910471/149‑1 (MQ=14)
             cGGGGTAAAAAACCCGCTGCAGGAAAAAGAGAGCGGAGAATAGAGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGaa                                                                          <  1:675544/149‑1 (MQ=14)
              gtgttAAAAACAAGGTTGCATGAAAAAGAGATCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑a                                                                     <  8:529123/147‑1 (MQ=17)
              ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑t                                                                       <  7:888879/151‑1 (MQ=255)
              ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑t                                                                       <  1:1035409/151‑1 (MQ=255)
              gggcTAACAAAAAGGCTGCATGAAAACGCGAGCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑a                                                                     <  1:1715437/151‑1 (MQ=17)
               gggTAAAAAAAAGGTTGCATGAACACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                     <  5:455079/151‑1 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                     <  2:951581/151‑1 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                     <  2:698749/151‑1 (MQ=255)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaatgaacaccctaccgg                                                                      >  6:999860/1‑136 (MQ=255)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aa                                                                    <  7:119748/151‑1 (MQ=255)
                ggTAAAAAAAAGGTTGCATGAAAACGAGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTAGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCAGCAACCAATTattgaacaaactaccg                                                                       >  8:234390/1‑135 (MQ=18)
                gctAAAATTAAGGATGCATGAAAAATCAATCGTAGTATTTTTCCTCTCCACTGACGCGGTGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaa                                                                   <  5:1762681/149‑1 (MQ=11)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                   <  1:472246/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                   <  3:530777/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                   <  3:354827/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                   <  7:405216/151‑1 (MQ=255)
                 gTAAAAAAAACGTTCCATGAAAACTCTAGCGGAGAATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaa                                                                  <  6:1070165/151‑2 (MQ=17)
                   aaaaaaacGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                 <  3:959265/151‑1 (MQ=255)
                   aaaaaaaaGGTTGGATGAAAACGCGAGCGGAGTATAGTGCGCATGCACGGACGCGGGGTGTAGCAGCCTGGTAGCTCGTCGGGCTCATTACCGTAAGGTCGTCGGTGCAAATCCGGCCGCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                 <  4:1412177/151‑1 (MQ=21)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAg                                                                  <  2:956109/150‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                 <  7:830337/151‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                 <  3:920472/151‑1 (MQ=255)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATTTTGTGGAAG‑A‑AAAGc                                                                 >  1:170609/1‑150 (MQ=255)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTAGATAGCGCATACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaagca                                                               <  6:1686625/151‑5 (MQ=14)
                    aaaaaaaGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCGCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATATGATGAAG‑T‑AAAGCa                                                                <  7:618714/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAg                                                               <  6:172269/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAg                                                               <  7:849145/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAg                                                               <  8:1044083/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAg                                                               <  4:1198645/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAATTT‑ATGAAGTT‑AAAGc                                                                 >  4:1608671/1‑149 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTTGATGAAG‑T‑aaaagc                                                                >  7:586067/1‑146 (MQ=255)
                      aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGt                                                              <  6:1463318/151‑1 (MQ=255)
                       aaaaCGTTGCAGGAAAACCCGAGCGGAGAATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTa                                                             <  1:950610/151‑1 (MQ=14)
                        taagcctgaaTGAAAAAGATATCGGAGAATAGTGCGAATCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTa                                                             <  6:1225700/141‑1 (MQ=1)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCCGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTTATTAAAG‑A‑AAAGCAGTa                                                             >  5:1230388/1‑150 (MQ=18)
                        aaaGGCTGCATGAAAACGAGAGCGTAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTAc                                                            <  3:586064/151‑1 (MQ=14)
                         ccggcgccATCAAAACCCCAGGGGAGAAGACTCCCCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACg                                                           <  1:1301879/144‑1 (MQ=11)
                            gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                           >  7:1086596/1‑123 (MQ=255)
                            gTTGAATGAAAACGAGAGCGTAGTATAGTGCGAATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTg                                                        <  8:1461895/151‑1 (MQ=11)
                            gTTCCATGAAAACCCGATCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                           <  3:1087066/134‑12 (MQ=14)
                              tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTagatagg                                                       >  4:1824015/1‑144 (MQ=17)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAcg                                                     <  4:880969/151‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAcg                                                     <  4:1452841/151‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAc                                                      <  8:136942/150‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAAATTGATGAAT‑TAAAAGCAGTAc                                                            >  5:1043700/1‑144 (MQ=17)
                               gCATGAAAAATCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTAc                                                            <  6:1041342/144‑1 (MQ=14)
                                 aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAcgcg                                                   <  6:64493/151‑1 (MQ=255)
                                   gaatacgcttggggtgTATAGTGCGCCTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGACGCggg                                                 <  2:172095/136‑1 (MQ=11)
                                   gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCgg                                                  <  8:586924/150‑1 (MQ=255)
                                    aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                   <  8:1794431/118‑4 (MQ=255)
                                    aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                   >  4:1794892/1‑115 (MQ=255)
                                     aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGt                                               <  1:599126/151‑1 (MQ=255)
                                     aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGt                                               <  3:857033/151‑1 (MQ=255)
                                      aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTg                                              <  7:1587975/151‑1 (MQ=255)
                                      aaCGCGAGCGGAGTATAGCGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTg                                              <  8:604336/151‑1 (MQ=255)
                                       aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTgg                                             <  1:280707/151‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGa                                            <  3:1302821/151‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGa                                            <  5:748016/151‑1 (MQ=255)
                                         gcgAGCGGAGTATAGTGGGCATCCACGGTCGCGGGGTGGAGCAGCCTGGTAGCTCGTAGGGCTCATAACCCGAAGGTCGTCGGTTCTAATCCGGCCGCCGCAAGCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGag                                           <  4:385928/151‑1 (MQ=18)
                                         gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGAGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGag                                           <  1:1328119/151‑1 (MQ=255)
                                          cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagc                                          <  1:1748487/151‑1 (MQ=255)
                                          cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagc                                          <  2:1389954/151‑1 (MQ=255)
                                           gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagca                                         <  8:1589495/151‑1 (MQ=255)
                                             gCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagcagc                                       <  4:1627526/151‑1 (MQ=255)
                                                gAGTATAGTGCTCATCCCCGGACGCTCGGTGGAGCAGGCTGGTAGCTCGGCGGTCTCATAACCCTAATGGCGTCGGTTCCAATCCGTCCCCCGCAACCAATTAAAATTTGAAGAAT‑T‑CACGCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                    <  4:1362036/151‑1 (MQ=17)
                                                    ttaGTGCGGATCGACGGACTCGGTGGGGTGCAGCCTGGTAGCTCTGCGGTATCATAACCCTAATGGCGGCGGGTCAAATCCTGCCCCCGCAACCAATAAAAATTCGATGTAT‑T‑AACCCATTACCGGGACGCGGGGTGGAGCAGCCTGGTa                                 <  3:1145117/149‑1 (MQ=9)
                                                       cggctcATCCCCGGACGAGGGGGGGAGCAGGCTGGTCGCTGTGCGGCCAATTAACCCGAATGGCGGCGCGGCAAAACCCGCCCCCGCAACCAAATACACTGCGCTGAAG‑A‑CACGCACTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCtc                             <  3:1024949/146‑1 (MQ=255)
                                                             aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagca                                         <  8:1602769/133‑1 (MQ=255)
                                                             aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTAGGGCTCATAACACGAAGGTCGTCGGTTCAAATCAGGCCCGCGGAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagca                                         >  4:1603420/1‑133 (MQ=17)
                                                               ccACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                   >  5:1434336/1‑88 (MQ=14)
                                                                  cGGACTCGGGGTGGAGCAGCCTGGTCTCTCGTCTGGCTCTTAACCCGAAGGGGGTCGGTGCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTc                    <  3:1461014/149‑1 (MQ=14)
                                                                  cGGACGCGGGGTGGAGCAGCCTGTTAGCTCTTCGCTAACATAACCCGAAGGTCGTCTGGGCAAAACACCCCGCCGCAACCAATTAAAATTCGATGAAG‑A‑AACCCACTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCAt                  <  8:77280/151‑1 (MQ=9)
                                                                                                            aCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAATAAAT‑A‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAcc              >  2:36878/1‑113 (MQ=11)
                                                                                                                          ggTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGtcgtcg  >  3:262509/1‑111 (MQ=255)
                                                                                                                                       cccccGCAACCAATTAAAATTTCATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAcc              <  1:967893/86‑1 (MQ=255)
                                                                                                                                                      aaaaTTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                 >  7:1095873/1‑52 (MQ=40)
                                                                                                                                                      |                                                                                    
CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCG  >  NC_000913/2947317‑2947549

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