Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,323,499	G→T	7.6%	P191Q (CCG→CAG)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,323,499	0	G	T	7.6%	201.2 / 4.3	61	P191Q (CCG→CAG)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base T (2/3);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCC  >  NC_000913/3323350‑3323645                                                                                                                                                      |                                                                                                                                                   cGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCg                                                                                                                                                    >  2:565272/1‑150 (MQ=255)  gAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCgg                                                                                                                                                   >  1:1007069/1‑150 (MQ=255)         aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaac                                                                                                                                            <  1:185110/150‑3 (MQ=255)         aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaac                                                                                                                                            <  1:1370053/150‑3 (MQ=255)           aCCGAGATCCGCTTTTTCCGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtgtg                                                                                                                                          <  1:1002448/150‑1 (MQ=255)            ccGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaactgg                                                                                                                                         <  1:197759/150‑6 (MQ=255)             cGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGAt                                                                                                                                        <  2:235152/150‑1 (MQ=255)                gaTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAg                                                                                                                                     <  1:904276/150‑1 (MQ=255)                gaTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAg                                                                                                                                     >  2:662842/1‑150 (MQ=255)                 aTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGt                                                                                                                                    >  1:1410053/1‑150 (MQ=255)                  tCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTc                                                                                                                                   >  1:768470/1‑150 (MQ=255)                           tCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGttt                                                                                                                          <  2:1007069/150‑1 (MQ=255)                               ccagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGcac                                                                                                                      <  1:5726/150‑1 (MQ=255)                                   cacaCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATc                                                                                                                  <  1:273887/150‑1 (MQ=255)                                    acacGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGCTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATcc                                                                                                                 <  1:995250/150‑1 (MQ=255)                                    acacGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGCTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATcc                                                                                                                 <  1:1000196/150‑1 (MQ=255)                                     cacGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGc                                                                                                                        <  2:1295176/142‑1 (MQ=255)                                     cacGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGc                                                                                                                        >  1:1295176/1‑142 (MQ=255)                                       cGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATccacc                                                                                                              >  2:654001/1‑150 (MQ=255)                                       cGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATccacc                                                                                                              <  1:390891/150‑1 (MQ=255)                                        gAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtg                                                                                                                                            >  1:460262/1‑119 (MQ=255)                                        gAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtg                                                                                                                                            <  2:460262/119‑1 (MQ=255)                                            cTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAAt                                                                                                         >  1:499637/1‑150 (MQ=255)                                              ggtgcgcgcGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCt                                                                                                       <  1:1213006/147‑1 (MQ=255)                                              ggAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCt                                                                                                       <  2:1046992/150‑1 (MQ=255)                                                aGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                                     <  2:1390705/150‑1 (MQ=255)                                                       gAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCAc                                                                                              <  1:193594/150‑1 (MQ=255)                                                             cGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGc                                                                                        <  2:1291554/150‑1 (MQ=255)                                                                   cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgtt                                                                                  >  1:1163895/1‑150 (MQ=255)                                                                    ccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttc                                                                                 >  1:1428833/1‑150 (MQ=255)                                                                      ccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACGATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcgt                                                                               >  1:865672/1‑150 (MQ=255)                                                                        tcTTCGGCATTGATGTTTACACCGTTTGGGCGAGAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCAGCGTTTGCACAATGCTGCACAATCTTCAGTTCTCTGAGGCTAAgttcgttc                                                                             <  2:304258/149‑1 (MQ=255)                                                                        aCTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCTCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcgttc                                                                             <  1:118417/150‑1 (MQ=255)                                                                           tCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGAGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCggg                                                                          <  1:763244/150‑1 (MQ=255)                                                                               cATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACg                                                                      <  2:1009498/150‑1 (MQ=255)                                                                                 ttGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTg                                                                    <  1:1102256/150‑1 (MQ=255)                                                                                         tACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCa                                                            >  2:1299317/1‑150 (MQ=255)                                                                                               gTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTc                                                      >  2:1281619/1‑150 (MQ=255)                                                                                                  tGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGat                                                   <  2:723028/150‑1 (MQ=255)                                                                                                    gCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGatat                                                 >  1:624441/1‑150 (MQ=255)                                                                                                     cTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGAtata                                                >  1:1079387/1‑150 (MQ=255)                                                                                                      tCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACGATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAg                                               <  2:865672/150‑1 (MQ=255)                                                                                                      tCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAg                                               >  1:154096/1‑150 (MQ=255)                                                                                                       cacaACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGc                                              >  2:912432/1‑150 (MQ=255)                                                                                                       cacaACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGc                                              >  1:430160/1‑150 (MQ=255)                                                                                                               aTAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCcgcg                                      <  2:1024287/150‑1 (MQ=255)                                                                                                                  gCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCa                                   <  1:1281619/150‑1 (MQ=255)                                                                                                                  gCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCa                                   >  1:1374109/1‑150 (MQ=255)                                                                                                                  gCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCa                                   >  1:28930/1‑150 (MQ=255)                                                                                                                       aaCGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAAc                              >  2:544160/1‑150 (MQ=255)                                                                                                                              gCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACgg                       >  2:348916/1‑150 (MQ=255)                                                                                                                               cccccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggc                      >  1:879018/1‑150 (MQ=255)                                                                                                                               cccccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggc                      <  1:662842/150‑1 (MQ=255)                                                                                                                                 cccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctg                    >  2:857979/1‑150 (MQ=255)                                                                                                                                   cAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctggc                  <  1:565272/150‑1 (MQ=255)                                                                                                                                    acatCGCGCCGAACGTGCTGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctggct                 >  1:510716/4‑150 (MQ=255)                                                                                                                                     catCGCGCCGAACGTGCTGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCtt                >  1:1127329/3‑150 (MQ=255)                                                                                                                                        cgcgCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTAc             >  1:628334/1‑150 (MQ=255)                                                                                                                                        cgcgCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTAc             >  2:672195/1‑150 (MQ=255)                                                                                                                                              gcaCGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCGACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGtt       <  2:768470/150‑1 (MQ=255)                                                                                                                                                  gTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGc   >  1:670438/1‑150 (MQ=255)                                                                                                                                                   ttCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGcc  <  2:430160/150‑1 (MQ=255)                                                                                                                                                   ttCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTGCTGcc  >  1:1074496/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCC  >  NC_000913/3323350‑3323645

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