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

GCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  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)
                           aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                             <  2:100990/150‑1 (MQ=255)
                           aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCTCCCGGCCGTTGTTTTGCAACCtt                                                                                                                                >  1:634263/1‑148 (MQ=255)
                                    gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCGGCCCCGGCTTCTCAACGTTCCAATACCCTcc                                                                                                                  >  2:1105404/1‑151 (MQ=255)
                                     cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCGCCACCCTCCCGGCCGTTGCTTCGCAACGTTCAAATcccccccc                                                                                                                 >  1:426282/1‑151 (MQ=255)
                                      aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGGCCCCCACCCCCCCGGCCGTTGCTTCGCTACGTTCTAAAACGCcccc                                                                                                                 >  2:315608/1‑150 (MQ=255)
                                       ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCCCCCCCCCGGCCGTTGCTTTGCAAACTTCAAAAcccccccc                                                                                                                 >  1:766052/1‑149 (MQ=255)
                                        gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGTGTCCTTCCCCCCACCCTCCGGGCGGCTGCCTTTCCAGGTTCAAATCCGCTccc                                                                                                                 >  1:48034/1‑148 (MQ=255)
                                          aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                             <  2:266308/150‑1 (MQ=255)
                                          aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCCCCCCCCCGGCCGTTGGTTTGCAACGTTCCAATCCGCTCCccgcgg                                                                                                            >  1:258591/1‑151 (MQ=255)
                                             ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCCGTGCTTCGCAACGTTCAAATCCGCTCCccgcgg                                                                                                            >  1:615531/1‑148 (MQ=18)
                                                ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                      <  2:1072651/150‑1 (MQ=255)
                                                ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCCCCCCCCGCGGGTTTGt                                                                                                     >  1:208441/1‑151 (MQ=18)
                                                ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTGCCTTGCAACGTTCAAATCCGCcccc                                                                                                                 >  2:326342/1‑140 (MQ=255)
                                                 tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCCCCCTCCCCCCACCCTCCGGCGCGTTCCCTCGCAACCTTCAAATCCACTCCccgctgtgtttt                                                                                                      >  1:1111890/1‑142 (MQ=255)
                                                      gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                <  1:503359/150‑1 (MQ=255)
                                                      gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                               <  1:1109087/151‑1 (MQ=255)
                                                       aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCGTTGCTTCGCCACGGTCAAAACCCCTCCCCGCGGGGTTGTTCTtctc                                                                                              >  1:640442/1‑151 (MQ=18)
                                                            tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                          <  1:599033/150‑1 (MQ=255)
                                                              gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGagag                                                                                       <  2:610728/151‑1 (MQ=255)
                                                               ttCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGCTGCCCGCCCCCCCCCCGGCCCCTGCTTCTCAACTTTCAAATCCGCCCCCGGCGGGGTTGTTCttatccgggg                                                                                         >  2:579318/1‑139 (MQ=255)
                                                                 cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                     <  1:586417/150‑1 (MQ=255)
                                                                   aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCGCAAACTTCAAAACCCCCCCCCGCGGGTTTGTTCTCCTCCGgag                                                                                         >  1:247817/1‑144 (MQ=255)
                                                                    aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCCCCCCCCCGGCCGTTGTTTTGCAACGGTCAAATCCGCCCCCCCCGGGGTTGTTCTTCTCCGGagaggg                                                                                     >  1:945153/1‑145 (MQ=255)
                                                                          aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACCTTCAAATCCCCCCCCCGCGGGTTTTTCCCACTCCGGAGAGGGTTTACCCac                                                                            >  1:1039955/1‑150 (MQ=14)
                                                                           aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGGCCTTTCTTTCCAACGTTCCAAACCCCCCCCCGCCGGGTTTTTCTACTCTGGggaggggtgatc                                                                               >  2:110193/1‑138 (MQ=9)
                                                                            aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGCCGCCCGCCCCCCCCCGGGCCGTTGCTTTGCAAACTTCACAATCCCCCCCCGCGGGCTTTTCCCCCTCCCGCGAGTGTTTACCGa                                                                             >  2:410860/1‑147 (MQ=9)
                                                                            aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTTCCTTCCAAACTTCCAATCCCCCCCCGGCGGGGTTGTTCTACTCGGGGGGGGGGTTACCGac                                                                            >  2:854936/1‑148 (MQ=12)
                                                                               gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAAAACCCCCCCCGCGGGTTTTTTCTCCTCCGGGGGGGGTTCAcc                                                                               >  1:886706/1‑142 (MQ=14)
                                                                                 cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                    <  2:169547/151‑1 (MQ=255)
                                                                                        tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                              <  1:2789/150‑1 (MQ=255)
                                                                                        tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCCCCCCCCGCGGGTTTTTCCTACTCGGGAGAGTGGTCACCCACAAACAACAGATaaa                                                              >  1:854364/1‑150 (MQ=14)
                                                                                         cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCCGGAGAGTGGTTACCGACAAACAACAGATaaaa                                                             >  1:521110/1‑150 (MQ=14)
                                                                                          ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGGCCCTGCTTCGCAACGTTTAAATCCGCTCCCCGCGGGTTTGTCCCACTCGGGAGAGGGTTTTCCCACAAACCAAAAATaaaa                                                             >  1:1023985/1‑149 (MQ=14)
                                                                                          ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAAAGTTCAAATCCCCCCCCCGCGGGTTTTTTCTACTCGGGAGAGAGGTCACCCACAAAAAAAAGATaaaa                                                             >  1:922691/1‑149 (MQ=14)
                                                                                               gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                     <  2:1018116/151‑1 (MQ=255)
                                                                                               gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                     <  1:621691/151‑1 (MQ=255)
                                                                                                     cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                <  2:157211/150‑1 (MQ=255)
                                                                                                      gCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCCCCTCGCCGGTGCTTTGCAAAGATCAAATACCCTCCCCGCGGGGTTTTCCTACTTGGGGGAGGGGTTTCCGGCAAAAAAAAGAAAAAAAAAAAAGGCCCGTc                                              >  2:1037756/1‑151 (MQ=255)
                                                                                                          gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                          <  2:767366/151‑1 (MQ=39)
                                                                                                          gatttgatGTCTGGCAGTTTATGGCGGGCGGGCTGCCCCCCCCCCTCCCGGCCGTTGATTTTCAAAGGTTAAAAACGCTCCCCGCGGGCCTTTCCTCCTCGCGGGAGGGTTCACCGACCAACAAAAGATAaaaaaaaaaaaccaacc                                              >  2:914612/1‑138 (MQ=255)
                                                                                                           atttgatGTCTGGCAGTTTATGGCGGGCGGCCCGGCCCCCCCCCTCCCGGGCGGTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGCGTTTTTCCTCTCCGGAGGGTGGTGAACGACAAAAAACAAAAAAAACAAAAAACCCCGACTTccc                                         >  2:891634/1‑150 (MQ=2)
                                                                                                            tttgatGCCTGGCAGTTTATGGCGGGCGGCCTCCCCGCCACCCTCCCCGCCCGTTCTTTGCCACGTTTCAATCCGCTCCCCGCGGGCTTTTCCTTCCCCAGGGAGGGGTCCCCCCCCAAAAAAAAATAAAACAAAAAGGCCCGGCCttt                                           >  2:603708/1‑148 (MQ=0)
                                                                                                                    cTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCCCCGGGCCGTGGCTTTGCAACGTTCCAATCCGCCCCCGGCGGGGTTGTCCTACTCCGCGGAGAGTTCACCCACAAAAAACAGATAAAACAAAAAACCCCGCCTTCCGACTGAGcc                                 >  2:66527/1‑150 (MQ=11)
                                                                                                                     tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACGCTCCGGGCCGTTGCTTCGCAACGTTGAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGCGAGTTTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                <  2:629141/150‑1 (MQ=255)
                                                                                                                       gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                             <  2:748118/151‑1 (MQ=37)
                                                                                                                        cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                            <  2:519473/151‑1 (MQ=37)
                                                                                                                                  cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                  <  1:130324/151‑1 (MQ=38)
                                                                                                                                  cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                  <  2:857411/151‑1 (MQ=38)
                                                                                                                                   gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGCTAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                 >  2:275663/1‑151 (MQ=38)
                                                                                                                                        tCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg            >  1:18023/1‑151 (MQ=40)
                                                                                                                                         ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg           <  1:559895/151‑1 (MQ=40)
                                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          >  2:542247/1‑149 (MQ=40)
                                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa         <  2:770781/150‑1 (MQ=40)
                                                                                                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc   >  1:735742/1‑151 (MQ=40)
                                                                                                                                                  cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt  <  1:1045536/151‑1 (MQ=40)
                                                                                                                                                   |                                                                                                                                                         
GCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  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: