Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC25.2% 349.2 / 135.6 151intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (13/100);  new base C (0/38);  total (13/138)
Fisher's exact test for biased strand distribution p-value = 3.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

Only 100 of 204 total aligned reads displayed.

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTG  >  NC_000913/2947317‑2947605
                                                                                                                                                       |                                                                                                                                           
ccAAATCCGTGATCGGGGTATAAAAACGGTTGCATGAAAACTCGATCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                             <  3:335485/151‑1 (MQ=255)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACTCTAGGTGAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                            <  2:720228/151‑2 (MQ=255)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                            <  1:1598151/151‑1 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCGCGCAACCAATTatt                                                                                                                                           >  7:1178319/1‑149 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatt                                                                                                                                           >  1:216043/1‑149 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                                                           <  8:779468/151‑1 (MQ=255)
   aaTCGGTGATCTGGTTAAAAAAAAGGTTGCATGAAAACGATAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                          <  6:1791083/151‑4 (MQ=21)
   aaTCCGGGCGCGGGCTAAAAACCAGCTGGAATCAAAACCCGAGCGGAGTATACTCCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                          <  3:138273/151‑4 (MQ=14)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAATCGCGAGCTGAGTATAGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                         <  4:1694297/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                         <  1:1145914/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                         <  1:402909/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                         <  8:1710248/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGCTGCATGAAAACGCTAGCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                         <  8:1149921/151‑5 (MQ=255)
     tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                                        >  1:633584/1‑151 (MQ=255)
     tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                        >  3:856044/1‑151 (MQ=255)
      ccGTGATCGGGTTAAAAAACAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                        <  1:1028385/150‑1 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                       <  4:1357833/151‑1 (MQ=255)
       ggTGATCGGGGTAAAAAAAAGGGGCAATGAAAACGCTATCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                             <  5:1489656/143‑1 (MQ=18)
       cGTGATCTGGCTAAAACAAAGGTTGCAAGAAAACTATATCGTAGTATACTTCGCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                                     <  6:1179217/151‑1 (MQ=14)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                                                                       >  6:1532123/1‑144 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaaa                                                                                                                                       >  6:26383/1‑144 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaatgaac                                                                                                                                       >  6:580942/1‑145 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                                      <  6:1584033/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAAttt                                                                                                                                       >  6:1014738/1‑150 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAAGCTCGCAAGCGCGTCCAATTAATAtttt                                                                                                                                      >  7:808201/1‑150 (MQ=25)
        gTGATGGTGTTAAAAAAAAGGCTGCATGAAAAATCGATCGTAGTATAGTTCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                                                                                    <  6:1500909/151‑1 (MQ=14)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacac                                                                                                                                     >  1:1751050/1‑143 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                     <  3:272513/151‑1 (MQ=255)
         ttatgtggctAAAAATAAGGTGGCATGAATTAGATATCGTTGTATTGTCCGAATACACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                                                                                    <  7:338435/141‑1 (MQ=2)
         ggATCGGCGTAAAAACAAGCCTGCATCAAAACCCGATCGCAGTATACTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatg                                                                                                                                   <  8:35756/150‑1 (MQ=14)
          gATCGTGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                                                                   <  7:850530/151‑1 (MQ=255)
          gATCGTGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                    <  1:926882/150‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                                                                   <  6:831541/151‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                                                                   <  4:172231/151‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                    <  8:1353384/150‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                    <  8:572065/150‑1 (MQ=255)
           gagcgtctaaaaaaaacCGCTGCAGGAATACCCAATCGCAGAATACTCCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatga                                                                                                                                  <  3:1658966/142‑1 (MQ=2)
           atgttgttAAAATTAAGGTTGCATGAAATATATATCGTTGTATTGTTCGATTCCTCTGTCGCGGGGTGGATCAGGCTTGTAGCTCGTCGGGCTCATAACGCGAAGGTCGTCTGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                                                                                                                 <  6:1673071/144‑1 (MQ=2)
           atggtgttAAAAAACAGGTTGCATGAAAACGCGAGCGCAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatga                                                                                                                                  <  3:1027691/143‑1 (MQ=18)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                  <  6:662580/151‑1 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTAAAATttgttg                                                                                                                                   >  1:994132/1‑150 (MQ=255)
            gcGGGGTAAAAAAAAGGTTGCATGAAAACTCGAGCGGAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAg                                                                                                                                <  2:1231721/150‑1 (MQ=255)
             gggggTAAAAAAACCGGTGCAGGAAAACGCCATCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt                                                                                                                               <  1:615931/150‑1 (MQ=17)
             gggcctAAAAAACAGGTTGAATGAAAACACGAGCGTAGTATACTTCGACTCCACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt                                                                                                                               <  8:117938/146‑1 (MQ=14)
             cGGGGTAAAAAAAAAGGTTGAATGAAAAAGCGAGCGGAGTATAGTGCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                                                                                  <  4:1064547/150‑12 (MQ=17)
             cGGGGTAAAAATAAGGGTGCATGAATACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                 <  4:1211167/150‑13 (MQ=17)
              ggggTAAAAAAAAGGTGGCATGAAAACGCGAGCGGTGTATAGTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaac                                                                                                                               <  8:846124/151‑15 (MQ=17)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                             <  5:1797785/151‑1 (MQ=255)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                                            <  3:1074535/151‑1 (MQ=255)
                ggTAAAAAAAAGGGTGAATGAAAAATCGATCGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                           <  7:359268/151‑1 (MQ=17)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                           <  4:1587658/151‑1 (MQ=255)
                 gTAAAAAAAAGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                           <  7:1541723/151‑1 (MQ=255)
                  ttataattaaGGTTGAATGAAATAGCTTGAGTAGAATATTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                           <  4:311379/142‑1 (MQ=11)
                    aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                                          <  1:1318009/150‑1 (MQ=255)
                    aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGc                                                                                                                         <  5:678433/151‑1 (MQ=255)
                    aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGc                                                                                                                         <  7:1362646/151‑1 (MQ=255)
                     aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                                        <  8:542579/151‑1 (MQ=255)
                     aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAATAAGTaaaagc                                                                                                                        >  2:1088459/1‑147 (MQ=255)
                     aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                                        <  4:451999/151‑1 (MQ=255)
                      aaaataGGTTGCATGAATACGCGAGGTGAGTATAGTGCGGATCCACGGACGCGGGGTGGAGCAGCCTGTTAGCTCGTCGGGCTCATAACGGGAAGGTCGTCGGTTCATATCCGGACCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                                       <  2:1636296/151‑1 (MQ=18)
                      aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGCCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                                       <  2:1021697/151‑1 (MQ=255)
                      aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaac                                                                                                                               >  5:313604/1‑130 (MQ=255)
                        aaaaGGTTGCATGAAAACGCGAGCGCAGTATAGTGCGACTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                                                                                                                     <  8:1200827/151‑1 (MQ=14)
                         aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                    <  1:357147/151‑1 (MQ=255)
                         aaaGGATGAATGAAAAATCGAGAGGAGTATAGAGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                    <  8:1215501/151‑1 (MQ=11)
                          aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                                   <  5:815082/151‑1 (MQ=255)
                            ggctGCATGAAATCGCGAGCGCAGTATAGAGGGATTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                                 <  6:1413467/148‑1 (MQ=1)
                            ggTTGCATGATAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATTACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                                                                                 <  5:1091545/151‑1 (MQ=255)
                            ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGTGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTGAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                                                                                 <  7:921063/151‑1 (MQ=255)
                             gTTCCATGAAAAAGCGAGCGGAGTATACTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                         <  5:1399497/127‑5 (MQ=17)
                               tgaaTGAAAAATAGATCGTAGTATAGTGCGAATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                               <  6:1638913/147‑1 (MQ=11)
                               tGCATGAAAACTATATCGTAGTATATTCCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                               <  4:223158/150‑1 (MQ=11)
                                gaaggaagtcgagtgCGTAGTATTGTGGGCTTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                             <  5:1138516/137‑1 (MQ=11)
                                gCATGAAAACTATAGCGTAGTATAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                             <  4:601378/151‑1 (MQ=11)
                                 cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                                                                       >  4:536843/1‑119 (MQ=255)
                                 cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                          >  2:934877/1‑119 (MQ=255)
                                  accgatcttagAGCGGAGTATTGTGCGACTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                           <  8:1500623/141‑1 (MQ=11)
                                    gAAATCGCGAGCGGAGTATAGGGCGCATCGACGGACGGGGGGGTGAGCAGGCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATGCGGCCCCCGCAACCAATTAAAATTTGATGAAGTATAGCAGTACGGTGACGCggg                                                                                                         <  8:892124/151‑1 (MQ=18)
                                    gAAAACGCGTTGGTAGTATTCTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCggg                                                                                                         <  3:759568/151‑1 (MQ=11)
                                    gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTTGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCggg                                                                                                         <  8:1380025/151‑1 (MQ=255)
                                    gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCgg                                                                                                          <  3:232605/150‑1 (MQ=255)
                                    gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCgg                                                                                                          <  6:1384594/150‑1 (MQ=255)
                                      aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                                       <  3:447654/151‑1 (MQ=255)
                                       aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattaaacac                                                                                                                                     >  1:950312/1‑113 (MQ=255)
                                       aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattaaacac                                                                                                                                     <  2:948561/121‑9 (MQ=255)
                                       aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                                                                      <  3:184405/151‑1 (MQ=255)
                                        aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                                                                     <  6:685440/151‑1 (MQ=255)
                                         cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                    <  5:1328260/151‑1 (MQ=255)
                                          gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGag                                                                                                   <  5:1728651/151‑1 (MQ=255)
                                          gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGag                                                                                                   <  4:351698/151‑1 (MQ=255)
                                            tagagtaGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                          <  2:747829/105‑2 (MQ=255)
                                            gAGCGGAGTATAGTGCGCATCCTCTGACGCGGGGTGTTGGAGCGTGGTAGCTGGTCGCGCTCATATCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCTCAACGAATTAAAATTTGATGAAGTACAGCAGTACGCTGACGCGGGGTGGagca                                                                                                 <  2:337645/151‑1 (MQ=17)
                                            gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                                                 <  5:473110/151‑1 (MQ=255)
                                             cgCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTAATAACCCGAAGGTCTTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCATTACGGTGACGCGGGGTGGagcag                                                                                                <  6:447107/150‑1 (MQ=255)
                                             aGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacc                                                                                                                                    <  3:1539435/116‑10 (MQ=255)
                                                  aGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGCGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGCTGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                           <  3:526951/151‑1 (MQ=255)
                                                   gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                          <  5:396704/151‑1 (MQ=255)
                                                     ataGTCCGCATACCCTGACGCGGGGTGGAGTAGGCTTGTAGGTCGTGGGTCTCATAACCCTAAGGTCGTCGGTTGTAATCCGTCCCCCGCAACCAATTAAACTTTGCTGAAGTAACGCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTa                                                                                         <  7:1287988/150‑1 (MQ=14)
                                                                                                ctgcgGCCTCATAACCCGAAGCTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAACTTTGATGAAGTAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATcc                                                <  7:504993/145‑1 (MQ=255)
                                                                                                                                           cctcAACGAATTACACTTTGATGAATTCTCGGAATACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg  <  5:1016694/148‑1 (MQ=14)
                                                                                                                                                       |                                                                                                                                           
CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTG  >  NC_000913/2947317‑2947605

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: