Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2970GC40.3% 98.2 / 27.8 63intergenic (+103/+126)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (10/27);  new base C (25/0);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.32e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  NC_000913/3423150‑3423443
                                                                                                                                                   |                                                                                                                                                      
aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                                                                                                                                     <  2:1053732/151‑1 (MQ=255)
aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCTTCCCCCCCCCCTCCGGGc                                                                                                                                                     >  1:376398/1‑151 (MQ=255)
 aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                                                                    <  1:247387/151‑1 (MQ=255)
         cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                            <  2:1099589/151‑1 (MQ=255)
            aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGcccgttct                                                                                                                                              >  2:188988/1‑140 (MQ=37)
               aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGCCGGTGCTTCGCaa                                                                                                                                        >  1:1018444/1‑149 (MQ=255)
                    aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑Ca                                                                                                                                 <  2:100990/150‑1 (MQ=255)
                    aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCTCCCGGCCGTTGTTTTGCAACCtt                                                                                                                                    >  1:634263/1‑148 (MQ=255)
                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCGGCCCCGGCTTCTCAACGTT‑CC‑AATACCCTcc                                                                                                                      >  2:1105404/1‑151 (MQ=255)
                              cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCGCCACCCTCCCGGCCGTTGCTTCGCAACGTT‑CA‑AATcccccccc                                                                                                                     >  1:426282/1‑151 (MQ=255)
                               aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGGCCCCCACCCCCCCGGCCGTTGCTTCGCTACGTT‑CT‑AAAACGCcccc                                                                                                                     >  2:315608/1‑150 (MQ=255)
                                ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCCCCCCCCCGGCCGTTGCTTTGCAAACTT‑CA‑AAAcccccccc                                                                                                                     >  1:766052/1‑149 (MQ=255)
                                 gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGTGTCCTTCCCCCCACCCTCCGGGCGGCTGCCTTTCCAGGTT‑CA‑AATCCGCTccc                                                                                                                     >  1:48034/1‑148 (MQ=255)
                                   aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCcggcg                                                                                                                 <  2:266308/150‑1 (MQ=255)
                                   aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCCCCCCCCCGGCCGTTGGTTTGCAACGTT‑CC‑AATCCGCTCCccgcgg                                                                                                                >  1:258591/1‑151 (MQ=255)
                                      ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCCGTGCTTCGCAACGTT‑CA‑AATCCGCTCCccgcgg                                                                                                                >  1:615531/1‑148 (MQ=18)
                                         ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTg                                                                                                          <  2:1072651/150‑1 (MQ=255)
                                         ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTT‑CA‑AATCCCCCCCCCGCGG‑GTTTGt                                                                                                         >  1:208441/1‑151 (MQ=18)
                                         ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTGCCTTGCAACGTT‑CA‑AATCCGCcccc                                                                                                                     >  2:326342/1‑140 (MQ=255)
                                          tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCCCCCTCCCCCCACCCTCCGGCGCGTTCCCTCGCAACCTT‑CA‑AATCCACTCCccgctgtgtttt                                                                                                          >  1:1111890/1‑142 (MQ=255)
                                               gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTAc                                                                                                    <  1:503359/150‑1 (MQ=255)
                                               gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACt                                                                                                   <  1:1109087/151‑1 (MQ=255)
                                                aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCGTTGCTTCGCCACGGT‑CA‑AAACCCCTCCCCGCGG‑GGTTGTTCTtctc                                                                                                  >  1:640442/1‑151 (MQ=18)
                                                     tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGga                                                                                              <  1:599033/150‑1 (MQ=255)
                                                       gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGagag                                                                                           <  2:610728/151‑1 (MQ=255)
                                                        ttCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGCTGCCCGCCCCCCCCCCGGCCCCTGCTTCTCAACTTT‑CA‑AATCCGCCCCCGGCGG‑GGTTGTTCttatccgggg                                                                                             >  2:579318/1‑139 (MQ=255)
                                                          cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAgtg                                                                                         <  1:586417/150‑1 (MQ=255)
                                                            aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCGCAAACTT‑CA‑AAACCCCCCCCCGCGG‑GTTTGTTCTCCTCCGgag                                                                                             >  1:247817/1‑144 (MQ=255)
                                                             aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCCCCCCCCCGGCCGTTGTTTTGCAACGGT‑CA‑AATCCGCCCCCCCCGG‑GGTTGTTCTTCTCCGGagaggg                                                                                         >  1:945153/1‑145 (MQ=255)
                                                                   aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACCTT‑CA‑AATCCCCCCCCCGCGG‑GTTTTTCCCACTCCGGAGAGGGTTTACCCac                                                                                >  1:1039955/1‑150 (MQ=14)
                                                                    aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGGCCTTTCTTTCCAACGTTCCA‑AACCCCCCCCCGCCGG‑GTTT‑TTCTACTCTGGggaggggtgatc                                                                                   >  2:110193/1‑138 (MQ=9)
                                                                     aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGCCGCCCGCCCCCCCCCGGGCCGTTGCTTTGCAAACTT‑CACAATCC‑CCCCCCGCGG‑GCTTTTCCCCCTCCCGCGAGTGTTTACCGa                                                                                 >  2:410860/1‑147 (MQ=9)
                                                                     aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTTCCTTCCAAACTT‑CC‑AATCCCCCCCCGGCGG‑GGTTGTTCTACTCGGGGGGGGGGTTACCGac                                                                                >  2:854936/1‑148 (MQ=12)
                                                                        gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTT‑CA‑AAAACCCCCCCCGCGG‑GTTTTTTCTCCTCCGGGGGGGGTTCAcc                                                                                   >  1:886706/1‑142 (MQ=14)
                                                                          cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                        <  2:169547/151‑1 (MQ=255)
                                                                                 tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                                  <  1:2789/150‑1 (MQ=255)
                                                                                 tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTT‑CA‑AATCCCCCCCCCGCGG‑GTTTTTCCTACTCGGGAGAGTGGTCACCCACAAACAACAGATaaa                                                                  >  1:854364/1‑150 (MQ=14)
                                                                                  cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCGGGCCGGTGCTTCGCAACGTT‑CA‑AATCCGCTCCCCGCGG‑GTTTTTCCTACTCCGGAGAGTGGTTACCGACAAACAACAGATaaaa                                                                 >  1:521110/1‑150 (MQ=14)
                                                                                   ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGGCCCTGCTTCGCAACGTT‑TA‑AATCCGCTCCCCGCGG‑GTTTGTCCCACTCGGGAGAGGGTTTTCCCACAAACCAAAAATaaaa                                                                 >  1:1023985/1‑149 (MQ=14)
                                                                                   ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAAAGTT‑CA‑AATCCCCCCCCCGCGG‑GTTTTTTCTACTCGGGAGAGAGGTCACCCACAAAAAAAAGATaaaa                                                                 >  1:922691/1‑149 (MQ=14)
                                                                                        gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAgg                                                         <  2:1018116/151‑1 (MQ=255)
                                                                                        gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAgg                                                         <  1:621691/151‑1 (MQ=255)
                                                                                              cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAg                                                    <  2:157211/150‑1 (MQ=255)
                                                                                               gCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCCCCTCGCCGGTGCTTTGCAAAGAT‑CA‑AATACCCTCCCCGCGG‑GGTTTTCCTACTTGGGGGAGGGGTTTCCGGCAAAAAAAAGAAAAAAA‑AAAAAGGCCCGTc                                                  >  2:1037756/1‑151 (MQ=255)
                                                                                                   gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTcc                                              <  2:767366/151‑1 (MQ=39)
                                                                                                   gatttgatGTCTGGCAGTTTATGGCGGGCGGGCTGCCCCCCCCCCTCCCGGCCGTTGATTTTCAAAGGT‑TA‑AAAACGCTCCCCGCGG‑GCCTTTCCTCCTCGCGGGAGGGTTCACCGACCAACAAAAGATAaaaa‑aaaaaaaccaacc                                                  >  2:914612/1‑138 (MQ=255)
                                                                                                    atttgatGTCTGGCAGTTTATGGCGGGCGGCCCGGCCCCCCCCCTCCCGGGCGGTGCTTCGCAACGTT‑CA‑AATCCGCTCCCCGCGGCGTTTTTCCT‑CTCCGGAGGGTGGTGAACGACAAAAAACAAAAAAAAC‑AAAAAACCCCGACTTccc                                             >  2:891634/1‑150 (MQ=2)
                                                                                                     tttgatGCCTGGCAGTTTATGGCGGGCGGCCTCCCCGCCACCCTCCCCGCCCGTTCTTTGCCACGTT‑TC‑AATCCGCTCCCCGCGG‑GCTTTTCCTTCCCCAGGGAGGGGTCCCCCCCCAAAAAAAAATAAAACAAAAAGGCCCGGCCttt                                               >  2:603708/1‑148 (MQ=0)
                                                                                                             cTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCCCCGGGCCGTGGCTTTGCAACGTT‑CC‑AATCCGCCCCCGGCGG‑GGTTGTCCTACTCCGCGGAGAGTTCACCCACAAAAAACAGATAAAAC‑AAAAAACCCCGCCTTCCGACTGAGcc                                     >  2:66527/1‑150 (MQ=11)
                                                                                                              tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACGCTCCGGGCCGTTGCTTCGCAACGTT‑GA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGCGAGTTTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCt                                    <  2:629141/150‑1 (MQ=255)
                                                                                                                gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCtttt                                 <  2:748118/151‑1 (MQ=37)
                                                                                                                 cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                <  2:519473/151‑1 (MQ=37)
                                                                                                                           cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                      <  2:857411/151‑1 (MQ=38)
                                                                                                                           cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                      <  1:130324/151‑1 (MQ=38)
                                                                                                                            gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGCTAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                     >  2:275663/1‑151 (MQ=38)
                                                                                                                                 tCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                >  1:18023/1‑151 (MQ=40)
                                                                                                                                  ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg               <  1:559895/151‑1 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc              >  2:542247/1‑149 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa             <  2:770781/150‑1 (MQ=40)
                                                                                                                                          ccACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc       >  1:735742/1‑151 (MQ=40)
                                                                                                                                           cACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt      <  1:1045536/151‑1 (MQ=40)
                                                                                                                                             cccTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa     <  1:295228/150‑1 (MQ=40)
                                                                                                                                                 ccGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc  <  1:579516/149‑1 (MQ=40)
                                                                                                                                                   |                                                                                                                                                      
AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  NC_000913/3423150‑3423443

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