Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC34.5% 162.1 / 151.1 114intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (12/62);  new base C (1/38);  total (13/101)
Fisher's exact test for biased strand distribution p-value = 3.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-02

Only 100 of 158 total aligned reads displayed.

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGG  >  NC_000913/2947317‑2947612
                                                                                                                                                      |                                                                                                                                                  
ccAATTCGGTGAGCGGGTTAAAAAAAAGGTTGCATTAAAAGGCGAGGGGAGAATAGTGCGCTTCCACGGAGGCGGGGTGGAGCAGCCGGGTAGCTGGTCGGGGTCATACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                                   <  8:1583447/151‑1 (MQ=18)
ccAAATGCGTGATCGTGTTAAAAAAAAGGTTGCATGAATAAGCTAGCGGAGTATAGTTCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                    <  7:790098/151‑1 (MQ=18)
ccAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                    <  2:63742/151‑1 (MQ=255)
ccAAATCCGTGATCGGGGTAAAAAAAAGGGTGCATGAAAACTCTATCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                    <  3:72276/151‑1 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatt                                                                                                                                                  >  4:912231/1‑149 (MQ=255)
   aaTGCGTTATCGTGGTAAAAAAAAGGTTGCATGAAAAATATATCGCAGTATATTGCGCATACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                                 <  1:1170793/151‑2 (MQ=17)
   aaTCGGTGATCGGCTTAAAACCAAGCTTGAATGAAAACCAGAGCGGAGTATACTCCGCATCCACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                                 <  3:1515395/151‑2 (MQ=17)
   aaTCCGTTCTCGGGGTAAAAAAAAGGTTGCATGAAAACTATATCGGAGTATAGTGCGCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                                 <  5:928633/151‑4 (MQ=18)
   aaTCCGTTATCGGGTTAAAAAAAAGGTTGCATGAAAACTCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                                 <  6:1708642/151‑2 (MQ=255)
   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattt                                                                                                                                                 >  8:1565415/1‑148 (MQ=255)
    atgagTGATGGGGTTAAAAATAAGGGTGAATGAAAAAGAGAGAGTAGTATTGTGCGCTTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                                                                <  5:804244/147‑3 (MQ=14)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCAGGCCCCCGCAACCAATTAAAAt                                                                                                                                                <  7:349323/151‑1 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGTTAGCGCTGCGGGCTCATAACCCGAAGGGCGGCGGTGCAAATCCGCCCCCCCCAACCAATTattcaa                                                                                                                                               >  2:1018957/1‑145 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                                               <  8:532070/150‑1 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAtttt                                                                                                                                              >  8:1030325/1‑151 (MQ=255)
       ggTGATCGGGGTAAAAAAAAGGGTGCAGGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                                            <  2:257132/150‑1 (MQ=255)
        gTGATCGGGGTAACAACACGGTTGAATGAAAACTCTAGTGGAGTATAGTGCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                                                                                           <  8:630844/151‑1 (MQ=17)
         tGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacc                                                                                                                                           >  1:1311420/1‑142 (MQ=255)
         tGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttttctt                                                                                                                                           >  7:1069689/1‑148 (MQ=255)
          gctCGGGGTAAAACAAACGCGGCATGAAAACGCGAGCGGAGTATTGTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatga                                                                                                                                         <  6:638163/149‑1 (MQ=17)
          gatggtgttAAAAGAAAGGGTGCATGAAATATCTAGCGCAGTATAGACCGCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTTGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatg                                                                                                                                          <  3:1350336/142‑1 (MQ=11)
          gATCGTGGTAAAAAAAAGGTTGAATGAAAAATCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccc                                                                                                                                          <  7:1227513/151‑11 (MQ=18)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                                                                          <  2:1477155/151‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                           <  5:1734772/150‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                                            >  1:1289846/1‑148 (MQ=255)
          gATCGGGGAAAAAAAACCGCTGCATGAAAACGCGAGCGGAGTAGACTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                                                                                           <  3:527303/149‑1 (MQ=17)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                         <  6:1177512/151‑1 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                         <  8:208280/151‑1 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                         <  5:407015/151‑1 (MQ=255)
            gcGGGCTAAAAAAAACGTTGCATCAAAACGCCAGCGGAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAg                                                                                                                                       <  3:730785/150‑1 (MQ=17)
             ggtgttAAAAATAAGGTTGAATGAAAACGAGATTGGTGTATATTGCGATTCCTCGGACGCGGGGTGGAGCAGCCTGGTATCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt                                                                                                                                      <  5:590296/146‑1 (MQ=14)
             ggggtTAAAATAAAGGTTGCATGAAAACGCGAGCGTAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaa                                                                                                                                       <  5:853152/150‑14 (MQ=17)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                        <  8:1312642/150‑13 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATGCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                                                                       <  3:104647/151‑1 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgatcaccctac                                                                                                                                       >  8:1431684/1‑138 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaaccta                                                                                                                                        >  5:1225890/1‑138 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                                        <  1:391607/150‑1 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                                                                       <  8:560050/151‑1 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTAAAGTGCGCATCCACGGACGCGGTGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTAGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                                        <  1:1524366/150‑1 (MQ=255)
               ggtaaaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                                    <  2:915320/147‑1 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaacg                                                                                                                                     >  4:1499166/1‑136 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaaac                                                                                                                                     >  5:485793/1‑136 (MQ=255)
                gctAAAATACAGCTTGCATGAATAAGCGAGCGCAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                                  <  8:110579/149‑1 (MQ=14)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                                  <  4:252587/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                                  <  4:1088878/151‑1 (MQ=255)
                 gTAAAAAAAAGGGTGCATGAATACCCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                 <  6:1003628/151‑2 (MQ=17)
                  tAAAAAAAAGGTTGCATGAAAACGCTAGCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                 <  8:25693/150‑2 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                                                 <  5:28962/150‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGc                                                                                                                                <  3:219763/151‑1 (MQ=255)
                     aaaaaaGGGTGAATGAAAAACCCAGGGTAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagcag                                                                                                                             <  1:1567372/151‑6 (MQ=11)
                      aacaaGGGTGAATGAAATAGCCATCGGAGTAGAGTCCGACTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                                                                                                                             <  8:391640/148‑1 (MQ=11)
                      aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                                                             <  4:297831/151‑1 (MQ=255)
                      aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                                                             <  2:2217/151‑1 (MQ=255)
                      aaaaaGGTTGAAGGAAAAATAGAGCGGAGTATAGTCCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                                                                                                                             <  4:193631/151‑1 (MQ=11)
                       aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaa                                                                                                                                               >  3:567242/1‑128 (MQ=255)
                       aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                                                            <  4:565105/151‑1 (MQ=255)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                           <  1:1696085/151‑1 (MQ=255)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                           <  7:314795/151‑1 (MQ=255)
                        aaaGGTTGCATGAAAAAGAGAGCGGTGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                           <  1:1028983/151‑1 (MQ=14)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTata                                                                                                                                                  >  6:1345240/1‑126 (MQ=255)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                                          <  2:1590606/151‑1 (MQ=255)
                          cggatcaaTGAAAAAGAGATGGGAGTATTTTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                                         <  4:396081/144‑1 (MQ=1)
                           ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                                                                                        <  6:1639750/151‑1 (MQ=255)
                            ggtGCACGAAATCGCGATCGTAGTATAGTGGGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                                       <  2:1517835/149‑1 (MQ=11)
                              tGCATGAAAACTCTAGTGTTGTATTTTGCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                      <  6:575058/150‑1 (MQ=11)
                              tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                                      <  2:1179951/150‑1 (MQ=255)
                              tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                                      <  8:623164/150‑1 (MQ=255)
                              ggCAGGAAAACGCGTGGGGAGTATTGTTCGCCTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                                     <  2:1584724/150‑1 (MQ=11)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTACATCAACAAGCGCAGTACTGTGAcg                                                                                                                    >  3:103107/1‑151 (MQ=18)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAc                                                                                                                     <  8:1379276/150‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAc                                                                                                                     <  8:1081260/150‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                                    <  8:1291116/151‑1 (MQ=255)
                                 atgaaaaatcttgCGGAGTATAGTCCGACTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                                  <  2:1743734/139‑1 (MQ=11)
                                   gAAAACCCGTGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg                                                                                                                 <  6:256442/150‑1 (MQ=11)
                                    aaaaCTCGAGCGGAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg                                                                                                               <  4:751870/151‑1 (MQ=14)
                                     aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                                              <  4:1405813/151‑1 (MQ=255)
                                      taCGCTAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                    <  4:1406145/112‑1 (MQ=255)
                                       aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                        >  8:1471954/1‑112 (MQ=255)
                                       aCGCGAGCGGAGTATAGTGCCCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                                                                            <  2:1626874/151‑1 (MQ=255)
                                        ctcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                                                                                                                                  <  1:90771/111‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                                                                                                                                  >  8:90675/1‑113 (MQ=255)
                                         gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTATTAAATAAAAGCAGTACGGTGAcg                                                                                                                    >  7:1700159/1‑141 (MQ=17)
                                          cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacac                                                                                                                                            >  4:142604/1‑109 (MQ=255)
                                          cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacac                                                                                                                                            <  6:142340/117‑9 (MQ=255)
                                            aGCGTAGTATTCTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTa                                                                                                                                    <  1:4906/121‑1 (MQ=11)
                                            aGCGGAGTATAGTTCGCATCCACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCTGTCTCATAACCCGAAGGTCGTCGGTTCAAATACGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                                       <  6:202412/151‑1 (MQ=255)
                                            aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTGGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                                       <  3:968020/151‑1 (MQ=255)
                                               ggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                                                    <  2:199790/151‑1 (MQ=255)
                                                 aGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                                            <  5:1515982/109‑1 (MQ=255)
                                                 aGTAGAGTGCTCATCCACTGAATCTGTTTGGTTTAGCGTGTTAGCTTGTGGCTCTCATAACCCGAATCTCGGCGGTTGATATCCTGCCCCCGGAACCAATTAACCTTTGATGAAGTCAACCACTACCCTGACGCGGGGTGGAGCAGCCTg                                                                                                   <  1:695891/150‑1 (MQ=9)
                                                    attgTGCTCATCCACGGAAGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGATCATAAGCCGAATTGCGTCGGGTCAAATCAGTCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                                <  4:1156151/147‑1 (MQ=17)
                                                    ataGTGCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacc                                                                                                                                           <  7:532726/108‑10 (MQ=11)
                                                         ccggctCCTCGGCCGCGGGTTGGATCAGCGTTGTATCTCGGCTGTCTAATAACCCTAAGGGCGGCGCTTCAAAACCGGCCCCCGCAACCAATTAAATATTGCAGAATTAACCGCCTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCtc                                                                                            <  3:1537193/144‑1 (MQ=2)
                                                           gcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                    >  3:1224797/1‑92 (MQ=40)
                                                                                                           aaCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaatt                                                                                                                                                 >  7:821388/1‑45 (MQ=14)
                                                                                                                                 aTCCGGCCCCCGCAACCAATTAAAATTTTTTTAATTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCgg                                                     <  5:703211/116‑1 (MQ=255)
                                                                                                                                 aTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCgg                                                     >  2:705000/1‑116 (MQ=255)
                                                                                                                                              aaCGAATTATATTTTGTTGAAGTTTCGCATTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc     <  1:283879/151‑1 (MQ=14)
                                                                                                                                                 cAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGCTCGTCGCTTCAAATCCGGCCCCCGCAACCAATCAAATAAGAAGAACAAAAAGCAGTAATGTGACGAggg  >  7:390548/1‑151 (MQ=17)
                                                                                                                                                   aTTAAATTTTTATTACGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc                                                <  5:1615905/103‑1 (MQ=11)
                                                                                                                                                      |                                                                                                                                                  
CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGG  >  NC_000913/2947317‑2947612

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: