Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC30.9% 242.8 / 146.7 136intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (8/86);  new base C (2/40);  total (10/126)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 189 total aligned reads displayed.

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGG  >  NC_000913/2947317‑2947613
                                                                                                                                                      |                                                                                                                                                    
ccAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                      <  8:1700290/151‑1 (MQ=255)
 cAAATGCGTGCTGGGGTTAAAAACACGGGTGCAGCAAAACCAGAGCGGAGTATAGTGCGCATCCCCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaa                                                                                                                                                     <  7:1237986/151‑1 (MQ=17)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                                     >  3:577907/1‑150 (MQ=255)
   aaTGCGTGATGGGGTTAAAAACAACGCTCAATGAAAACGCGAGAGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                                   <  7:1804374/151‑2 (MQ=17)
   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                                   >  1:1009091/1‑148 (MQ=255)
   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaatt                                                                                                                                                   >  2:456442/1‑149 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                                   >  4:986377/1‑147 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                                  <  7:229604/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                                                                  >  6:1515980/1‑149 (MQ=255)
     tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAAttt                                                                                                                                                >  6:215759/1‑151 (MQ=255)
      gcGTGATCGGGGTAAAATTAAGGTTGCATGAAATCGCGAGCGGAGTATAGTGCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGTCCCCGCAACCAATCAAAttt                                                                                                                                                 <  5:1397494/149‑1 (MQ=18)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                                                 <  2:1506615/150‑1 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGAATGAAAACTCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                                 <  5:215992/150‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAACGGGGGAATGAAAAAGCGAGCGGAGTATACTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtga                                                                                                                                              <  6:1476157/151‑1 (MQ=18)
        gTGATCTGGGTAAAAAAAAGCGTGCAGGAAAACGCGAGCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgat                                                                                                                                             <  6:358598/151‑1 (MQ=18)
        gTGATCGGTTT‑AAAAAAAGGTTGCATGAATACTATATCGGAGTATATTGCGAATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgatg                                                                                                                                            <  2:198834/151‑1 (MQ=14)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacac                                                                                                                                              <  4:1099270/151‑9 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacac                                                                                                                                              >  1:363278/1‑143 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaaca                                                                                                                                               >  3:1098512/1‑143 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                              <  7:1581678/151‑1 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                              <  6:1148971/151‑1 (MQ=255)
         tGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacc                                                                                                                                             >  5:539421/1‑142 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAACCGAAGGACGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTtga                                                                                                                                              >  6:563370/1‑148 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                           <  4:431593/151‑1 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                           <  7:92778/151‑1 (MQ=255)
           aTCGGGGAAAAATACAGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                                                                                           <  7:839504/151‑12 (MQ=17)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                          >  2:155758/1‑138 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                                     >  3:447584/1‑138 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTagataggaagagcc                                                                                                                                         >  7:1453134/1‑138 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                                          <  3:826167/150‑1 (MQ=255)
             cGGGCTAAAAAAAAGGGGGAATGAAAAAGCTATGGGAGTATAGTGCGCATCCCCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                          <  1:156147/150‑13 (MQ=14)
              gggtTAAAAAAAAGGTTGCATGAAATATCTAGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑a                                                                                                                                     <  2:190021/151‑1 (MQ=17)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑a                                                                                                                                     <  1:1228472/151‑1 (MQ=255)
                gttAAAAACAAGGTTGCATGAAAACCCGAGCGTAGTATATTTGGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  5:682558/149‑1 (MQ=14)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aa                                                                                                                                    <  4:1548986/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  8:1131306/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  4:1345237/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  6:127764/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  7:731604/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  4:707130/151‑1 (MQ=255)
                  tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                                   <  6:503709/150‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAg                                                                                                                                  <  1:1851377/150‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGc                                                                                                                                 <  1:820042/151‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACCCGAGCAGAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagc                                                                                                                                <  6:1679522/151‑4 (MQ=17)
                    aaaccacgggggCAGGAAAACCCGAGCGCAGTATAGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagca                                                                                                                               <  6:1006755/140‑5 (MQ=14)
                    aaaataaGGTTGCATGAAAAATCGATCGTTGTATTTTTCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagca                                                                                                                               <  8:1188797/151‑5 (MQ=11)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                                                <  4:1453647/151‑1 (MQ=255)
                    aaaaaaaGCGTGCATGAAAAAGCGAGCGCAGTATAGTGCCCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagca                                                                                                                               <  6:468671/151‑5 (MQ=14)
                     aatcaaGGTTGCATGAAAACGATAGTGTAGTATAGTGGGAATCCTCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagcag                                                                                                                              <  6:1139321/147‑6 (MQ=11)
                     aaataaGGTTGAATGAAAACCATAGCGTAGTATAGTCGGAATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagcag                                                                                                                              <  5:1207617/151‑6 (MQ=11)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCATCCTGGTAGCTCGTCGGGCTCATAACCAGAAGGTCGTCGGTTCAAATACGGACCCCGCAACCAATTAAAATTTAATGAA‑GT‑AAAGc                                                                                                                                 >  5:1839480/1‑149 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAg                                                                                                                               <  2:1607006/151‑1 (MQ=255)
                       aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGAGGCGTGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTa                                                                                                                             <  5:368807/151‑1 (MQ=255)
                       aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTa                                                                                                                             <  6:681101/151‑1 (MQ=255)
                       aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTa                                                                                                                             <  3:1427519/151‑1 (MQ=255)
                       aaaaGGTTGAATGAAAAATAGATAGTAGTATAGTTCGCTTACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTa                                                                                                                             <  3:761016/151‑1 (MQ=11)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTAc                                                                                                                            <  5:778386/151‑1 (MQ=255)
                        aaaGGGTGCAGGAAAAAGCGAGCGCAGTATAGTCCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTAc                                                                                                                            <  4:255631/151‑1 (MQ=14)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACg                                                                                                                           <  5:207244/151‑1 (MQ=255)
                         aaGGTTGCATGAAAAAGCGATGGGAGTATACTCCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACg                                                                                                                           <  5:1517456/151‑1 (MQ=11)
                          acggtcaatcaaaaAGCGTTCGTAGTATTGAGGGCATCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACgg                                                                                                                          <  7:1334045/141‑1 (MQ=14)
                           ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                          >  4:1867948/1‑124 (MQ=255)
                           ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGt                                                                                                                         <  1:176679/151‑1 (MQ=255)
                           ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATTAATTT‑AAAGCAGTACg                                                                                                                           >  8:363512/1‑149 (MQ=17)
                           ggTTGAATGAAAACGCGATCGGAGTATAGTGCGCATACTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                          <  3:1867862/136‑13 (MQ=14)
                            ggtGCAGGAATACGCGTGCGGAGTATAGTCCGCCTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTg                                                                                                                        <  1:951649/149‑1 (MQ=11)
                            gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTg                                                                                                                        <  8:931369/151‑1 (MQ=255)
                            gTTGCAGGAAAACTCGAGGGCTGTATAGAGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTg                                                                                                                        <  3:665781/151‑1 (MQ=11)
                              tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaacca                                                                                                                                           >  8:261393/1‑121 (MQ=255)
                              tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGa                                                                                                                       <  7:1198469/150‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGAcg                                                                                                                     <  7:307241/151‑1 (MQ=255)
                               caaggAAAAAGCCATCGGAGAATAGAGCCCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcg                                                                                                                     <  5:1681917/147‑1 (MQ=11)
                                 aTGAAAACGCGAGCGGAGTATAGTGCGCATGCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGAcgcg                                                                                                                   <  1:1238117/151‑1 (MQ=255)
                                 aTGAAAACGAGATCGTAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcgcg                                                                                                                   <  7:1609489/151‑1 (MQ=11)
                                   gNAAACGCGAGNGGAGTATANTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCggg                                                                                                                 <  2:1278744/151‑1 (MQ=255)
                                   gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCggg                                                                                                                 <  4:568458/151‑1 (MQ=255)
                                   gAAAAAGCTAGCGGTGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGACGCggg                                                                                                                 <  1:119329/151‑1 (MQ=11)
                                     aaactatagggtaGTATTGTGCGAATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGACGCGGGGt                                                                                                               <  4:978051/139‑1 (MQ=11)
                                     aaaCGCGAGCGGAGTATAGTGTGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                               <  8:316427/151‑1 (MQ=255)
                                     aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                               <  7:375662/151‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  5:1726649/151‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  4:332946/151‑1 (MQ=255)
                                           gAGCGGAGTATAGTGCTCATTCACTGACGAGTGTTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGagca                                                                                                         <  3:92182/151‑1 (MQ=17)
                                           gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaaca                                                                                                                                               <  4:178869/115‑8 (MQ=255)
                                           gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGagca                                                                                                         <  4:497893/151‑1 (MQ=255)
                                            aggggAGTATAGTGAGCATCCACGGACGCGCGTTGGAGCAGGCTTGTAGCTCGTCGGTCTCATAACCCGAAGGTCTGAGGTTCAAATCCGGCCCCCGCAACCAAATAAAAATTGATGAA‑GT‑AAAGCAGTACTGTGACGCGGGGTGGagcag                                                                                                        <  3:652103/148‑1 (MQ=17)
                                             gCGGAGTATACTCAGCTTCTACGGACTCCGGTTGGAGCAGCCTTTTATCTTGTCTGTCTCATATGCCGAAGGGCGTCGGTTCATATAATACCCGCGGAACCAAATAAATTTTGTTGAA‑GT‑TAAGCATTACGGTGACGCGGGGTGGagcagc                                                                                                       <  2:1537261/151‑1 (MQ=9)
                                             gCGGAGAATAGTCCGCATCCTCGGACTCTGTGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTGGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAAATTGATGAA‑GT‑AACGCAGTACGGTGACGCGGGGTGGagcagc                                                                                                       <  3:1747802/151‑1 (MQ=17)
                                              ggtaGTATAGTGCTCATCCACGGACGCGGGGTGGAGCAGCCTGTTAGCTCGTCGGTCTCATAACCCGAAGGTGGTCGGTTCTAATCCTTCCCCGGCAACGAATTAAAATTTGATGAA‑GT‑AAAGGAGTACGGTGACGCGGGGTGGAGCAGcc                                                                                                      <  8:708412/148‑1 (MQ=17)
                                               gtaGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaa                                                                                                                                                 <  5:38541/107‑6 (MQ=255)
                                                gAGTTTAGTCCCAATCCACGGACTCGGGGTGGAGAAGCGTGGTTGGTGGTCGGTCTAATAACCCTAAGGGCGTCGGTTCAAAACCGGCCCCCGCAACCAATTAAAATTTGATGAA‑TT‑AAACCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                                                                                    <  1:577382/151‑1 (MQ=14)
                                                gAGTATAGAGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                                                                                    <  8:1382289/151‑1 (MQ=255)
                                                   tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAAT‑AA‑GTAAAAGCAGTAc                                                                                                                            >  5:1474654/1‑124 (MQ=14)
                                                   gataGTCCGCATCCACGGACGCGGGGGGGAGCAGCCTGTTAGCTCGGCGGGCTCATAACCCGAAGGTCGGCGGGTCAAATCACCCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AACGCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTa                                                                                                 <  1:332530/150‑1 (MQ=17)
                                                    ataGTGCGCATGGACTGACGCGGTGTGGAGCAGCGTGGTTGCTGGGGGGGATCATAACCCGAAGGTCGTCGGTTCAAATCCCGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑CATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                                                <  6:301771/151‑1 (MQ=14)
                                                    ataGTGAGCATCCACCGACGCGGGGGGGATCAGCCTGGTAGCTCGTCGGGATCAAAACCCGAAGGGCGGCGGTGCAAAAACTCCCCCCGCAACCAATTAAAAATTGATGAA‑GT‑AACCCAGTCCGCTGACGCGGGGTGGAGCAGCCTGGTAg                                                                                                <  6:853034/151‑1 (MQ=17)
                                                               gcacgcaacacgtgTGGAGCAGCCTGGTCGCTCGGCTCTATAACAACCCGAATGTCTGCTCGTCAAATCCCCCCCCCGCAAGCAAAAAAAAATCGCAGAA‑GA‑CACGCCCTACGCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCt                                                                                     <  4:1302709/138‑1 (MQ=2)
                                                                      accgGTGTGGAGTAGGGTGGTTGGTCGGCGGGATAACAACCCGAAGGGCGGCTGTGGAAATCAGCCCCCCTGAACCAATAAAATATTGCAGTA‑GT‑CACCCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAcc                                                                              <  7:675114/149‑1 (MQ=2)
                                                                            ttGGAGCAGCCTGTTAGCTCTTCGGTCTCATAACCCTAAGGTCGTCGGGTCAAATCCCCACCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGGCGTACCCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                                         <  1:554907/149‑1 (MQ=11)
                                                                                                                                        ccgcggaaCCAATTACATTTTGTTGAA‑TT‑CTCGGATTCCTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg            <  5:607190/145‑1 (MQ=14)
                                                                                                                                                      |                                                                                                                                                    
CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGG  >  NC_000913/2947317‑2947613

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