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

AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGC  >  NC_000913/3423150‑3423447
                                                                                                                                                   |                                                                                                                                                          
aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCgg                                                                                                                                                           <  2:1190548/149‑1 (MQ=255)
              gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCCCCCGGCCGGTGCTTCGCaa                                                                                                                                           >  2:336906/1‑150 (MQ=37)
              gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGGTGCTTCGCaa                                                                                                                                           >  2:1149705/1‑150 (MQ=255)
              gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCGGCCGGTGCCTCGCaa                                                                                                                                           >  1:699431/1‑150 (MQ=255)
                      ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCCCAAAGTTCaa                                                                                                                                  >  2:355181/1‑149 (MQ=255)
                      ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                   <  2:803520/148‑1 (MQ=255)
                        ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCCCCCGCCCCCCCCCCGGCCCGTGCGTCCCAACGTTCtaa                                                                                                                                 >  1:466347/1‑148 (MQ=255)
                         tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                              <  1:340077/150‑1 (MQ=255)
                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                             <  2:963512/150‑1 (MQ=255)
                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCGCCCCCCCCCCGGGCGGTGCTTTCCAACGGTCCAAACCGc                                                                                                                            >  2:54426/1‑151 (MQ=255)
                            cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGGTGGTTTGCAACGGTTAAAtaccgcc                                                                                                                          >  1:62042/1‑145 (MQ=255)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGCTGCTTCGCAACGGTCAAAtacccccc                                                                                                                         >  1:213208/1‑144 (MQ=255)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACCTTCAAAACCGcccc                                                                                                                         >  2:958671/1‑151 (MQ=21)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGCCGTTGTCTTGCAACGTTTAAAACCCCTc                                                                                                                          >  2:1172145/1‑150 (MQ=21)
                                  tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTGGGCCCGCCCGCCCCCCTCCCGGCCGGTGCCTCGTTAACTTCAAAttcaccccccgc                                                                                                                     >  2:667396/1‑139 (MQ=255)
                                  tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                     <  1:182718/150‑1 (MQ=255)
                                  tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGTTTTTCCACGCTCAAAATCGCccccc                                                                                                                       >  2:352529/1‑147 (MQ=21)
                                   aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                   <  2:1125651/151‑1 (MQ=255)
                                    aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCCCCCGGCCCCTGCTTTTCCACCTTCCAATcccccccc                                                                                                                        >  1:761545/1‑145 (MQ=21)
                                     tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCTGCCCCCCCCCCCCCCGGCCGGTTCTTCGCAACCTTCAAAACCGCTCCccgcgg                                                                                                                   >  2:453591/1‑149 (MQ=21)
                                       cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGGttt                                                                                                               >  2:135194/1‑151 (MQ=18)
                                         ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGCCCCCCCCCCCCCCGGCCGCTGCTTCGCCACGTTCAAAACCGCTCCCCGCGGGttt                                                                                                               >  1:310514/1‑149 (MQ=18)
                                           aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                           <  2:817606/151‑1 (MQ=255)
                                              tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCCCCCCCCCCCCGGCCGGTGCTTTGCCACGTTCAAATCACCCCCccgcggggtttttc                                                                                                           >  1:819321/1‑140 (MQ=255)
                                                aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                       <  1:130494/150‑1 (MQ=255)
                                                 aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCCGTGCTTCCCACCGCTCTAAACCGCTCCCCGCGGACTTGTCCCACTc                                                                                                      >  2:109509/1‑150 (MQ=17)
                                                  tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGGTTCTTCTCCAACTTCAAATCCCCCCCCCGCGGGGTTTTTCTACTc                                                                                                      >  1:171777/1‑149 (MQ=255)
                                                     tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGGTGCTTTTCAAACGTTAAAAACCCCCCCCGCGGGGTTTTTCCtctccgggg                                                                                                  >  1:95226/1‑143 (MQ=12)
                                                      aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                               <  2:986976/151‑1 (MQ=255)
                                                           gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTTTCAACGTTTCAAACCCCCCCCCGCCGGTTTTTCCTCCTCggggg                                                                                                 >  1:1203511/1‑143 (MQ=17)
                                                                    aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGTGTCCTGCCCGCCCCCCCCCGGGCCGTTGCTTCGCCAAGCTCCAATCCGCCCCCCGCGGGGTTGTTCTACTCCGGGGCGGGTTTAcc                                                                                      >  2:884447/1‑146 (MQ=255)
                                                                    aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                                 <  1:1197515/151‑1 (MQ=255)
                                                                      aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGGCGGTGCTTCGCAAAATTTAAATACGCCCCCCGCGGGGTTTTTCTACTCCGGAGAGGGTTCACCCccaaa                                                                                >  2:723175/1‑150 (MQ=9)
                                                                      aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                               <  2:935345/151‑1 (MQ=255)
                                                                           aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAg                                                                          <  2:1006259/151‑1 (MQ=255)
                                                                           aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGGTGCTTCGCAACGTTCAAAAACCCCCCCCGCGGGGTTTTTCTTCTCCGGGGAGTGTTTCCCCacaaacaaaa                                                                           >  2:852615/1‑148 (MQ=12)
                                                                             ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa                                                                         <  2:415672/150‑1 (MQ=255)
                                                                             ccTTTCGGATGGCCTTTCGCTTGATTTGAGGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa                                                                         <  2:722177/150‑1 (MQ=255)
                                                                              cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTGCTTTGCAAAATTCAAAACCCCCCCCCGCGGGGTTTTCCTACTCCGGAGAGTGTTTACCCCCAAaaaaca                                                                           >  2:51958/1‑147 (MQ=255)
                                                                              cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCCCCACCCCCCCGGCCGTTGCTTCCCCACGTTTTAAAACCCTCCCCGCGGGGTTGTCCCACCCCGGAGAGAGTTTACCCacaaacaaaaaaa                                                                        >  2:1100447/1‑145 (MQ=255)
                                                                               tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCTGCCCGCCCCCCCCCCGGGCGGTGCTTTGCCACCTTTAAAACCGCCCCCCGCGGGGTTTTTCTTCTCCGGGGGGTGTTTACCGacaaacaaaaaaaa                                                                       >  1:540736/1‑144 (MQ=255)
                                                                                ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa                                                                      <  2:342596/150‑1 (MQ=255)
                                                                                 tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                                     <  1:266908/150‑1 (MQ=255)
                                                                                  cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGGTGCTTCGCAACGTTCAAAACCGCCCCCCGCGGGTTTGTCCTTCTCCGGGGAGGGGTTACCGCCAaaaaaaaaaaaaaa                                                                    >  1:48551/1‑150 (MQ=255)
                                                                                   ggATGGCCTCCCGCTTCCTTGGATGTCTGGCCGTTTCTGGCGGCCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCACTCCCGGCGGATTTGTCCTACTCGGGAGAGAGTTCACCGACAAACAACAGATaaaa                                                                    <  1:31727/149‑1 (MQ=11)
                                                                                     aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                                <  2:792738/151‑1 (MQ=255)
                                                                                     aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCGCCCCCCCCCCGGGCCTTTCCTCTCCACCTTTAAACCCCCCCCCCGCCGGTGTTTTCCCCTCCGGTGGGTGTTTACCCCCCAACAACAGAaaaaaaaa                                                                 >  1:402247/1‑150 (MQ=255)
                                                                                       ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGGTGCTTCGCCACGTTCCAATCCGCTCCCCGCGGCGTTGTTCTTCCCGGGGGAGGGTTTTCCGCAAAAAAACAaaaaaaaaaaaa                                                               >  1:162300/1‑150 (MQ=255)
                                                                                         ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCCGATGTTCTCGCCACCCCCCCGGCCGTTGCTTCTCCATTCTCAAAATCCCTTCCCCCCGGATTTTTCCTCCTCTGATAGGGTTCTCCGGCAAACCAAAAATAAAAAAAAAggg                                                            >  2:1250587/1‑150 (MQ=255)
                                                                                         ccTTTCGCTTGATTTGATGTCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                             <  1:386194/150‑1 (MQ=255)
                                                                                          cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCGTTGCTTTTCAACCTTCAAAACCCCTCCCCGCGGGATTGTTCTACCCGGGAGAGGGTTTACCGACAAAAAAAAGAAAAAACaaaaaggc                                                           >  2:456111/1‑147 (MQ=14)
                                                                                           tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCCGGCCGTTTCTTCGCAACGTTCAAATCCGCCCCCGGCCGGGTTTTCCCTCTCCGGAGAGGGTTTACCGCCCAACAACACATAaaaaaaaaa                                                              >  2:732199/1‑146 (MQ=14)
                                                                                               gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                      <  1:528955/151‑1 (MQ=40)
                                                                                               gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                      <  2:659310/151‑1 (MQ=40)
                                                                                                cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCTGCCCCCCCCCCCCCGGGGCGCTGCTTCGCAAAGCTCAAATCCCCCCCCCGCGGGTTTGTCCTACTCCGGAGAGGGTTTCCCCCCAAAAAACAGAAAAAAAAAAAAGCCCaagc                                                      >  2:461029/1‑147 (MQ=255)
                                                                                                   gatttgatGTCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                   <  1:993765/150‑1 (MQ=255)
                                                                                                    atttgatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCTCCCGGCCGTTGCTTTGCAAAGCTCAAATCCCCCCCCGGCGGGATTTTCCCCCTCCGGAGAGTGTTCCCCGCCACACAACAGATAAAACaaaaaaacccggctt                                                    >  2:1054892/1‑137 (MQ=2)
                                                                                                       tgatGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCGGGCGGTGCTTTGCCAACCTCAAAAACGCCCCCCGCCGGTTTTTTCTTCTCCGGAGGGGGGTTACCCCCCAAACAAAAACAAAAAAAAAAGCCCCAGCCTTcc                                                  >  1:1184420/1‑148 (MQ=255)
                                                                                                               ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                        <  1:372118/149‑1 (MQ=35)
                                                                                                               ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                        <  1:471890/149‑1 (MQ=35)
                                                                                                                gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                     <  1:212598/151‑1 (MQ=37)
                                                                                                                    tttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                  <  1:299412/150‑1 (MQ=38)
                                                                                                                    tttATGGCAGGCATCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                  <  2:628141/150‑1 (MQ=255)
                                                                                                                       aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                               <  1:53833/150‑1 (MQ=38)
                                                                                                                       aTGGCGGGCGCCCTCCCCGCCGCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAGCAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                               <  1:940708/150‑1 (MQ=255)
                                                                                                                         ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                             <  2:1231915/150‑1 (MQ=38)
                                                                                                                          gcgggcgTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                           <  2:212297/151‑1 (MQ=38)
                                                                                                                            gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                         <  2:299184/151‑1 (MQ=38)
                                                                                                                            gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                         >  1:1207992/1‑151 (MQ=38)
                                                                                                                             ggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg                        <  2:708966/151‑1 (MQ=38)
                                                                                                                              gcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                       <  2:699456/151‑1 (MQ=39)
                                                                                                                                gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                      <  1:626837/150‑1 (MQ=40)
                                                                                                                                  ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                    <  1:1048614/150‑1 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                  <  2:109520/149‑1 (MQ=40)
                                                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt              >  2:1215447/1‑151 (MQ=40)
                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt              >  2:522685/1‑150 (MQ=40)
                                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc             <  2:711194/151‑1 (MQ=255)
                                                                                                                                         gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc            <  1:735982/151‑1 (MQ=40)
                                                                                                                                             cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa         <  2:540752/150‑1 (MQ=40)
                                                                                                                                                 ccGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct     <  2:116756/150‑1 (MQ=40)
                                                                                                                                                  cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCg   <  1:1150127/151‑1 (MQ=40)
                                                                                                                                                   gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGc  <  1:1087253/151‑1 (MQ=40)
                                                                                                                                                   |                                                                                                                                                          
AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGC  >  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: