Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	45.3%	94.2 / 10.4	65	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (11/24);  new base G (29/0);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 1.73e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133154‑3133451                                                                                                                                                     |                                                                                                                                                         aTGTATTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtgt                                                                                                                                                        <  2:477689/151‑1 (MQ=255)      gTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGccc                                                                                                                                                   >  1:537075/1‑151 (MQ=255)           gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCa                                                                                                                                              <  2:895485/151‑1 (MQ=255)             tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGt                                                                                                                                            >  2:87696/1‑151 (MQ=255)               cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTcc                                                                                                                                               >  2:484820/1‑146 (MQ=255)                ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCa                                                                                                                                         <  2:1004364/151‑1 (MQ=255)                    cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                     <  1:806968/151‑1 (MQ=255)                    cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCACCCAGTTCAAg                                                                                                                                       >  2:32665/1‑149 (MQ=255)                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑t                                                                                                                                   <  2:665513/151‑1 (MQ=255)                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑t                                                                                                                                   <  1:196655/151‑1 (MQ=255)                      gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGA‑Tc                                                                                                                                  >  2:269091/1‑151 (MQ=255)                      gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAgg                                                                                                                                      >  1:27396/1‑148 (MQ=255)                       ttGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCAAGGA‑TCt                                                                                                                                 >  1:928956/1‑151 (MQ=255)                             tATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGA‑TCttt                                                                                                                               >  1:619954/1‑147 (MQ=255)                               ttAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGGCCCTCCCGGTCAAGGA‑ACTTTTCGcc                                                                                                                          >  2:764613/1‑150 (MQ=255)                                 aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGG‑TCTTTGCGCCAtt                                                                                                                       >  2:218933/1‑151 (MQ=255)                                 aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCCCTTCCAGGG‑TCtttttgccc                                                                                                                         >  2:335316/1‑143 (MQ=255)                                   cAGTTTGAGTTTGTAGTGCATCCATTCGTTACATCCTCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCCAGGA‑TTTTTTCTCCCTTgg                                                                                                                     >  1:912305/1‑150 (MQ=255)                                    agtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGaa                                                                                                                    <  2:1098714/151‑1 (MQ=255)                                     gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAAt                                                                                                                   <  1:567554/151‑1 (MQ=255)                                       ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAGGA‑TCTTTTCGCCCTTTAATGt                                                                                                                 >  1:1079048/1‑151 (MQ=255)                                             ttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGA‑TCTTTTCGCCATTGAATGTCAg                                                                                                              >  1:11977/1‑148 (MQ=255)                                              tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGGTTAAGGA‑TCTTTTCGCCCTTGAAAGTCAGtt                                                                                                            >  1:122761/1‑149 (MQ=255)                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTgcg                                                                                                         <  1:358443/151‑1 (MQ=255)                                                  gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTgcggc                                                                                                       <  1:1088007/150‑1 (MQ=255)                                                   tGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGA‑TCTTTGCGCCCTTGGATGTCAGTTgcggcg                                                                                                      >  2:185742/1‑150 (MQ=255)                                                    gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTa                                                                                                    <  2:938987/151‑1 (MQ=255)                                                     cATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCCTCCCACCAGTTCAAAGA‑TCTTTTCGCCACTTGATGTCGGATgcgg                                                                                                        >  1:283469/1‑146 (MQ=255)                                                        ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                                <  2:527054/151‑1 (MQ=255)                                                        ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                                <  2:2297/151‑1 (MQ=255)                                                            tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGGGTAACGGCg                                                                                              >  1:105699/1‑149 (MQ=255)                                                            tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGACGCCCGCCCCGTTTAAAGG‑TCTTTTCTCCCCTCAATGTCATCTgcggggcgaagcggg                                                                                             >  2:703076/1‑139 (MQ=255)                                                               ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGA‑TCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAACGGCGGtttt                                                                                         >  2:915880/1‑151 (MQ=255)                                                                     catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCCGGTCAAAGG‑TCTTTTCGCCCTTGGATGTCAGTTGCGGGGTAACGGCGAtttttt                                                                                       >  1:717053/1‑147 (MQ=255)                                                                         cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAAc                                                                               <  1:843750/151‑1 (MQ=255)                                                                         cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTGACGGCGATTTTTTTCTGTaa                                                                                >  1:202364/1‑150 (MQ=255)                                                                             aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCt                                                                            <  2:177189/150‑1 (MQ=255)                                                                             aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGA‑TCTTTTCGCCCTTTAATGTCAGTTTCGGCGGAACGGCGGTTTTTTTCTGTAAAGc                                                                             >  1:781093/1‑149 (MQ=255)                                                                               ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGA‑TCTTTGCGCCCTTGAATGTCCGGTGCGGCGGAACGGCGATTTTTTTCTGTTACGCTGcc                                                                         >  2:474641/1‑151 (MQ=255)                                                                                cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGA‑TCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTAACGCTGccc                                                                        >  2:225335/1‑150 (MQ=255)                                                                                    ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAAGG‑TCTTTTCGCCCTTTAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTAAAGCTGcccggt                                                                     >  1:261846/1‑146 (MQ=255)                                                                                        tGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTCAAGGA‑ACTTTTCGCCCTTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAAAGCTGCCCGGtggggg                                                                >  1:436262/1‑151 (MQ=255)                                                                                         gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtggtgg                                                                <  2:318150/150‑1 (MQ=255)                                                                                               gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGA‑TCTTTTCGCCATTGAATGTCAGTTGGGGCGGAAAGGCGGTTTTTTTCTGTTAAGCTTCCCGTTGGGGGTCCgggt                                                         >  2:890566/1‑149 (MQ=255)                                                                                                      cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAGTTCAAGGG‑TCTTTTCGCCACTTTATGTCCGTTTCGGCGGAAAGGCGGTTTTTTTCTTTAAAGCTGCCAGTTGGGGGGTAGGGTT‑TCgcgc                                                 >  1:339192/1‑151 (MQ=255)                                                                                                      cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGGTTAAGGA‑GCTTTTCGCCCTTGTATGGTAGTTGCGGCGTAACGGCGGGTTTTTTCTGTAAAGCCGCCCGCTGGGGGTGAGGGTT‑Tcgcg                                                  >  2:908090/1‑150 (MQ=255)                                                                                                         cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCCTCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTTTTAAGCTGCCAGTTGGTGGGCAGGGTT‑ACGGGCaa                                               >  2:266801/1‑150 (MQ=255)                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCCGTTTAAGGA‑TCTTTTCGCCATTTAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTTTAACGCTGCCCGTTGGGGGGCAGGGTT‑ACGCGCaaa                                              >  1:1023650/1‑148 (MQ=255)                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAAcgc                                           <  2:865567/151‑1 (MQ=255)                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAAcg                                            <  1:474746/150‑1 (MQ=255)                                                                                                                   gcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGGTCAAGGA‑TCTTTTCGCCATTGAATGTCAGTTGCGGCGGAAAGGCGGTTTTTTTCTGTAAAGCTGCCAGTTTGGGGGCAGGGTT‑ACGCGCAAACGCGGGTTGt                                    >  1:449955/1‑151 (MQ=255)                                                                                                                    cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGtggggt                                   <  1:651142/151‑6 (MQ=255)                                                                                                                     gTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTtc                                  <  2:453391/151‑1 (MQ=255)                                                                                                                      tAGTGACGCTTGCAGAATTTCATCGGCGGGGAGGAGCCCATCCAGGTGCAGGA‑GCTTTTCGCCGTTTGATTTGCGGTGTGGGGGAAAGGTGATTTTTTTTTGCAAAGCTTCCTTTTGGAGGTTTGGTTTAACGGCCAAACCGGGTTTGt                                    >  2:982124/1‑149 (MQ=255)                                                                                                                       aGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCACCCAGTTCAAAGG‑TCTTTTCGCCATTGAATGTTAGTTTCGGCGGAACGGCGGTTTTTTTCCGTAACGCTGCCAGGTGGTGGGCAGGGTT‑ACGCGCAAACGCGGgttgtt                                   >  2:950524/1‑148 (MQ=255)                                                                                                                        gTGACGCTTGCAGAATTTCATCGGCGGTGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACTCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTtcat                                <  2:125749/150‑1 (MQ=255)                                                                                                                        gTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCCGATCAAGGA‑TCTTTGCGCCCTTGAATGTCAGTTGCGGCGGAAAGGCGATTTTTTCCTGTAACGCCGTCAGGTGGTGGTTAGGGTT‑TCGCGCCAAAGCGGGTTGTtca                                 >  1:662142/1‑149 (MQ=255)                                                                                                                             gCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGG‑GCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTAACGCTGCCAGGTGGGGGGCCGGGTT‑ACGCGCAAACGCGGGTTGTttat‑cca                            >  2:860871/1‑146 (MQ=255)                                                                                                                                 gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCaa                      <  2:871950/150‑1 (MQ=255)                                                                                                                                   aGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTTAAGGA‑TCTTTTCGCCATTGAATGTCAGTTGCGGCGGAAAGGCGATTTTTTTCTGTAACGCTGCCAGTTGGTGGGCAGGTTT‑TAGCGCCAAAGCGGTTTGTTCAT‑CAATGCCAAGa                    >  2:201109/1‑150 (MQ=255)                                                                                                                                     aaTTTCATCGGCGGGGTGTCGCCCCTCCAGTTTCAGGATTCTTTTC‑CCATTTAATGTCAGGTGGGGCGGAACGGGGACTTTTTTCTTTAACTCTGCCCGCTGGGGGGGAGGTTT‑ACGGGCCAACGCCGGTTGTTTATCCAATGCCa                       >  2:1072767/1‑146 (MQ=255)                                                                                                                                      aTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAGCGCGGTTTGTTCAT‑CAATGCCAAGACGaa                <  1:401032/151‑1 (MQ=255)                                                                                                                                      aTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGaa                <  2:566877/151‑1 (MQ=255)                                                                                                                                          cATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATca             <  2:472471/150‑1 (MQ=255)                                                                                                                                                 ggggggTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATCAACAATACg     >  2:281303/1‑151 (MQ=255)                                                                                                                                                 ggggTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATCAACAATAc      >  1:599813/1‑150 (MQ=255)                                                                                                                                                 ggggTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATCAACAATACg     >  2:844192/1‑151 (MQ=255)                                                                                                                                                     tgtCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATCAACAATACGGGc  >  2:1051315/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                         ATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGA‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTT‑ACGCGCAAACGCGGTTTGTTCAT‑CAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133154‑3133451

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