Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA43.0% 144.4 / 52.7 107E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (16/45);  new base A (46/0);  total (62/45)
Fisher's exact test for biased strand distribution p-value = 7.96e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508519‑3508810
                                                                                                                                                 |                                                                                                                                                     
gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTcacac‑a                                                                                                                                               <  1:2086166/151‑1 (MQ=255)
 tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATcaca                                                                                                                                                  >  2:345067/1‑148 (MQ=255)
   tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGCGGGAAAAAttaaaca                                                                                                                                                 >  1:1189530/1‑143 (MQ=255)
   tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACg                                                                                                                                             <  1:1132175/150‑1 (MQ=255)
      aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACGGTa                                                                                                                                          >  1:1202713/1‑150 (MQ=255)
      aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTa                                                                                                                                          <  1:903855/150‑1 (MQ=255)
       aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGc                                                                                                                                        <  1:1334215/151‑1 (MQ=255)
        aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGcc                                                                                                                                       <  1:987869/151‑1 (MQ=255)
        aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGCGTAACCCGTCGGATCAAAACTCATAGGCCGTCTCCTGCAAACAGGAGGGAGGACAGAAAAGAGCTCCCCTTTTGTGGC‑GGCCGCGCGC‑GGAGAAAAACATAC‑ACGCTGGcc                                                                                                                                       >  2:1119041/1‑151 (MQ=255)
           cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGGCG‑GC‑GGGAAAAAAcacac‑aacggaggcccc                                                                                                                                    >  2:2125252/1‑140 (MQ=255)
           cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAATTTACAA‑AAAGGAGccccc                                                                                                                                    >  1:774560/1‑151 (MQ=255)
           cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAAACACAC‑ACGGTAggcgcc                                                                                                                                    >  2:602288/1‑146 (MQ=255)
           cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACGGTAGgcccc                                                                                                                                    >  1:797131/1‑151 (MQ=255)
            ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCca                                                                                                                                   <  1:1985784/151‑1 (MQ=255)
             gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACGGTAGccccccc                                                                                                                                  >  1:250516/1‑149 (MQ=255)
              agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGG‑GGGAAAAATCAAACAACGGAAccccccc                                                                                                                                   >  2:1560119/1‑149 (MQ=255)
               gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATCACAC‑ACCGGAGccccccc                                                                                                                                  >  1:708503/1‑147 (MQ=255)
               gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCcaca                                                                                                                                 <  2:675059/150‑1 (MQ=255)
                agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAAATCACAA‑AACGGAGccccc                                                                                                                                    >  1:1962140/1‑146 (MQ=255)
                agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGAGGCGGGG‑GGAAAAAATAACAC‑AAGGTAGCCCCcccacc                                                                                                                               >  1:1966717/1‑149 (MQ=255)
                  aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATCAAAA‑ACGGTAGccccccc                                                                                                                                  >  2:733646/1‑144 (MQ=255)
                  aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAATTCACAA‑ACGGGAGGCCCCCCATc                                                                                                                               >  2:2053604/1‑149 (MQ=255)
                  aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCa                                                                                                                              <  2:726328/150‑1 (MQ=255)
                  aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAg                                                                                                                             <  1:1494043/151‑1 (MQ=255)
                     aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGGGG‑GGGAAAAAAAACAC‑AAAGAAGACCCCACCTCa                                                                                                                              >  2:215201/1‑147 (MQ=255)
                     aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGaaa                                                                                                                          <  1:1353447/151‑1 (MQ=255)
                       tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATCAAAC‑ACCGGAGCCCCCCCACCAGaaa                                                                                                                          >  2:1986561/1‑149 (MQ=255)
                        atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGt                                                                                                                                           <  1:1360129/131‑1 (MQ=255)
                          gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                     <  2:202338/151‑1 (MQ=255)
                            aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATTACAA‑ACGGGAACCCCCCCAACAGgaaagggggg                                                                                                                   >  2:622792/1‑144 (MQ=255)
                            aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGGGGGGCGC‑GGGAAAAATCACAC‑ACGGTAGCCCCCACATCaacaagagggg                                                                                                                    >  2:1679631/1‑140 (MQ=255)
                             gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCgtggtg                                                                                                                  <  2:1442762/151‑1 (MQ=255)
                              aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGGGGCGCGC‑GGGAAAAATCAAAC‑ACCGTAGCCCCCACCCCAGAAAgggggg                                                                                                                    >  2:894283/1‑143 (MQ=255)
                                 cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATCACAC‑ACCGTAGCCCCCCCATCAGAAAACGTGGTgcg                                                                                                                >  1:195472/1‑149 (MQ=255)
                                 cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATAACAC‑ACGGTAGCCCCCCCCTCAGAAAGCgtgggggg                                                                                                                >  2:1167630/1‑145 (MQ=255)
                                 cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATAACAC‑ACCGACGACCCCACATCACaaaaagcgggggg                                                                                                                >  1:2118889/1‑139 (MQ=255)
                                  tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATTACAA‑ACAGGAGACCCCCCACCCGAGAGCGTGGGgcg                                                                                                                >  2:1348757/1‑148 (MQ=255)
                                   cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGG‑GGAGAAATTCACAC‑ACGGGAACCCCCCCATCAGAAAACGGGGGGCGgac                                                                                                             >  2:344083/1‑150 (MQ=255)
                                   cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAATTCACCC‑ACCGTAGCCCCCACCTCCGAAAGCGTGGGgcgc                                                                                                               >  1:1335024/1‑148 (MQ=255)
                                      gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACAcc                                                                                                          <  1:1738662/150‑1 (MQ=255)
                                        gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                                       <  1:1185140/151‑1 (MQ=255)
                                          gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACCGTAGCCCCCCCATCAGAAAACGGGGGGCGGACAcc                                                                                                          >  1:1780406/1‑146 (MQ=255)
                                          gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑AAGGGAGCCCCCACAACAGAAAACGGGGGGCGCACAcc                                                                                                          >  1:883915/1‑146 (MQ=255)
                                          gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                     <  1:567692/151‑1 (MQ=255)
                                          gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                     <  1:703801/151‑1 (MQ=255)
                                          gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                     <  2:759331/151‑1 (MQ=255)
                                          gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                     <  1:1415157/151‑1 (MQ=255)
                                           aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                     <  2:1041094/150‑1 (MQ=255)
                                              gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                                                                 <  2:2009836/151‑1 (MQ=255)
                                                gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAATTTAAAC‑ACCGGAGACCCCACATCAGAAAACGGGGGGGGCACACCCAGccg                                                                                                    >  2:1377621/1‑144 (MQ=255)
                                                  aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                              <  1:987367/150‑1 (MQ=255)
                                                  aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAAATCACAC‑ACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAg                                                                                               >  1:230421/1‑149 (MQ=255)
                                                      cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                                                                         <  2:444649/151‑1 (MQ=255)
                                                         cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGAAAAAATCAAAC‑ACCGGAGCCCCCCCATCAGAAAACGGGGGGCGGACACCCAGCTGCGTAGCGatca                                                                                         >  1:1422883/1‑148 (MQ=255)
                                                          ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGAGGGGCGC‑GGGAAAAATTAAAC‑ACGGGAGCCCCCAAATTCCAAAGCGGGGGGGGGACACAGGGGCGCGTAGCGATCATccc                                                                                     >  2:1688317/1‑151 (MQ=255)
                                                           gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                                                                    <  2:1602718/151‑1 (MQ=255)
                                                               aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGAAAAAATCACAC‑ACGGGAGCCCCCACATCAGAAAAAGTGGGGCGCACACCCAGCCGCGGAGCGATCATCCCtt                                                                                   >  1:898001/1‑148 (MQ=255)
                                                                aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                                >  1:1633387/1‑150 (MQ=255)
                                                                  aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAAATCACAC‑ACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTggt                                                                             >  1:1288989/1‑149 (MQ=255)
                                                                   cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGAAAAAATCACAC‑ACGGTAGCCCCCCCATCCGAAAGCGGGGTGCGCACACCGAGCTGCGTAGCGACCCTCCCTTGCTGAt                                                                             >  2:647477/1‑150 (MQ=255)
                                                                    tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATTACAA‑ACCGGAGCCCCCACATCAGAAAGCGGGGGGCGCACACCGAGCTGCGGAGCGATCATCCCt                                                                                    >  2:1635747/1‑142 (MQ=255)
                                                                       tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACGGGAGCCCCCCCAACAGAAAACGGGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                          >  1:2033773/1‑149 (MQ=255)
                                                                         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAg                                                                      <  2:658321/151‑1 (MQ=255)
                                                                             gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACGGTAGCCCCCCCATCAGAAAGCGTGGTGGGCACACCGAGCTGCGTAGCGATCATCCCTTTCTGGTAAAACAGCAc                                                                   >  1:125593/1‑150 (MQ=255)
                                                                                 ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                                 <  1:1806679/148‑1 (MQ=255)
                                                                                 ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGGGGGGCGG‑GGGAAAAATCACAA‑ACGGGAGCCCCCCCCTCAGAAAGGGGGGGGCGCACAACGAGCCGGGGGGCGAACATCCCCTGCTGGTAAAACAGCACTgg                                                                >  1:1651518/1‑149 (MQ=255)
                                                                                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac                                                             <  1:1634353/149‑1 (MQ=255)
                                                                                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAAATCACAC‑ACGGGAGCCCCCACAACAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGcaccac                                                          >  1:801419/1‑151 (MQ=255)
                                                                                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCAAAA‑AACGGAGCCCCCACACCAGAAAAAGCGGGGTGCCCACCGAGCTGCGTAGCGATCATTCCCTGCTGGTAAACCAACACTGGcacca                                                           >  1:975461/1‑149 (MQ=255)
                                                                                         ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATCACAC‑ACGGTAGCCCCCACATCAGAAAGAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAACACTGGcacc                                                            >  1:2138640/1‑145 (MQ=255)
                                                                                          cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                                     <  2:353050/151‑1 (MQ=255)
                                                                                          cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                                     <  2:1711944/151‑1 (MQ=255)
                                                                                           ggagggagGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGGAAAAATTAAAC‑ACCGTAACCCCCCCATCAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGGTAAAACAGCACTGGCCCCACCTcgc                                                     >  2:1864474/1‑150 (MQ=255)
                                                                                             agggaggAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGCGC‑GGAAAAAATCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCggg                                                  >  1:304519/1‑151 (MQ=255)
                                                                                              gggaggAAAGAAAAGAGCGGCCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                 <  2:1497429/151‑1 (MQ=255)
                                                                                              gggaggAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATCACAC‑ACCGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCCGCGTAGCGGTCATCCCTTGCTGGTAAAACAGCACTGGCACCACCTCGCgggg                                                 >  1:2044895/1‑150 (MQ=255)
                                                                                               ggaggaAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGGGGGC‑GGGAAAAATTACAC‑ACGGGAGCCCCCCCAACAGAAAACGTGGGGCGCACACCGAGCTGCGGAGCGGTCAACCCCTGCTGGTAAAACAGCACTGGCACCACTTCCCggg                                                  >  1:650256/1‑149 (MQ=255)
                                                                                                   gaaagaaaAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                            <  2:1965609/151‑1 (MQ=255)
                                                                                                         aaagaGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAATTCACAC‑ACGGGAGCCCCCCCCTCCGAAAGCGTGGGGCGCACACCCAGCAGCGTAGCGATCATCCCTTGCTGATAAACCAACACTGGCACCACCTCGCGGGGGTGGGGGcgg                                       >  2:806239/1‑148 (MQ=255)
                                                                                                          aagaGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                     <  1:1722061/151‑1 (MQ=255)
                                                                                                          aagaGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCGCGCTGCGTAGCGATCATCCCTTTCTGATAAACCAGCACTGGCACCACCTCGCGGGGATGGTGGgtgt                                      >  1:722729/1‑150 (MQ=255)
                                                                                                               cGACCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAAATTACAC‑ACGGTAGCCCCCCCATCAGAAAAAGTGGTGCGCACACCCAGCCGCGTAGCGATCATCCCCTTCTGATAAAACCGCACCGGCACCACTTCCCGGGGATTGTGGGTGTGAcc                                  >  1:1868946/1‑149 (MQ=255)
                                                                                                                gACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGGAGCCCCCCCATCCGAAAGCGGGGGGGGCACACCCAGCTGCATAGCGATAATCCCTTGCTGATAAACCAACACTGGCACCACCTCGCGGGGATGGTGGCtgtg                                     >  2:192898/1‑145 (MQ=255)
                                                                                                                 aCCATTTTGTGGC‑GGTGGCGCGC‑GGGAAAATTTACAC‑ACGGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGAcc                                  >  1:2080924/1‑147 (MQ=255)
                                                                                                                   cATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                            <  1:438898/151‑1 (MQ=255)
                                                                                                                   cATTTTGTGGC‑GGTGGCGCGC‑GGAAAAAATCACAC‑CCGGTAGCCCCCACACCAGAAAGCGTGGGGCGCCCACCGCGCTGCCTAGCGATCATCCCTTGCCGATAAACCAGCACTGGCACCCCTTCGCGGGGATGGTGGccgggatc                                  >  2:771647/1‑138 (MQ=255)
                                                                                                                      tttGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                          <  1:756524/150‑1 (MQ=255)
                                                                                                                         gtggc‑ggtggcgCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                      <  1:940207/151‑1 (MQ=255)
                                                                                                                           ggc‑ggtggcgCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCTCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                     <  1:968548/150‑1 (MQ=255)
                                                                                                                            gc‑ggtggcgCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                   >  1:2112819/1‑151 (MQ=255)
                                                                                                                             c‑ggtggccCGC‑GGAACAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                  >  2:1501020/1‑151 (MQ=255)
                                                                                                                               ggtggcgCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc                 <  2:883794/151‑1 (MQ=255)
                                                                                                                                  ggcgCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg               <  2:1213821/150‑1 (MQ=255)
                                                                                                                                    cgcgc‑gGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGa            <  1:1907631/151‑1 (MQ=255)
                                                                                                                                     gcgc‑gGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt           <  1:2143356/151‑1 (MQ=255)
                                                                                                                                       gc‑gGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         >  2:1159222/1‑151 (MQ=255)
                                                                                                                                        c‑gGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg        <  2:15171/151‑1 (MQ=255)
                                                                                                                                          ggAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       >  1:1697412/1‑151 (MQ=255)
                                                                                                                                          ggAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCACCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       >  1:1687887/1‑151 (MQ=255)
                                                                                                                                           gAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc      >  1:1901371/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGC‑GGTGGCGCGC‑GGAAAAATTCACAC‑ACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508519‑3508810

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