Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,215	T→A	6.8%	K386I (AAA→ATA)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,215	0	T	A	6.8%	178.2 / 2.9	68	K386I (AAA→ATA)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base T (34/28);  new base A (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACG  >  NC_000913/3728066‑3728363                                                                                                                                                      |                                                                                                                                                     ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                                                                                                                      <  2:1588961/150‑1 (MQ=255) ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                                                                                                                      <  1:1306954/150‑1 (MQ=255)     gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                                                                                                                  <  2:1051594/150‑1 (MQ=255)       tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc                                                                                                                                                <  2:1096164/150‑1 (MQ=255)        cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAcc                                                                                                                                               <  1:139892/150‑1 (MQ=255)              cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGc                                                                                                                                         <  1:1411825/150‑1 (MQ=255)              cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaag                                                                                                                                         <  2:999183/150‑4 (MQ=255)               aGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaagt                                                                                                                                        <  1:1383947/150‑5 (MQ=255)                gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCTATCAACGTAACACTTTTCGGTATTAAACCGCaaagtt                                                                                                                                       <  1:1477103/150‑6 (MQ=255)                          gcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATc                                                                                                                             <  2:426496/150‑1 (MQ=255)                                  cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt                                                                                                                     <  2:1193967/150‑1 (MQ=255)                                         cacaTCCCCCCCCGTACTGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAg                                                                                                              <  2:1310478/150‑1 (MQ=255)                                             tCCCCCCCCGTACGGTAGTCGAGCTGCTGACCGCGGGTGTCCGCCCGCCTCGGCCGCCCGTACTCCCGACGCGCCCCGCCCCCAATCAACGTAACACTTTGGGGTTTTATCCCGCAGCGCTGGACGCCATACATGCCAGCTGCCAgcggg                                                                                                          >  2:1557932/1‑148 (MQ=255)                                                ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCACCGTAACCCTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCagtc                                                                                                       >  1:687406/1‑147 (MQ=255)                                                  ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc                                                                                                     <  1:712258/150‑1 (MQ=255)                                                        aCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGcc                                                                                               <  2:188250/150‑1 (MQ=255)                                                             aaTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCa                                                                                          >  2:622674/1‑150 (MQ=255)                                                              aTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAg                                                                                         >  1:1560124/1‑150 (MQ=255)                                                                cGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCCGGCATGCACGACATCCGTGCCATCTGCCAGCGCATAACCCACGCCGTCCAGCa                                                                                       >  2:1497293/1‑150 (MQ=255)                                                                 gAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCAc                                                                                      >  2:1021779/1‑150 (MQ=255)                                                                   gctgctGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTg                                                                                    >  1:991269/1‑150 (MQ=255)                                                                         gACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACCTCCATGCCATCTGCCCGCGCATAACCCACGCCTTCCAGCACTGCTcgcg                                                                              >  1:58202/1‑150 (MQ=255)                                                                         gACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTcgcg                                                                              >  1:1260185/1‑150 (MQ=255)                                                                             gCTGATATCCGCCAGCATCTGACGCCAATACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAg                                                                          >  2:100486/1‑150 (MQ=255)                                                                              cTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACGCTTTGCGGTATACTACCGCAGGCATGCCCGACATCCATGCCATCTGCCCGCGCATAACCCACGCCTTCCAGCCCTGCGCGCGCCAgg                                                                         >  1:339358/1‑149 (MQ=255)                                                                              cTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGt                                                                         >  2:1649541/1‑150 (MQ=255)                                                                               tGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGtt                                                                        >  2:719831/1‑150 (MQ=255)                                                                                 atatCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCa                                                                      >  2:1613735/1‑150 (MQ=255)                                                                                 atatCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCa                                                                      >  2:1558724/1‑150 (MQ=255)                                                                                  tatCCGCCAGTATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGGTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAt                                                                     <  1:313210/150‑1 (MQ=255)                                                                                  tatCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCCTCACCCACGCCTTCCCGCACTGCTCGCGCCAGTTCAt                                                                     >  1:547312/1‑150 (MQ=255)                                                                                  tatCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAt                                                                     >  1:554092/1‑150 (MQ=255)                                                                                   atCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAtt                                                                    >  1:1279143/1‑150 (MQ=255)                                                                                    tCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTg                                                                   >  2:1476937/1‑150 (MQ=255)                                                                                       gCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTgggg                                                                <  2:228540/150‑1 (MQ=255)                                                                                       gCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTgggg                                                                <  2:436407/150‑1 (MQ=255)                                                                                        ccAGCATCTGCCGCCAGTGCTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCTCGACATCCATGCCATCGGCCAGCGCATAACCCCCGCGTTCCAGCCCTGCTCGCGCCAGGTCATTGGGGc                                                               >  2:161216/1‑150 (MQ=255)                                                                                          aGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCa                                                             <  1:1021779/150‑1 (MQ=255)                                                                                             aTCTGACTCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGt                                                          <  2:1684360/150‑1 (MQ=255)                                                                                             aTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGt                                                          <  1:868300/150‑1 (MQ=255)                                                                                                    gCCAGGACTCACTCCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATTCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGGGCATTGGGGCCATGTTGATGAg                                                   <  2:1147876/150‑1 (MQ=255)                                                                                                      cAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCCGGCCAGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTc                                                 >  1:614180/1‑150 (MQ=255)                                                                                                        gTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCaa                                               >  1:391770/1‑150 (MQ=255)                                                                                                          gcTCACTACGCGCCCCGCCCCCAATCATCGTAACACTTTGCGGTGTAATACCGCAGGCATGCACGACTTCCATGCCATCTGCCAGCGCATAACCCACGCATTCCAGCACTTCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAAc                                             <  2:1553816/149‑1 (MQ=255)                                                                                                           cTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTACAACCGCAGGCATGCACGACATCCATGCCATCTGCCTGCGCATAACCCACGCCTTCCCGCACTGCTCGCGCCAGTTCCTTGGGGCCATGTTGATGAGTCCAAca                                            >  1:1553816/1‑149 (MQ=255)                                                                                                           cTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCCCGACATCCATGCCATCTGCCAGCGCATACCCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAAcc                                            >  1:1172823/1‑150 (MQ=255)                                                                                                           cTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCATTTCATTGGGGCCATGTTGATGAGTCAAAcc                                            <  1:1658413/150‑1 (MQ=255)                                                                                                              aCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGCCATGCACGACATCCATGCCATCTGCCAGCGCATACCCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaaa                                         >  1:36022/1‑150 (MQ=255)                                                                                                                 aCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGaa                                      <  1:200284/150‑1 (MQ=255)                                                                                                                 aCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTAAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGaa                                      <  1:205896/150‑1 (MQ=255)                                                                                                                    cgccccgcccCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAc                                   <  1:1380317/150‑1 (MQ=255)                                                                                                                     gccccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAcc                                  >  2:1422136/1‑150 (MQ=255)                                                                                                                          gccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCtt                             >  2:1236417/1‑150 (MQ=255)                                                                                                                           cccccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCTGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTc                            >  1:929783/1‑150 (MQ=255)                                                                                                                               cAATCAACGTAACGCTTTGCGGTTTAATACCGCAGGCAGGCACGACATCCATGCCATCTGCCAGCGAATAACCCCCGCCTTACAGCACTGCTCGCTCCAGTTCAGTGGGGCCAGGTTGATTAGGCAAACCAAAGAAATCCCCCTGCGccg                        >  2:883889/1‑149 (MQ=255)                                                                                                                               cAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCt                        <  2:1560124/150‑1 (MQ=255)                                                                                                                                aaTCCACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTg                       >  1:1292129/1‑150 (MQ=255)                                                                                                                                aaTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTg                       <  1:1558724/150‑1 (MQ=255)                                                                                                                                 aTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgg                      <  1:622674/150‑1 (MQ=255)                                                                                                                                 aTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgg                      <  1:1356210/150‑1 (MQ=255)                                                                                                                                   cAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgggg                    <  1:1641700/150‑1 (MQ=255)                                                                                                                                     aCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGat                  >  2:1715136/1‑150 (MQ=255)                                                                                                                                      cGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGACAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGatt                 >  2:1339766/1‑150 (MQ=255)                                                                                                                                       gTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATTCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGatta                <  1:1557932/150‑1 (MQ=255)                                                                                                                                        tAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGattat               >  1:958994/1‑150 (MQ=255)                                                                                                                                           cacTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgt            <  1:1613735/150‑1 (MQ=255)                                                                                                                                              tttGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgtg         <  1:810963/150‑1 (MQ=255)                                                                                                                                              cctGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCAGGTTGATGAGTCAACCCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgtg         >  2:316227/3‑150 (MQ=255)                                                                                                                                                   ggTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAACACCCCCTTCGCCTGGGGATTATTGTGTGGCGTa    >  2:922838/1‑150 (MQ=255)                                                                                                                                                    gTTTAATACCGCAGGCATGCACGACATCCATGCCATCTACCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTAc   >  2:1063804/1‑150 (MQ=255)                                                                                                                                                     tttAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACg  >  1:806793/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACG  >  NC_000913/3728066‑3728363

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