Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	103,584	A→T	6.6%	T144S (ACC→TCC)	ftsQ →	cell division protein FtsQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	103,584	0	A	T	6.6%	189.9 / 3.4	61	T144S (ACC→TCC)	ftsQ	cell division protein FtsQ
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base T (2/2);  total (32/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.93e-01

TCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGA  >  NC_000913/103435‑103731                                                                                                                                                      |                                                                                                                                                    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                                                                                                                   <  2:254565/150‑1 (MQ=255)                                   gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                                   >  1:753512/1‑150 (MQ=255)                                   gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                                   <  1:438448/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)                                                                                                                                                      |                                                                                                                                                    TCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGA  >  NC_000913/103435‑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: ‑