Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	103,580	A→T	6.5%	G142G (GGA→GGT)	ftsQ →	cell division protein FtsQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	103,580	0	A	T	6.5%	194.0 / 2.2	61	G142G (GGA→GGT)	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 = 1.00e+00

AACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGC  >  NC_000913/103431‑103726                                                                                                                                                      |                                                                                                                                                   aaCATCATCCAGACGCAAATAGATCAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaagga                                                                                                                                                    <  2:997994/150‑1 (MQ=255) aaCATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaagga                                                                                                                                                    >  1:1079982/1‑150 (MQ=255)  aCATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCggaaggaa                                                                                                                                                   >  1:1315011/1‑150 (MQ=255)    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:254565/150‑1 (MQ=255)                                       gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                              <  1:438448/150‑1 (MQ=255)                                       gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                              >  2:32550/1‑150 (MQ=255)                                       gATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagcaagca                                                                                                              >  1:753512/1‑150 (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:1268056/1‑150 (MQ=255)                                                                                                                             aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                        >  1:212747/1‑150 (MQ=255)                                                                                                                             aTCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAgg                        >  2:1181573/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)                                                                                                                                                      |                                                                                                                                                   AACATCATCCAGACGCAAATAGAACAACGCCTGCCGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGC  >  NC_000913/103431‑103726

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