Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2970GC42.5% 120.5 / 53.3 80intergenic (+103/+126)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (5/41);  new base C (34/0);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.31e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGC  >  NC_000913/3423150‑3423447
                                                                                                                                                   |                                                                                                                                                          
aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCgg                                                                                                                                                           <  2:1190548/149‑1 (MQ=255)
              gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCCCCCGGCCGGTGCTTCGC‑aa                                                                                                                                           >  2:336906/1‑150 (MQ=37)
              gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGGTGCTTCGC‑aa                                                                                                                                           >  2:1149705/1‑150 (MQ=255)
              gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCGGCCGGTGCCTCGC‑aa                                                                                                                                           >  1:699431/1‑150 (MQ=255)
                      ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCCC‑AAAG‑TTC‑aa                                                                                                                                  >  2:355181/1‑149 (MQ=255)
                      ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑a                                                                                                                                   <  2:803520/148‑1 (MQ=255)
                        ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCCCCCGCCCCCCCCCCGGCCCGTGCGTCCC‑AACG‑TTC‑taa                                                                                                                                 >  1:466347/1‑148 (MQ=255)
                         tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATcc                                                                                                                              <  1:340077/150‑1 (MQ=255)
                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCg                                                                                                                             <  2:963512/150‑1 (MQ=255)
                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCGCCCCCCCCCCGGGCGGTGCTTTCC‑AACG‑GTC‑CAAACCGc                                                                                                                            >  2:54426/1‑151 (MQ=255)
                            cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGGTGGTTTGC‑AACG‑GTT‑AAAtaccgcc                                                                                                                          >  1:62042/1‑145 (MQ=255)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGCTGCTTCGC‑AACG‑GTC‑AAAtacccccc                                                                                                                         >  1:213208/1‑144 (MQ=255)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGC‑AACC‑TTC‑AAAACCGcccc                                                                                                                         >  2:958671/1‑151 (MQ=21)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGCCGTTGTCTTGC‑AACG‑TTT‑AAAACCCCTc                                                                                                                          >  2:1172145/1‑150 (MQ=21)
                                  tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTGGGCCCGCCCGCCCCCCTCCCGGCCGGTGCCTCGT‑TAAC‑TTC‑AAAttcaccccccgc                                                                                                                     >  2:667396/1‑139 (MQ=255)
                                  tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCcggc                                                                                                                     <  1:182718/150‑1 (MQ=255)
                                  tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGTTTTTC‑CACG‑CTC‑AAAATCGCccccc                                                                                                                       >  2:352529/1‑147 (MQ=21)
                                   aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCcggcgg                                                                                                                   <  2:1125651/151‑1 (MQ=255)
                                    aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCCCCCGGCCCCTGCTTTTC‑CACC‑TTC‑CAATcccccccc                                                                                                                        >  1:761545/1‑145 (MQ=21)
                                     tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCTGCCCCCCCCCCCCCCGGCCGGTTCTTCGC‑AACC‑TTC‑AAAACCGCTCCccgcgg                                                                                                                   >  2:453591/1‑149 (MQ=21)
                                       cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCCGCGGGttt                                                                                                               >  2:135194/1‑151 (MQ=18)
                                         ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGCCCCCCCCCCCCCCGGCCGCTGCTTCGC‑CACG‑TTC‑AAAACCGCTCCCCGCGGGttt                                                                                                               >  1:310514/1‑149 (MQ=18)
                                           aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTcc                                                                                                           <  2:817606/151‑1 (MQ=255)
                                              tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCCCCCCCCCCCCGGCCGGTGCTTTGC‑CACG‑TTC‑AAATCACCCCCccgcggggtttttc                                                                                                           >  1:819321/1‑140 (MQ=255)
                                                aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                       <  1:130494/150‑1 (MQ=255)
                                                 aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCCGTGCTTCCC‑ACCG‑CTC‑TAAACCGCTCCCCGCGGACTTGTCCCACTc                                                                                                      >  2:109509/1‑150 (MQ=17)
                                                  tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGGTTCTTCTC‑CAAC‑TTC‑AAATCCCCCCCCCGCGGGGTTTTTCTACTc                                                                                                      >  1:171777/1‑149 (MQ=255)
                                                     tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGGTGCTTTTCAAACG‑TTA‑AAAACCCCCCCCGCGGGGTTTTTCCtctccgggg                                                                                                  >  1:95226/1‑143 (MQ=12)
                                                      aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGga‑ga                                                                                               <  2:986976/151‑1 (MQ=255)
                                                           gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTTTC‑AACGTTTC‑AAACCCCCCCCCGCCGG‑TTTTTCCTCCTCggggg                                                                                                 >  1:1203511/1‑143 (MQ=17)
                                                                    aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGTGTCCTGCCCGCCCCCCCCCGGGCCGTTGCTTCGC‑CAAG‑CTC‑CAATCCGCCCCCCGCGGGGTTGTTCTACTCCGGG‑GCGGGTTTAcc                                                                                      >  2:884447/1‑146 (MQ=255)
                                                                    aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGacaa                                                                                 <  1:1197515/151‑1 (MQ=255)
                                                                      aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGGCGGTGCTTCGC‑AAAA‑TTT‑AAATACGCCCCCCGCGGGGTTTTTCTACTCCGGA‑GAGGGTTCACCCccaaa                                                                                >  2:723175/1‑150 (MQ=9)
                                                                      aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGacaaac                                                                               <  2:935345/151‑1 (MQ=255)
                                                                           aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAg                                                                          <  2:1006259/151‑1 (MQ=255)
                                                                           aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGGTGCTTCGC‑AACG‑TTC‑AAAAACCCCCCCCGCGGGGTTTTTCTTCTCCGGG‑GAGTGTTTCCCCacaaacaaaa                                                                           >  2:852615/1‑148 (MQ=12)
                                                                             ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGa                                                                         <  2:415672/150‑1 (MQ=255)
                                                                             ccTTTCGGATGGCCTTTCGCTTGATTTGAGGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGa                                                                         <  2:722177/150‑1 (MQ=255)
                                                                              cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTGCTTTGC‑AAAA‑TTC‑AAAACCCCCCCCCGCGGGGTTTTCCTACTCCGGA‑GAGTGTTTACCCCCAAaaaaca                                                                           >  2:51958/1‑147 (MQ=255)
                                                                              cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCCCCACCCCCCCGGCCGTTGCTTCCC‑CACG‑TTT‑TAAAACCCTCCCCGCGGGGTTGTCCCACCCCGGA‑GAGAGTTTACCCacaaacaaaaaaa                                                                        >  2:1100447/1‑145 (MQ=255)
                                                                               tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCTGCCCGCCCCCCCCCCGGGCGGTGCTTTGC‑CACC‑TTT‑AAAACCGCCCCCCGCGGGGTTTTTCTTCTCCGGG‑GGGTGTTTACCGacaaacaaaaaaaa                                                                       >  1:540736/1‑144 (MQ=255)
                                                                                ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATaa                                                                      <  2:342596/150‑1 (MQ=255)
                                                                                 tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATaaa                                                                     <  1:266908/150‑1 (MQ=255)
                                                                                  cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGGTGCTTCGC‑AACG‑TTC‑AAAACCGCCCCCCGCGGGTTTGTCCTTCTCCGGG‑GAGGGGTTACCGCCAaaaaaaaaaaaaaa                                                                    >  1:48551/1‑150 (MQ=255)
                                                                                   ggATGGCCTCCCGCTTCCTTGGATGTCTGGCCGTTTCTGGCGGCCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCACTCCCGGCGGATTTGTCCTACTCGGGA‑GAGAGTTCACCGACAAACAACAGATaaaa                                                                    <  1:31727/149‑1 (MQ=11)
                                                                                     aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACaaa                                                                <  2:792738/151‑1 (MQ=255)
                                                                                     aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCGCCCCCCCCCCGGGCCTTTCCTCTC‑CACC‑TTT‑AAACCCCCCCCCCGCCGGTGTTTTCCCCTCCGGT‑GGGTGTTTACCCCCCAACAACAGAaaaaaaaa                                                                 >  1:402247/1‑150 (MQ=255)
                                                                                       ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGGTGCTTCGC‑CACG‑TTC‑CAATCCGCTCCCCGCGGCGTTGTTCTTCCCGGGG‑GAGGGTTTTCCGCAAAAAAACAaaaaaaaaaaaa                                                               >  1:162300/1‑150 (MQ=255)
                                                                                         ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCCGATGTTCTCGCCACCCCCCCGGCCGTTGCTTCTC‑CATT‑CTCAAAATCCCTTCCCCCCGGATTTTTCCTCCTC‑TGA‑TAGGGTTCTCCGGCAAACCAAAAATAAAAAAAAAggg                                                            >  2:1250587/1‑150 (MQ=255)
                                                                                         ccTTTCGCTTGATTTGATGTCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                             <  1:386194/150‑1 (MQ=255)
                                                                                          cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCGTTGCTTTTC‑AACC‑TTC‑AAAACCCCTCCCCGCGGGATTGTTCTACCCGGGA‑GAGGGTTTACCGACAAAAAAAAGAAAAAACaaaaaggc                                                           >  2:456111/1‑147 (MQ=14)
                                                                                           tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCCGGCCGTTTCTTCGC‑AACG‑TTC‑AAATCCGCCCCCGGCCGGGTTTTCCCTCTCCGGA‑GAGGGTTTACCGCCCAACAACACATAaaaaaaaaa                                                              >  2:732199/1‑146 (MQ=14)
                                                                                               gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                      <  1:528955/151‑1 (MQ=40)
                                                                                               gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                      <  2:659310/151‑1 (MQ=40)
                                                                                                cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCTGCCCCCCCCCCCCCGGGGCGCTGCTTCGC‑AAAG‑CTC‑AAATCCCCCCCCCGCGGGTTTGTCCTACTCCGGA‑GAGGGTTTCCCCCCAAAAAACAGAAAAAAAAAAAAGCCCaagc                                                      >  2:461029/1‑147 (MQ=255)
                                                                                                   gatttgatGTCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                   <  1:993765/150‑1 (MQ=255)
                                                                                                    atttgatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCTCCCGGCCGTTGCTTTGC‑AAAG‑CTC‑AAATCCCCCCCCGGCGGGATTTTCCCCCTCCGGA‑GAGTGTTCCCCGCCACACAACAGATAAAACaaaaaaacccggctt                                                    >  2:1054892/1‑137 (MQ=2)
                                                                                                       tgatGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCGGGCGGTGCTTTGC‑CAAC‑CTC‑AAAAACGCCCCCCGCCGGTTTTTTCTTCTCCGGAGGGGGGTTACCCCCCAAACAAAAACAAAAAAAAAAGCCCCAGCCTTcc                                                  >  1:1184420/1‑148 (MQ=255)
                                                                                                               ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                        <  1:372118/149‑1 (MQ=35)
                                                                                                               ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                        <  1:471890/149‑1 (MQ=35)
                                                                                                                gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                     <  1:212598/151‑1 (MQ=37)
                                                                                                                    tttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                  <  1:299412/150‑1 (MQ=38)
                                                                                                                    tttATGGCAGGCATCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                  <  2:628141/150‑1 (MQ=255)
                                                                                                                       aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                               <  1:53833/150‑1 (MQ=38)
                                                                                                                       aTGGCGGGCGCCCTCCCCGCCGCCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAGCAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                               <  1:940708/150‑1 (MQ=255)
                                                                                                                         ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                             <  2:1231915/150‑1 (MQ=38)
                                                                                                                          gcgggcgTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                           <  2:212297/151‑1 (MQ=38)
                                                                                                                            gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                         <  2:299184/151‑1 (MQ=38)
                                                                                                                            gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                         >  1:1207992/1‑151 (MQ=38)
                                                                                                                             ggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg                        <  2:708966/151‑1 (MQ=38)
                                                                                                                              gcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                       <  2:699456/151‑1 (MQ=39)
                                                                                                                                gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                      <  1:626837/150‑1 (MQ=40)
                                                                                                                                  ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                    <  1:1048614/150‑1 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                  <  2:109520/149‑1 (MQ=40)
                                                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt              >  2:1215447/1‑151 (MQ=40)
                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt              >  2:522685/1‑150 (MQ=40)
                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑GTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc             <  2:711194/151‑1 (MQ=255)
                                                                                                                                         gCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc            <  1:735982/151‑1 (MQ=40)
                                                                                                                                             cccTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa         <  2:540752/150‑1 (MQ=40)
                                                                                                                                                 ccGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct     <  2:116756/150‑1 (MQ=40)
                                                                                                                                                  cGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCg   <  1:1150127/151‑1 (MQ=40)
                                                                                                                                                   gggCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGc  <  1:1087253/151‑1 (MQ=40)
                                                                                                                                                   |                                                                                                                                                          
AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGC  >  NC_000913/3423150‑3423447

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