Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,375,123	0	T	G	43.9%	170.5 / 20.3	98	G85G (GGT→GGG)	dcuD	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (4/51);  new base G (43/0);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.45e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCA  >  NC_000913/3374976‑3375265                                                                                                                                                    |                                                                                                                                               gatgTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTAt                                                                                                                                             <  1:499445/151‑1 (MQ=255)     tgCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGccc                                                                                                                                         >  2:282008/1‑151 (MQ=255)       tggTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTGTGccccc                                                                                                                                       >  1:261158/2‑149 (MQ=255)       cGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGc                                                                                                                                       <  2:806738/151‑1 (MQ=255)        ggTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCt                                                                                                                                      <  1:846978/151‑1 (MQ=255)               ggTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACAt                                                                                                                                  <  2:299823/148‑1 (MQ=255)               ggTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATgg                                                                                                                                <  1:678408/150‑1 (MQ=255)                gTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAg                                                                                                                              <  1:1401462/151‑1 (MQ=255)                   aGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGGTTGCCCGCCACAttg                                                                                                                                >  1:1383406/1‑144 (MQ=255)                      gTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATa                                                                                                                        <  1:1210920/151‑1 (MQ=255)                         cGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAggg                                                                                                                     <  2:717846/151‑1 (MQ=255)                         cGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAggg                                                                                                                     <  2:545893/151‑1 (MQ=255)                             ggtgTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCa                                                                                                                 <  2:875656/151‑1 (MQ=255)                              gtgtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTTCCCCCTACAAggggggccgaggg                                                                                                                     >  1:978087/1‑135 (MQ=255)                               tgtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTTCCGCCACATGGAGCGCACAgggg                                                                                                                    >  1:731626/1‑146 (MQ=255)                                gtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGACCCCCACATGGAGCGCCTGgggggccggc                                                                                                              >  1:239300/1‑145 (MQ=255)                                  aCTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTGTTCCCCCCCCCTGGAGCGGATAgggggccggc                                                                                                              >  1:772061/1‑143 (MQ=255)                                   cTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCCCCTGGAGGGCGTAAGGGGGGGTcgcg                                                                                                           >  1:849351/1‑151 (MQ=255)                                   cTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCCCCCCATGGAGGGCATa                                                                                                                        >  1:1372926/1‑138 (MQ=255)                                     cGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTcgcgcg                                                                                                         <  2:499263/151‑1 (MQ=255)                                     cGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTcgcgc                                                                                                          <  1:1294528/150‑1 (MQ=255)                                          ccaccaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGt                                                                                                    <  2:885712/151‑1 (MQ=255)                                            accaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGa                                                                                                  <  2:1386207/151‑1 (MQ=255)                                             ccaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGGCCGCTACCTGGGGCGCCTAGGGGGCAGGcgcgcggcgggggg                                                                                                 >  1:409765/1‑143 (MQ=255)                                              caAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCTACATGGAGGGGGTAGGGGGCAGGCGCGCGCTGGGGAg                                                                                                 >  2:843362/1‑150 (MQ=255)                                               aaaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCt                                                                                               <  1:321684/151‑1 (MQ=255)                                               aaaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGTCCGCTACATGGAGCGCATAAGGGCCCGCCCCGCGCTGGTGAg                                                                                                 >  2:1247523/1‑149 (MQ=255)                                                aaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCtt                                                                                              <  1:1262220/151‑1 (MQ=255)                                                  agcagcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTTTGCCCCCCACAATGAGCGCATAGGGGCCCGTCGCGCGATGGGGAGGtt                                                                                              >  2:1481758/1‑149 (MQ=255)                                                    cagcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTGCCGCTACCTGGAGGGCATAGGGGGGAGGcgcgcgggggggggctttttt                                                                                          >  1:605987/1‑136 (MQ=255)                                                    cagcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTa                                                                                          <  1:323823/151‑1 (MQ=255)                                                     agcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGGCCGCTCCATGGGGCGCATAGGGGCCCGGCGCGCGAAgggggggtttgt                                                                                           >  2:16462/1‑139 (MQ=255)                                                     agcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGTCCGCCACATGGAGCGCATAGGGGGCCGTCGCGCGGTgggggggttttt                                                                                           >  2:40880/1‑139 (MQ=255)                                                       cGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTAATCCCCCCACATGGGGGGCACAGGGGGCAGTCGGGCGGTGGGGGGGTTGt                                                                                            >  1:1475734/1‑146 (MQ=255)                                                        ggCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAgcc                                                                                      <  2:879429/151‑1 (MQ=255)                                                                gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGtt                                                                               <  2:1114921/150‑1 (MQ=255)                                                                 aTCGTCGATATTTATAATGAAATCCTGCGCCTGCTGTCCATCCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGtt                                                                               <  2:558950/149‑1 (MQ=255)                                                                 aTCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTTTGGCCGCCTCATGGAGGGCATAGGGGGCAGCAGCGCGGTGGGGAGCTTTTTTAACCGccccct                                                                               >  2:330389/1‑146 (MQ=255)                                                                  tcgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTTCCGCCCCCTGGAGGGCATAGGGGGCAGTCGCGGGGGGGGGAGCGTGTTTAGCCGccc                                                                                  >  1:1267953/1‑145 (MQ=255)                                                                   cgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCCCATGGGGCGCCTAGGGGGCAGGCGCGCGGTGGGGAGGTTGTTTAGCCGCCCgcttaa                                                                            >  2:1456893/1‑145 (MQ=255)                                                                   cgtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTaaaa                                                                           <  2:936532/151‑1 (MQ=255)                                                                       gATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGCGGGTTTGCCCCCTCCACTGAGGGCATCAGGGGTAGGCGCGCGAGGGGGGGGCTGTTAAGACCCCCGCTaaaa                                                                           >  1:934114/1‑147 (MQ=255)                                                                        atatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCCCCCCATGGGGCGCCTAGGGGGCCGGCGCGCGATGGGGAGCTTGTTAAACCCCCCGTTtaaa                                                                           >  2:1183955/1‑146 (MQ=255)                                                                        atatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCa                                                                       <  1:37647/150‑1 (MQ=255)                                                                        atatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCAt                                                                      <  1:1198611/151‑1 (MQ=255)                                                                         tatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCCACATGGGGCGCCTAAGGGGCCGTCCAGCGGCGGGGGGGTTGTTTATCCCCCCCCTAAAACACAtt                                                                     >  1:1188540/1‑151 (MQ=255)                                                                         tatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCCCCCCCTGGGGGGCACAAGGGGCCGGCGCGGGCTGGTGGGGTTGTTTAACCCCCCGCTAAAAATCa                                                                       >  1:709491/1‑149 (MQ=255)                                                                          atTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCGCCACCTGGGGCGCATAGGGGGCCGTCGCGCGGTGGGGGGCTTTTTTAACCGcccccttaaaatcctt                                                                     >  2:1171719/1‑137 (MQ=255)                                                                            ttATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcg                                                                   <  1:1506986/150‑1 (MQ=255)                                                                             tataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgc                                                                  <  1:1010034/150‑1 (MQ=255)                                                                             tataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgc                                                                  <  2:1000053/150‑1 (MQ=255)                                                                              ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgctc                                                                <  1:211355/151‑1 (MQ=255)                                                                               taATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTATTCCCCCCCCATGGGGGGCATAGGGGCCAGGCGCGCGGTGGGGGGGTTGTTTAGCCGCCCCCTTAAAATCAtt                                                                     >  1:1380168/1‑145 (MQ=255)                                                                                  tGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCCCCACATGGAGGGGCTAAGGGGCAGGCGCGCGGGGGGGGGGTTGTTTAGCCGCCCCGTAAAAATCATtctcccccc                                                             >  1:675849/1‑143 (MQ=255)                                                                                   gAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGTTATTGCCGCCACATGGGGGGCATAGGGGGCAGCCGCGCGGTGGGGGGCTTGTTTAGACGCCCGTTTAAAATCATTTGCTcccc                                                             >  2:1305415/1‑149 (MQ=255)                                                                                   gAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACCAGGTGGGCAAAAGGGGCCCACGCGCCAAGGTGGGCTTTTTTAATCGCCCCTTAAAAATCATttcctcgtc                                                             >  2:512192/1‑147 (MQ=255)                                                                                   gAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCg                                                            <  1:375669/150‑1 (MQ=255)                                                                                    aaaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTGTGGCCGCCACATGGGGGGCATAGGGGCCAGTCGCGCGGCGGGGGGCTTTTTTAACCGCCCCCTTAAAAACATTTGCTTGcc                                                             >  2:529758/1‑148 (MQ=255)                                                                                     taTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGGCCGCGACATGGAGGGCATAGGGGGCCGTCGCGCGGTGGGGGGGTTGGTATACCGCCCGGTTAAAATCATTCGCTCGCCGtat                                                         >  1:789045/2‑151 (MQ=255)                                                                                     aaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGtat                                                         <  1:1089479/151‑1 (MQ=255)                                                                                      aTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCCCTACATGGAGCGCCTAGGGGGCCGTCGCGCGGTGGTGAGCTTTTTTAACCGCCCCCTTAAACTCATTCCCTccccctct                                                         >  2:618050/1‑146 (MQ=255)                                                                                      aTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCCCTACATGGAGCGCATAGGGGGCCGGCGCGCGATGGGGGGCTTGTTTAGCCCCCCGTTAAAACTCATTCCCTCCCCCtat                                                         >  2:322261/1‑150 (MQ=255)                                                                                       tCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGACCCCCACATGGGGGGCATAGGGGGCAGTCGCGCGATGGGGGGGTTTTTTAACCGCCCGTTAAAAATCATTcgcccgcc                                                             >  2:332795/1‑145 (MQ=255)                                                                                        ccTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGGCCCCCACATTGAGGGCCTAGGGGGCAGTCGCGCGGTGGTGGGGTTGTTTAGCCGCCCCTTTAAAATCATTTGCTCGCCGGAtatt                                                      >  2:1217588/1‑151 (MQ=255)                                                                                        ccTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCTGTatat                                                       <  2:1432596/150‑1 (MQ=255)                                                                                          tGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGGATGCCCGCCACAATGAGGGCCGAGGGGGCAGGCGCGCGGTGGGGGGGTTGTTTAACCGCCCCTTTAAAATCCTTTGCTCGCCCTAtatttt                                                    >  1:447438/1‑149 (MQ=255)                                                                                          tGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGGCCGCCACAAGGGGCGCCTAGGGGGCGGGCGCGCGGGGGTGAGCGTGTTTAACCGCCCGCTTAAACTCATTCGCTCGCCgcgtcttt                                                     >  1:1364554/1‑143 (MQ=255)                                                                                            cgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCCCTCCATTGAGGGGATAAGGGGCAGGCGCGCGAAGGGGAGGGTTTTTAGACGCCCCGTAAAAAACATTCCCTcccccggtattttt                                                   >  1:680702/1‑140 (MQ=255)                                                                                            cgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTCCCCCCACAAGGAGCGCCTAAGGGGCAGGCGCGCGGTGGGGAGGTTGTTTAGACCCCCCGTAAAAAACATTCGCCCCCCGCGTattattt                                                  >  1:1457736/1‑150 (MQ=255)                                                                                             gcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACACATTCGCTCGCCGTATATTATTCt                                                 <  2:294119/151‑1 (MQ=255)                                                                                                tGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGGCCCCTACATGGAGGGCCTAGGGGGCAGGCGCGCGGTGGGGAGGTTGTTTAGCCGCCCGCTTAAAATCATTTGCTCGCCGGATATTATTTTGt                                               >  2:322420/1‑150 (MQ=255)                                                                                                tGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGt                                               <  1:1468116/150‑1 (MQ=255)                                                                                                    gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCGCTACAAGGGGCGCATAGGGGGCGGGCGCGGGATGGGGGGGGTTTTAAACCCCCCCGTAAAAATCATTTTCTCGCCGCATattatttt                                                 >  2:840500/1‑142 (MQ=255)                                                                                                    gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGc                                           <  2:707616/150‑1 (MQ=255)                                                                                                      ccAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCaa                                         <  2:1381919/150‑1 (MQ=255)                                                                                                       cAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATTACCCCCACCTGGAGGGCATACGGGCCAGTCGCGCGATGGTGAGACTGTTAAACCGCCCCCTTAAAATCATTCTCCCACCGGATATTTTTTTGTTGGCAAc                                        >  2:659952/1‑150 (MQ=255)                                                                                                              aTTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTc                                 <  1:533268/150‑1 (MQ=255)                                                                                                                tGCCGGATTGGGGCTGTCGATTATGGCGGTTGGGGGTTATGTCCCCTCCAGGGAGGGCATAGGGGCCCGGCCCCCGGTGGGGAGGTTGTTTAGCCGCCCCTCTAAACACAATCGCTCGCCGGTTATTATTTTTTCgggcacctacctc                                 >  1:950240/1‑137 (MQ=255)                                                                                                                 gCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCCCCACATAGGGCGCCCAGGGGGGAGGGGCGCGGTGGGGGGGTTGTTAAACCGGCCGGTAAAAAATAATCGCCCGCCGGATATTATTTTTTTCGGAACTTAca                                   >  1:691754/1‑144 (MQ=255)                                                                                                                   cGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGc                           <  1:803698/151‑1 (MQ=255)                                                                                                                   cGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGc                           <  1:1456086/151‑1 (MQ=255)                                                                                                                    ggATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGcc                          <  2:1520038/151‑1 (MQ=255)                                                                                                                    ggATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGcc                          <  1:341567/151‑1 (MQ=255)                                                                                                                    ggATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGc                           <  2:1016707/150‑1 (MQ=255)                                                                                                                     gATTGGGGCTGTCGATTATGGCGGTTGGGGGGTTTGTCCGCTACATGGAGCGCCTAGGGGGCCGTCGCGCGAGGGTGGGGTTGTTTAGCCGCCCCTTTAAACTCATTCCCTCGCCGGATATTTTTTTTTCGGGCACTTACGTCATCGGcc                          >  2:534353/1‑150 (MQ=255)                                                                                                                      aTTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCaa                        <  1:979776/151‑1 (MQ=255)                                                                                                                      aTTGGGGCTGTCGATTATGGCGGCGGGGGGGTATTCCCGCCTCATGGGGCGCATAGGGGGCAGGCGGGCGGTGGGGAGGTTGTTAAGCCGCCCGCTTAAACTCATTCGCTCGCCGCATATTTTTCTTTCGGCAACTTTCGTCATCCGcc                          >  1:1193930/1‑149 (MQ=255)                                                                                                                          gggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCa                    <  2:1153584/151‑1 (MQ=255)                                                                                                                          gggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCa                    <  2:586483/151‑1 (MQ=255)                                                                                                                          gggCTGTCGATTATGGCGGTGGGCGGGTTTGCCCGCCTAATTGAGCGGATAGGGGGCACGGGCGCGGGGGGGAGCCTGTTTAACCGCACCTTAAAACTCATTTGCTCTCCCTATATTTTTTTTTTCGCGACTTACGTCATCTGGCAAATc                     >  1:800866/1‑150 (MQ=255)                                                                                                                               gTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATgg                 <  2:1138774/149‑1 (MQ=255)                                                                                                                               gTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGgc                <  1:1061105/150‑1 (MQ=255)                                                                                                                                 cGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGgcgc              <  1:1026500/150‑1 (MQ=255)                                                                                                                                     tATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGtt          <  1:1245514/150‑1 (MQ=255)                                                                                                                                        ggCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttatt      <  2:262530/151‑1 (MQ=255)                                                                                                                                        ggCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttat       <  2:261205/150‑1 (MQ=255)                                                                                                                                         gCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttatta     <  1:1374443/151‑1 (MQ=255)                                                                                                                                            gTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCa  <  2:1021323/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCA  >  NC_000913/3374976‑3375265

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