Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,375,152	0	C	G	41.7%	142.9 / 21.4	72	A95G (GCC→GGC)	dcuD	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (6/36);  new base G (30/0);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 2.37e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAGCGCCTCCGGTCTGGGTATGTTGC  >  NC_000913/3375007‑3375289                                                                                                                                                  |                                                                                                                                           gtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTGCCGCCACATGGAGCGCATAAGGGGCAg                                                                                                                                         >  2:876247/1‑149 (MQ=255)     tCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCCCTACATGGAGCGCATAGGGGGCAGGcg                                                                                                                                      >  1:228580/1‑148 (MQ=255)       gATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGTTATGCCCGCCACATGGAGGGCAAAAGGGGCAggcgggcg                                                                                                                                  >  2:286556/1‑150 (MQ=255)          cccACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCACCTGGAGCGCCTAGGGGGCCGGCGCGCGGTgg                                                                                                                              >  1:533563/1‑151 (MQ=255)          cccACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATg                                                                                                                               <  2:234068/150‑1 (MQ=255)             accaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCCCCACATGGAGGGCATAGGGGGCCGGCGCGCGATgg                                                                                                                              >  1:295822/1‑148 (MQ=255)                 aaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGTTATTCCCGCTACATGGAGGGGCTAGGGGGCAGTCGCGCGGTGGGGAGCt                                                                                                                        >  1:679006/1‑150 (MQ=255)                 aaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTACATGGGGCGCATAGGGGGCAGTCGCGCGATGGGGAGCtt                                                                                                                       >  1:149847/1‑151 (MQ=255)                  aaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGTTATTCCCCCTACATGGGGGGCATAAGGGGCCGGcgcggggagggggggtt                                                                                                                       >  1:58653/1‑137 (MQ=255)                   agcagcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTTCATTGGGCGGATAGGGGGCAATCGCGCGGTGGTGAGCTTGt                                                                                                                     >  1:166679/1‑151 (MQ=255)                    gcagcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGtt                                                                                                                    <  2:920341/151‑1 (MQ=255)                        cGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAgc                                                                                                                <  1:606998/151‑1 (MQ=255)                             aCTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGGGGGGGAGCTTGTTAAACCGccc                                                                                                           >  1:52500/1‑151 (MQ=255)                             aCTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCGCCACATGGGGCGCGTTGGGGGCAGTCGCGCGGGGGGGAGGTTGTTTAAGCGccc                                                                                                           >  2:193397/1‑151 (MQ=255)                             aCTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAgccgc                                                                                                             <  1:170260/149‑1 (MQ=255)                               ttGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCACATGGAGCGCGTTGGGGGCAGTCGCGCGACGGGGGGCTTTTTTAGACGccc                                                                                                           >  2:4654/1‑149 (MQ=255)                                tGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCg                                                                                                          <  1:817515/149‑1 (MQ=255)                                tGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCg                                                                                                          <  1:679747/149‑1 (MQ=255)                                tGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGt                                                                                                         <  1:525454/150‑1 (MQ=255)                                 gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGTCCGCCACATGGGGGGCCTAGGGGGCAGTCGCGCGGTGGGGAGGTTTTTAAGACGccc                                                                                                           >  1:819483/1‑147 (MQ=255)                                 gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGtt                                                                                                        <  2:166844/150‑1 (MQ=255)                                 gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTa                                                                                                       <  1:381077/151‑1 (MQ=255)                                 gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCGCCACATAGAGCGCATAAGGGGCGGGCGCGCGGTGGTGGGCTTGTTTAGCCGccc                                                                                                           >  2:716324/1‑147 (MQ=255)                                   tcgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCCCTACATGGAGCGCCTAGGGGGCCGTCGCGCGGGGGGGAGCTTTTTAAACCGCCCGCTaa                                                                                                      >  1:895249/1‑150 (MQ=255)                                   tcgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTaa                                                                                                      <  2:122316/150‑1 (MQ=255)                                   tcgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCGCGACATGGAGGGGCAAGGGGGCAGGCGCGCGGGGGGGGGGTTGTTAAGACGCCCGtt                                                                                                        >  2:64158/1‑148 (MQ=255)                                    cgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTaaa                                                                                                     <  1:346226/150‑1 (MQ=255)                                      tcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCGCCACATTGGGCGCATAGGGGCCAGTCGGGCGATGGTGAGCGTGGTTAGCCGCCCCTTaaaaa                                                                                                   >  2:655247/1‑149 (MQ=255)                                       cgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTc                                                                                                 <  1:259159/151‑1 (MQ=255)                                        gATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCa                                                                                                <  1:779294/151‑1 (MQ=255)                                         atatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCACCTGGAGCGCCTAGGGGGCAGGCGCGCGGTGGTGAGCTTGTTAAGCCGCCCGTTaaaa                                                                                                    >  1:258255/1‑146 (MQ=255)                                           atTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtc                                                                                             <  1:923519/151‑1 (MQ=255)                                             ttATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcg                                                                                            <  2:525120/150‑1 (MQ=255)                                               ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgct                                                                                          <  2:250593/150‑1 (MQ=255)                                               ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCCACAAGGGGCGCGTAGGGGGCAGTCGCGCGGTGGGGGGGTTGTTTAGACGCCCCGTAAAAATCCttttctc                                                                                         >  2:953239/1‑146 (MQ=255)                                                     aaaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTAAGCCCCCGCCACGGCGAGCCCCAGGGGCACTCGCGCGATGGGGGGGTTGGTTAGACCCCCCTGAAAACACAtttcct                                                                                          >  1:854596/1‑140 (MQ=255)                                                      aaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGtat                                                                                  <  2:923640/151‑1 (MQ=255)                                                         ccTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCCACATGGGGGGCCTAGGGGGCAGTCGCGCGGTGGGGGGGTTTTTTAGGCGCCCGCTTAAAATCATTTCCTCCCCCGAtatt                                                                               >  2:866280/1‑151 (MQ=255)                                                            gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTCCCGCCACATGGGGGGCATAGGGGGCCGGCGCGCGGTGGGGGGCTTGTTAAACCGCCCGGTTAAAATCATTTGCTCGCCGGATattatt                                                                            >  1:647223/1‑151 (MQ=255)                                                            gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTAtattat                                                                             <  1:320463/150‑1 (MQ=255)                                                            gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTAtattat                                                                             <  2:132495/150‑1 (MQ=255)                                                             cgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGGCGCGCGGTGGGGGGCTTGTTAAACCGCCCCTTAAAAATCATTCCCTCCCCGTATattatt                                                                            >  1:527001/1‑150 (MQ=255)                                                              gcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCACATGGAGCGCATAGGGGGCAGTCGCGCGATGGGGAGGGTGTTAAACCGCCCGCTAAAAATCATTCGCCCGCCGGATattatttt                                                                          >  1:787435/1‑149 (MQ=255)                                                               cATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGCTATGTTCGGCACCTTGTGGGGGTAGGGGGCGGTCCCGCGGTGGGGAGCTTGTTAAGCGGCCCGATCATACTTCATCGCCCGCCGtttttt                                                                               >  2:338525/1‑140 (MQ=255)                                                                  gCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCg                                                                      <  2:35100/151‑1 (MQ=255)                                                                  gCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTTTTCCCGCCACATGGGGCGCATAGGGGGCACGGGCGCGGTGGGGGGCTTGTTAAAACGCCCCTTTAAAAACATTTGCCCGCCGCGTattattttctt                                                                       >  2:423733/1‑145 (MQ=255)                                                                     gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCTGTATATTATTCTGTCgg                                                                     <  2:624818/149‑1 (MQ=255)                                                                     gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCa                                                                   <  1:268476/151‑1 (MQ=255)                                                                     gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTGCCGCCACCTGGAGCGCCTAGGGGGCAGGCGCGCGATGGGGAGGTTGTTAAGCCGCCCGGTAAAAATCAATCGCTCGCCGGATATTATTTTGTCgg                                                                     >  1:145208/1‑149 (MQ=255)                                                                        cAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGTTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCACCt                                                                <  2:243455/151‑1 (MQ=255)                                                                         aaCCGCATTGCCGTAGTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACtt                                                               <  1:636359/151‑1 (MQ=255)                                                                                  gCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCACATGGAGCGCATAGGGGGCCGTCGCGCGATGGGGGGCTTGTTAAACCGCCCGTTAAAAATCATTCCCTCGCCGGATATTTTTTTGTCGGCCAccttcctcctc                                                       >  1:346048/1‑140 (MQ=255)                                                                                  gCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATc                                                       <  1:217416/150‑1 (MQ=255)                                                                                    cGGATTGGGGCTGTCGATTATGGCGGGGGGCGGGTTTGCCCGCCACATTGTGGGCGTAGGGGGCAGTCCCGCGAGGGGGGGGGTGTTTAGCCGCCCGCTAAAAATCATTCCCCCCCCGTATATTTTTCTTTCGGCGACTTTCGTCATc                                                       >  2:487823/1‑148 (MQ=255)                                                                                      gATTGGGGCTGTCGATTATGGCGGTTGGGGGGTTTGCCCGCTACATGGAGCGGCCAGGGGGCGGGCGCGCGAGGGGGGGGTTGGTTTGCCGCCCGTTTAAACACATTTCGTCGCCGGTTATTATTCTTTGGGCGACATACGTCATC‑gg                                                    >  2:525892/1‑148 (MQ=255)                                                                                      gATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCa                                                 <  2:678901/151‑1 (MQ=255)                                                                                      gATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCa                                                 <  2:775752/151‑1 (MQ=255)                                                                                      gATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTACATGGAGCGCCTAGGGGGCAGTCGCGCGGTGGGGAGCTTGTTTAGCCGCCCGTTAAAAATCATTCGCTCGCCGTATATTATTTTTTCGGCAAATTACGTCATC‑GGcc                                                  >  1:574172/1‑150 (MQ=255)                                                                                       aTTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCaa                                                <  2:464720/151‑1 (MQ=255)                                                                                       aTTGGGGCTGTCGATTATGGCGGTGGGCGGGTTTTCCCGCCACATGGAGCGCCTAGGGGGCAGTCGCGCGGTGGGGAGGTTGTTAAGCCGCCCGTTTAAAATCATTTGCTCCCCGTATATTTTTTTTTCGGGAACTTTAGTCATC‑CGcc                                                  >  1:34605/1‑149 (MQ=255)                                                                                          ggggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCaaa                                               <  2:419597/149‑1 (MQ=255)                                                                                            ggCTGTCGATTATGGCGGTGGGCGGGTTTTGCCGCCACAAGGGGGGCCTATGGGGCCGGCGCCCGGGGGGGGGGTTGTTAAAGCGCCCGTTAAAAATCCTTCGCTCCCCGGATAATTTTTTTTCGGCCACTTAAGGCATCGGGCCaaaca                                             >  2:784827/1‑148 (MQ=255)                                                                                             gCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATg                                          <  2:92474/151‑1 (MQ=255)                                                                                               tGTCGATTATGGCGGTGGGCGGGTATGGCCGCCCCATGGAGGGCATAGGGGCCAGGCGCGCGGTGGGGGGGTTGTTTAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTTTTTTTTCGGCAACTTACGTCATC‑gggccaa                                               >  2:281447/1‑139 (MQ=255)                                                                                                         ggCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTttat                               <  1:512914/150‑1 (MQ=255)                                                                                                           cGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTAc                            <  1:446440/151‑1 (MQ=255)                                                                                                             gTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCa                          <  1:876371/151‑1 (MQ=255)                                                                                                               gggCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAgc                        <  2:713162/151‑1 (MQ=255)                                                                                                               gggCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAgc                        >  1:563592/1‑151 (MQ=255)                                                                                                                 gCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAgcgc                      <  1:815862/151‑1 (MQ=255)                                                                                                                                    ggAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAGCGCCTCCGGTCTGGGTATGTTg   <  1:286596/151‑1 (MQ=255)                                                                                                                                     gAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAGCGCCTCCGGTCTGGGTATGTTGc  >  2:561191/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                           GTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATC‑GGCCAAATCATGGCGCAGTTTATTACCAGCGCCTCCGGTCTGGGTATGTTGC  >  NC_000913/3375007‑3375289

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