Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,931,372	G→A	6.3%	G164G (GGG→GGA)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,372	0	G	A	6.3%	194.7 / 3.4	63	G164G (GGG→GGA)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base A (2/2);  total (29/34)
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.48e-01

GGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGG  >  NC_000913/1931226‑1931514                                                                                                                                                   |                                                                                                                                               ggTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAg                                                                                                                                             <  2:1011457/150‑1 (MQ=255) ggTATCCGTCGCCTGGCGGCAGAAGAGCGGCAGCTGCCCACTTCCACTTAGCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAg                                                                                                                                             <  2:1011628/150‑1 (MQ=255)  gTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGa                                                                                                                                            <  2:131990/150‑1 (MQ=255)       cGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGttt                                                                                                                                       >  1:1518652/1‑150 (MQ=255)       cGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGttt                                                                                                                                       >  1:1512444/1‑150 (MQ=255)          cGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGACAGACGTTTAtt                                                                                                                                    >  2:273279/1‑150 (MQ=255)          cGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGACAGACGTTTAtt                                                                                                                                    >  2:279262/1‑150 (MQ=255)           gCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTc                                                                                                                                   <  1:834315/150‑1 (MQ=255)           gCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTc                                                                                                                                   <  1:1267939/150‑1 (MQ=255)              tGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGtt                                                                                                                                >  1:1393432/1‑150 (MQ=255)                gcggcAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCt                                                                                                                              <  2:1552095/150‑1 (MQ=255)                     agaagaGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCaga                                                                                                                         <  2:389781/150‑1 (MQ=255)                        agagCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCa                                                                                                                      <  1:399455/150‑1 (MQ=255)                              gcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGCGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGc                                                                                                                >  2:560368/1‑150 (MQ=255)                              gcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGc                                                                                                                <  2:518928/150‑1 (MQ=255)                              gcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGACAGACGTTTATTCGTTCTGCAGAGCAACTTGc                                                                                                                <  1:273279/150‑1 (MQ=255)                              gcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGACAGACGTTTATTCGTTCTGCAGAGCAACTTGc                                                                                                                <  1:279262/150‑1 (MQ=255)                                   tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAgg                                                                                                           <  2:935772/150‑1 (MQ=255)                                    gcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGc                                                                                                          <  1:1574169/150‑1 (MQ=255)                                    gcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGc                                                                                                          <  2:604496/150‑1 (MQ=255)                                       cacttccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATg                                                                                                       <  1:1121086/150‑1 (MQ=255)                                            ccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGt                                                                                                  <  1:1313456/150‑1 (MQ=255)                                             cacttATCGTTTTCCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTa                                                                                                 <  2:1503618/150‑1 (MQ=255)                                                 tATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCt                                                                                             >  2:1282708/1‑150 (MQ=255)                                                  aTCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTc                                                                                            >  1:183717/1‑150 (MQ=255)                                                   tCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCa                                                                                           >  2:1249986/1‑150 (MQ=255)                                                     gTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAGGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGc                                                                                         <  1:104559/150‑1 (MQ=255)                                                     gTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGc                                                                                         <  2:1006950/150‑1 (MQ=255)                                                          gCCGATAGCTAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGc                                                                                    <  2:1548102/150‑1 (MQ=255)                                                          gCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGc                                                                                         <  1:111011/145‑1 (MQ=255)                                                          gCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGc                                                                                         >  2:111011/1‑145 (MQ=255)                                                          gCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGc                                                                                    <  1:1619173/150‑1 (MQ=255)                                                           ccGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCg                                                                                   >  2:1617175/1‑150 (MQ=255)                                                            cGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCgg                                                                                  <  2:506544/150‑1 (MQ=255)                                                                 gCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTcgcg                                                                             >  2:1251422/1‑150 (MQ=255)                                                                 gCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGAAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTcgcg                                                                             >  1:235686/1‑150 (MQ=255)                                                                  cGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCgcgc                                                                            <  1:595158/150‑1 (MQ=255)                                                                     aaGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCgg                                                                         <  2:1393432/150‑1 (MQ=255)                                                                         cTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCg                                                                     >  2:177509/1‑150 (MQ=255)                                                                          ttttCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCgg                                                                    >  1:245062/1‑150 (MQ=255)                                                                            ttCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGc                                                                  <  2:1518652/150‑1 (MQ=255)                                                                            ttCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGc                                                                  <  2:1512444/150‑1 (MQ=255)                                                                               cgcgAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgc                                                               >  1:547294/1‑150 (MQ=255)                                                                                  gAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTa                                                            <  2:345538/150‑1 (MQ=255)                                                                                              gACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGc                                                >  2:1407522/1‑150 (MQ=255)                                                                                               aCATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcg                                               <  1:1282708/150‑1 (MQ=255)                                                                                                 aTTGGCTATCCCTGCATTGTAAAACGGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcgtc                                             >  2:209962/1‑150 (MQ=255)                                                                                                       tATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCTGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAg                                       <  2:1465832/150‑1 (MQ=255)                                                                                                       tATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAg                                       <  2:392302/150‑1 (MQ=255)                                                                                                        aTCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGt                                      >  2:925897/1‑150 (MQ=255)                                                                                                             tGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGAc                                 <  1:1249986/150‑1 (MQ=255)                                                                                                                 ttGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCg                             >  1:1544214/1‑150 (MQ=255)                                                                                                                   gTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGaa                           <  2:323846/150‑1 (MQ=255)                                                                                                                   gTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGaa                           >  1:58221/1‑150 (MQ=255)                                                                                                                      aaaCCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAAtt                        >  2:778938/1‑150 (MQ=255)                                                                                                                                  aGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAAcc            >  2:805501/1‑150 (MQ=255)                                                                                                                                   gCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCg           >  1:976737/1‑150 (MQ=255)                                                                                                                                     tctTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTc         >  1:545057/1‑150 (MQ=255)                                                                                                                                      ctTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCa        <  2:336090/150‑1 (MQ=255)                                                                                                                                      ctTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCa        >  2:46273/1‑150 (MQ=255)                                                                                                                                      ctTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCa        >  2:1054739/1‑150 (MQ=255)                                                                                                                                      ctTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCa        >  2:545539/1‑150 (MQ=255)                                                                                                                                            gCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCgg  >  2:602786/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                               GGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGG  >  NC_000913/1931226‑1931514

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