Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC34.0% 150.7 / 137.2 108intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (10/60);  new base C (0/36);  total (11/97)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-02

Only 100 of 153 total aligned reads displayed.

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGG  >  NC_000913/2947317‑2947613
                                                                                                                                                      |                                                                                                                                                   
ccAAATCCGTGATCTGGTTAAAAACAAGGGTGCATGAAAAATCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                     <  2:1431210/151‑1 (MQ=21)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                                    <  8:1197098/151‑1 (MQ=255)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                                    <  5:259364/151‑1 (MQ=255)
 cAAATCCGTGATCGGCTTAAAAATAAGGGTGAAAGAAAACGCGAGCGGAGAATACTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                                    <  5:445234/151‑2 (MQ=18)
  aaaTGCGTGAGCGGGTTAAAATTCCGGTTGCATGAAAAATGGATGGTTTTATATTTCGCATACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAATGTCTTCGGTTCAAATCCGGCCACCGCAACCAATCaaa                                                                                                                                                   <  4:1563105/151‑1 (MQ=2)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                                                                   <  1:1399428/151‑1 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaa                                                                                                                                                    >  6:1170959/1‑150 (MQ=255)
   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaaa                                                                                                                                                  <  4:142119/151‑1 (MQ=255)
   aaTCCGTGAGCGGGGTAAAAAACCGGGGGCACGAAAAAGCGAGCGGAGTATACTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                                  <  7:354170/151‑4 (MQ=18)
    atgcggtttgggggTAAAAAAAAGGTTGCATGAAAAACCGTTGGGTGTATATTGCGCATCCAGGGACGCGGGGTGGATGATGCTGGTACGTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                                  <  6:203695/140‑4 (MQ=17)
    atcgctGATCGTGTTAAAAAAAAGGTTGCATGAAAAATAGAGCGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                                                                 <  2:1612035/146‑3 (MQ=18)
    aTCCGTGATCTGGTTATAAAAAAGGTTGAATGAAAACTCTAGGGGTGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                                 <  1:1326125/151‑5 (MQ=17)
    aTCCGTGATCGGGGTAACACAAAGGGGGCATGAAAACGCGATCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                                                                 <  4:1172513/151‑3 (MQ=18)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgt                                                                                                                                                 >  4:204063/1‑147 (MQ=255)
     tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAttt                                                                                                                                                >  6:763317/1‑151 (MQ=255)
      cggTGATCTGGTTAAAAACAAGGGTGCATGAAAACGATAGAGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                                <  6:350749/148‑1 (MQ=17)
      cgctGATCGGGGAAAAAAACAGGGGGCATGAAAACGCGAGGGTAGTATAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                                <  5:14730/147‑1 (MQ=17)
      ccgttatcttgttAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATATTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtttg                                                                                                                                               <  2:502319/139‑2 (MQ=18)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                               <  6:420178/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                                                                               >  1:1603315/1‑144 (MQ=255)
       cGTGATCGGGGGAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                                              <  3:1673456/151‑1 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                                             <  4:595854/150‑1 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGTGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTatttaaaa                                                                                                                                              >  5:1636368/1‑143 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                             <  6:664711/151‑1 (MQ=255)
        gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAAttga                                                                                                                                              >  6:595104/1‑148 (MQ=255)
         ggctggtgttAAAAAAAAGGTTGCATGAAAACGCGTGGGGTGTATTGTGCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCATTCAAATTtgatg                                                                                                                                           <  8:204834/142‑1 (MQ=17)
          gagcGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                            <  6:846051/147‑1 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGTAGCAGCCTGGTAGCTCGTCGGGGTCATATCGCGAAGGTCGTCGGTTCAAATCCGGGCCCCGCAACCAATTAAAATTtgtt                                                                                                                                            <  5:1554347/150‑3 (MQ=25)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgatcaccc                                                                                                                                           >  3:1167438/1‑141 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccc                                                                                                                                           >  4:1258614/1‑141 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                                                                                          >  4:1032428/1‑140 (MQ=255)
            tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                                                                                          >  4:119753/1‑139 (MQ=255)
            gcGGGGTAAAAAAAAGGGTGAACGAAAACGCGAGCGGAGTATAGAGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAg                                                                                                                                        <  5:1328155/150‑1 (MQ=18)
             tgtgttAAAATTAAGGTTGAATGAAATCGAGAGAGTAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAg                                                                                                                                        <  1:156318/145‑1 (MQ=14)
             ggggtTAAAAATAAGGTTGAATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaa                                                                                                                                        <  2:1090317/150‑14 (MQ=17)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                                                                        <  7:1612615/151‑1 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAATTatagaagacccta                                                                                                                                         >  6:1580044/1‑138 (MQ=255)
              ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaac                                                                                                                                       >  3:1480330/1‑137 (MQ=255)
              ggcgTAAAAAAACGGGTGCATGAAAACTCCTGGGTAGTATTGTCGCCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTa                                                                                                                                     <  5:78560/148‑1 (MQ=11)
               gggTAACAAAACCCTTGCATGAAAAACCGAGCGGAGTATAGTCCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaacg                                                                                                                                      <  3:680249/151‑16 (MQ=14)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                                    >  4:762383/1‑136 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                                     <  6:1080575/151‑1 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                                     <  1:715204/151‑1 (MQ=255)
               gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAAGAAGTa                                                                                                                                     >  2:1024140/1‑150 (MQ=255)
               gggTAAAAAAAAGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                                     <  4:1365016/151‑1 (MQ=255)
               ggctAAAAACAAGCGTCCATGAAAACTACATCGTAGTATATTCCCCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                                                                                                                                    <  4:875356/148‑1 (MQ=1)
                ggTAAAACAAACGGTGCAGGAAAACCCGTGCGGAGAATAGTCCCCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                                   <  5:1151998/151‑1 (MQ=11)
                ggTAAAAAAACGGCTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                                   <  7:1538602/151‑1 (MQ=255)
                 gTAACAACAAGGCTGAATGAAAACAAGATCGGAGTATACTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                  <  3:1456643/151‑2 (MQ=14)
                 gTAAACAAAAGGCTGCATCAAATCACGAGCGGAGTATACTACGAATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                  <  5:82883/151‑2 (MQ=14)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                                   <  8:1664086/151‑1 (MQ=255)
                 gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                                   <  5:1166939/151‑1 (MQ=255)
                  tAAAAATACGGTGGCATGAAAAAGCGAGCGGAGTATAGTGCGCCTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                  <  3:598397/150‑2 (MQ=17)
                  tAAAAAAACCGTTGCATGAAAACTCTAGCGGAGTATTGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                  <  1:466772/150‑2 (MQ=17)
                   aaaataaaGGGTGCATGAAAAACCGTTCGGAGTATACTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaag                                                                                                                                 <  4:764116/150‑3 (MQ=14)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                                                  <  4:542221/150‑1 (MQ=255)
                   aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGc                                                                                                                                 <  4:393722/151‑1 (MQ=255)
                    aaatcaaGGTGGCATGAAAAATCCAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagca                                                                                                                               <  6:1018834/146‑5 (MQ=14)
                    aaaaaacGGCTGAATGAAAACTCTAGCGGAGTATATTCGGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagca                                                                                                                               <  8:1570572/151‑5 (MQ=11)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtttgcagaaagaaaagc                                                                                                                                  >  2:761173/1‑135 (MQ=18)
                      atacaccttcaaTGAAAACGAGATGGCAGTATACTGCGCCTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                                                                                                                              <  3:761932/140‑1 (MQ=11)
                      aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                                                              <  6:60150/151‑1 (MQ=255)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                            <  5:9444/151‑1 (MQ=255)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                            <  3:1617109/151‑1 (MQ=255)
                        aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                            <  8:1480358/151‑1 (MQ=255)
                        aaaGGTTGAATGAAAACTCGATCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                            <  7:625152/151‑1 (MQ=14)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                                           <  2:942551/151‑1 (MQ=255)
                         aaGGGTGCATGAATAAGCGAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                                           <  8:1039296/151‑1 (MQ=14)
                          cggTTGCATGAAAACGCTTGTGTTGTATATTCGGATTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGGTCGGCGGGATCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                                          <  1:580622/150‑1 (MQ=11)
                            gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTg                                                                                                                        <  3:387958/151‑1 (MQ=255)
                             gtGCATCAAAAAGAGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                       <  4:1658782/150‑1 (MQ=14)
                              ggCATGAAAACGAGTGGGTAGTATTTTGCGATTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                                      <  2:1701989/150‑1 (MQ=11)
                              gcaatgaattctctaGCGGTGTATACTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                       <  3:525295/136‑1 (MQ=11)
                               gaaTGAAAACGCAAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                                     <  7:110221/149‑1 (MQ=14)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcg                                                                                                                     <  7:870164/151‑1 (MQ=255)
                               caaggAAAACTCCAGCGTAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                                      <  1:1411224/146‑1 (MQ=11)
                                 aTGAAAACTCTAGCGTATTATAGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                                         <  6:1261633/130‑13 (MQ=11)
                                 aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAGATCCGGACCCCGCACCCAATTAATAATTGATGAATTAAAGCAGTACGCTGCcgcg                                                                                                                   >  7:1387942/1‑151 (MQ=18)
                                 aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaaccta                                                                                                                                         >  4:1263093/1‑118 (MQ=255)
                                    tcctatatatcGGAGTATAGTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc                                                                                                                    <  4:936632/137‑1 (MQ=11)
                                    aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCgggg                                                                                                                <  2:138498/151‑1 (MQ=255)
                                       aCGCGAGCGGATTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                                                                             <  1:393146/151‑1 (MQ=255)
                                        cgggAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGGAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  2:659824/148‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  4:495595/151‑1 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  3:842731/151‑1 (MQ=255)
                                          cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacac                                                                                                                                             >  2:819629/1‑109 (MQ=255)
                                            aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAATAAATAAAAGCAGTACGGTGAcgc                                                                                                                    >  6:935324/1‑139 (MQ=17)
                                              cGGACTATAGTGCTCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGATCATAACCCGAAGGTCGGCGGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGAAGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGcc                                                                                                      <  1:1276032/151‑1 (MQ=17)
                                               ggAGTATAGTTATCATCTCCGGACTCTGTGTGGAGCAGTCTTGTATCTCGTTTGTCTCTTAACCCGAATGTGGTCTGTTCTTATAATTCAGCCTCAACGAATTAAAATTTGATGAAGTCAAGGAGTACTGTGACGCGGGGTGGAGCAGCCt                                                                                                     <  6:171843/151‑1 (MQ=9)
                                               ggAGTATAGTGCTCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                                                     <  6:642141/151‑1 (MQ=255)
                                               ggAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                                         <  8:622343/131‑1 (MQ=11)
                                                    ataGTGCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                                 <  4:1592405/103‑5 (MQ=255)
                                                    ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                                 >  6:1591554/1‑99 (MQ=255)
                                                          cgattcgcctgACGCGGTGTGGATCAGCCTGGTATCTCGGCGGCCTAATAACCCGAATGGCGGCGGTGCTAATCATTCCCCCGCAACCAATTACATTTTGATGAAGTAACGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                                          <  4:321431/141‑1 (MQ=2)
                                                                                  cagcCTGGTAGCTCGTCGCTCTCATAACCCGAAGGTCGTCGGTGCAAATCCCGCCCCCGCAACCAAATAAAAATTGATGAAGTATCGCACTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGtcgtcg                                                                  <  3:1225276/151‑1 (MQ=11)
                                                                                                                gAAGCTCGTCTGTTCAAATACGGCCCCCTCAACCAATTAAAAATGGATGAAGTAAAGCACTCCCCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                    <  7:1108262/151‑1 (MQ=14)
                                                                                                                     gcggcgcgtccaaTCCCCCCCCCGCAACCAATTATACTTTGATGTAGAAAACCCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                               <  3:951730/140‑1 (MQ=14)
                                                                                                                                           ggccagcAATTATAATTTGATGAAGTTTCGGCATACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa         <  2:66793/145‑1 (MQ=17)
                                                                                                                                              aaCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTTAGATGAAGGAAAAGccgaggggtgcgg       >  2:95197/1‑138 (MQ=14)
                                                                                                                                                  aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAATTTTCAGGAAGTAAAGGCTGTACTGTGACGCgggg  >  3:89573/1‑151 (MQ=17)
                                                                                                                                                      |                                                                                                                                                   
CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGG  >  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: