Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,219	T→A	6.3%	I385F (ATT→TTT) ‡	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,219	0	T	A	6.3%	208.2 / 2.3	73	I385F (ATT→TTT) ‡	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base T (40/28);  new base A (2/3);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGC  >  NC_000913/3728070‑3728368                                                                                                                                                      |                                                                                                                                                      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)                                                                                                                                                   tAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCt     >  2:1264290/1‑150 (MQ=255)                                                                                                                                                   tAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCt     >  1:1556537/1‑150 (MQ=255)                                                                                                                                                    aaTACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAATCCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTc    >  2:442246/1‑150 (MQ=255)                                                                                                                                                    aaTACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTc    >  2:752818/1‑150 (MQ=255)                                                                                                                                                     aTACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGAAGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCg   >  2:1310218/1‑150 (MQ=255)                                                                                                                                                      tACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGc  >  1:1295500/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGC  >  NC_000913/3728070‑3728368

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