Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2900AC38.5% 70.7 / 18.1 65intergenic (+96/+133)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (14/26);  new base C (25/0);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.17e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  NC_000913/3423143‑3423439
                                                                                                                                                   |                                                                                                                                                         
gcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCGCCAc                                                                                                                                                          >  2:229934/1‑149 (MQ=255)
   tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGGCCGCCCCCCTcc                                                                                                                                                     >  2:140965/1‑151 (MQ=255)
    ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                                                                                     <  1:259638/150‑1 (MQ=255)
     cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCccccccccccccc                                                                                                                                                    >  1:748353/1‑149 (MQ=255)
       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)
                   aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGTCCCCCCCCCCAccccccccggc                                                                                                                                                >  2:799295/1‑135 (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                  <  1:130324/151‑1 (MQ=38)
                                                                                                                                  cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                  <  2:857411/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)
                                                                                                                                                   |                                                                                                                                                         
GCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTT‑CA‑AATCCGCTCCCGGCGG‑ATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAC‑AAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  NC_000913/3423143‑3423439

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: ‑