Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	103,582	A→T	6.4%	N143I (AAT→ATT)	ftsQ →	cell division protein FtsQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	103,582	0	A	T	6.4%	197.3 / 3.3	62	N143I (AAT→ATT)	ftsQ	cell division protein FtsQ
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base T (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGA  >  NC_000913/103434‑103731                                                                                                                                                     |                                                                                                                                                      atcatcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAAt                                                                                                                                                      >  1:445568/1‑150 (MQ=255)  tcatcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATa                                                                                                                                                     <  2:5900/150‑1 (MQ=255)   catcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATAc                                                                                                                                                    <  2:1315011/150‑1 (MQ=255)    atcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATAcc                                                                                                                                                   <  1:552943/150‑1 (MQ=255)     tcCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCt                                                                                                                                                  >  1:629526/1‑150 (MQ=255)       cAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTc                                                                                                                                                >  1:1291386/1‑150 (MQ=255)           cGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCTGAAGGTATTTCCTTCcgcg                                                                                                                                            <  2:528739/150‑1 (MQ=255)               aaTAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCTGAAGGTATTTCCTTCCGCGtcta                                                                                                                                        <  1:62461/150‑4 (MQ=255)                aTAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTgccg                                                                                                                                       <  2:445568/150‑1 (MQ=255)                        acGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACg                                                                                                                               <  1:293138/150‑1 (MQ=255)                           ccTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCAc                                                                                                                            >  2:257981/1‑150 (MQ=255)                                cGTGGATTAAGCAGGTGATCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagca                                                                                                                       <  2:396394/150‑1 (MQ=255)                                    gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                                   >  2:32550/1‑150 (MQ=255)                                    gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                                   >  1:753512/1‑150 (MQ=255)                                    gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                                   <  1:438448/150‑1 (MQ=255)                                    gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                                   <  2:254565/150‑1 (MQ=255)                                      ttAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAgg                                                                                                                 >  2:214333/1‑150 (MQ=255)                                        aaGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTg                                                                                                               <  2:1291386/150‑1 (MQ=255)                                          gCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCt                                                                                                             <  2:629526/150‑1 (MQ=255)                                             ggTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                                                                                          <  2:1175797/150‑1 (MQ=255)                                               tGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                                                                                        >  2:533521/1‑150 (MQ=255)                                                gAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAAt                                                                                                       <  2:862877/150‑1 (MQ=255)                                                  gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                                     <  2:393271/150‑1 (MQ=255)                                                   cGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCt                                                                                                    >  1:1034274/1‑150 (MQ=255)                                                    gTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGAAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTg                                                                                                   <  2:17860/150‑1 (MQ=255)                                                     tCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGt                                                                                                  >  2:711002/1‑150 (MQ=255)                                                      cAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCTTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTa                                                                                                 >  2:901004/1‑150 (MQ=255)                                                                 ggCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAg                                                                                      >  2:161978/1‑150 (MQ=255)                                                                  gCCTTATGAATTAAAGATTAATCTGGTTGAATATGTGAAGATTGCGCGGTGGAATGATCAACATATGGTGGAATCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAgg                                                                                     >  2:838738/1‑150 (MQ=255)                                                                     tgatgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAg                                                                                  >  1:1132967/1‑150 (MQ=255)                                                                      gatgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAAAATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCAGAAGGCAgc                                                                                 >  2:1041798/1‑150 (MQ=255)                                                                         gaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGcc                                                                              <  1:533521/150‑1 (MQ=255)                                                                         gaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGcc                                                                              <  2:1327797/150‑1 (MQ=255)                                                                           aTTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCaa                                                                            >  2:683394/1‑150 (MQ=255)                                                                                aGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAg                                                                       <  2:346390/150‑1 (MQ=255)                                                                                 gATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt                                                                      <  2:906022/150‑1 (MQ=255)                                                                                 gATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt                                                                      <  1:683394/150‑1 (MQ=255)                                                                                        ctcgTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAACGAAATACCTTCAGCGTGCCGCCACAACGCCCAAGCAAGCAGGTGCTTCCAATCCTGTATGGCCCGGAAGGCACCGCCAATCAAGTGTTGCAg                                                               >  2:112772/4‑150 (MQ=255)                                                                                          ggTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCgg                                                                                         <  2:500839/122‑1 (MQ=255)                                                                                          ggTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCgg                                                                                         >  1:500839/1‑122 (MQ=255)                                                                                             tGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTa                                                          <  1:711002/150‑1 (MQ=255)                                                                                               aaTATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATc                                                        >  1:71289/1‑150 (MQ=255)                                                                                                      gCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAAt                                                 >  1:713678/1‑150 (MQ=255)                                                                                                      gCCGATTGCGCGGTGGAATGATCAACATATCGTAGACGCGGAAGTAAATACATAAAGCGTGACGCCATAACGCACCAGCAAGCAGGTGCTTACAATGCTGTATTGCCCGGAAGCAAGCGCAAATGAAGTGTTGCAGGGCTATAGCGAAAt                                                 >  2:1394561/1‑150 (MQ=255)                                                                                                       ccGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATg                                                <  2:1368923/150‑1 (MQ=255)                                                                                                        cGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATgg                                               <  1:214333/150‑1 (MQ=255)                                                                                                           ttGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                            <  2:753512/150‑1 (MQ=255)                                                                                                                         gATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAg                              >  2:1402160/1‑150 (MQ=255)                                                                                                                          aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                            >  1:1051473/1‑135 (MQ=255)                                                                                                                          aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                            <  2:1051473/135‑1 (MQ=255)                                                                                                                          aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                             >  1:212747/1‑150 (MQ=255)                                                                                                                          aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                             >  2:1181573/1‑150 (MQ=255)                                                                                                                          aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                             >  1:1268056/1‑150 (MQ=255)                                                                                                                             aaCATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACa                          >  1:1316080/1‑150 (MQ=255)                                                                                                                               cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                             <  2:445023/145‑1 (MQ=255)                                                                                                                               cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                             >  1:445023/1‑145 (MQ=255)                                                                                                                                   ggcggcaCGCTGAAGGTATTTCCTTCCGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTa                    >  1:1353551/7‑150 (MQ=255)                                                                                                                                     gTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTAct                  >  1:324085/1‑150 (MQ=255)                                                                                                                                      tAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActc                 <  2:270464/150‑1 (MQ=255)                                                                                                                                       gcaCGCTGAAGGTATTTCCTTCCGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActct                >  1:195432/3‑150 (MQ=255)                                                                                                                                          cgcgGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTgaa             >  2:861853/1‑150 (MQ=255)                                                                                                                                                aggaaATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAgc       >  2:1427644/1‑150 (MQ=255)                                                                                                                                                     aTACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGa  <  1:310934/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGA  >  NC_000913/103434‑103731

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