Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,937	A→T	5.1%	I99F (ATC→TTC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,937	0	A	T	5.1%	241.7 / 3.8	79	I99F (ATC→TTC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base T (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

AAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGT  >  NC_002947/1549793‑1550086                                                                                                                                                 |                                                                                                                                                       aaGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGtt                                                                                                                                                   >  1:1746996/1‑149 (MQ=255)  aGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCa                                                                                                                                                  <  1:1294504/150‑1 (MQ=255)        gCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCCTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                                                                                                            <  2:1847831/150‑1 (MQ=255)               cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGag                                                                                                                                     <  2:1746996/150‑2 (MQ=255)                gccggcAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCGGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCATGCGATCGACAACGAAGGCCTGaaa                                                                                                                                    >  2:967046/6‑150 (MQ=255)                gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                    >  1:1413267/1‑150 (MQ=255)                gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                    <  1:784864/150‑1 (MQ=255)                      aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCTCCCAGTTGGTCAAGGAGGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCAt                                                                                                                                                       <  1:1770674/125‑1 (MQ=255)                      aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCAt                                                                                                                                                       >  2:1770674/1‑125 (MQ=255)                      aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                              >  1:1817851/1‑150 (MQ=255)                      aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                              >  2:1805454/1‑150 (MQ=255)                        aGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGc                                                                                                                            >  1:1794837/1‑150 (MQ=255)                           aaTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGc                                                                                                                         >  2:650394/1‑150 (MQ=255)                            aTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGcc                                                                                                                        >  1:324720/1‑150 (MQ=255)                                 gCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCAt                                                                                                                   <  2:1794837/150‑1 (MQ=255)                                 gCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCAt                                                                                                                   >  2:52794/1‑150 (MQ=255)                                      aaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAAcc                                                                                                              <  1:1263762/150‑1 (MQ=255)                                        gACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCg                                                                                                            <  1:650394/150‑1 (MQ=255)                                          cGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAt                                                                                                          <  1:530579/150‑1 (MQ=255)                                           gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑g                                                                                                        <  1:263348/150‑1 (MQ=255)                                           gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑g                                                                                                        <  2:1413267/150‑1 (MQ=255)                                            cTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAttg                                                                                                        <  2:1662640/150‑3 (MQ=255)                                             tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGa                                                                                                      >  1:47859/1‑150 (MQ=255)                                              ttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGAc                                                                                                     <  1:1818370/150‑1 (MQ=255)                                              ttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGAc                                                                                                     <  2:1530543/150‑1 (MQ=255)                                                 gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGagcc                                                                                                                                   >  2:365550/1‑115 (MQ=38)                                                 gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGagcc                                                                                                                                   <  1:365550/118‑4 (MQ=38)                                                 gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCaa                                                                                                                                                 <  1:1889345/104‑1 (MQ=255)                                                 gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCaa                                                                                                                                                 >  2:1889345/1‑104 (MQ=255)                                                 gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTg                                                                                                  >  2:218163/1‑150 (MQ=255)                                                 gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTg                                                                                                  >  1:1440783/1‑150 (MQ=255)                                                  aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGa                                                                                                 >  2:461705/1‑150 (MQ=255)                                                    aaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAg                                                                                               >  2:226373/1‑150 (MQ=255)                                                            cgcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTAt                                                                                       <  1:1309203/150‑1 (MQ=255)                                                              cccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCg                                                                                     >  2:210411/1‑150 (MQ=255)                                                                  gCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACaa                                                                                 <  1:1805454/150‑1 (MQ=255)                                                                    tGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAgg                                                                               <  1:210411/150‑1 (MQ=255)                                                                     ggTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGc                                                                              <  1:132429/150‑1 (MQ=255)                                                                     ggTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACAACCCCCGCTACCGTTCTGGCGCAGGCCACCCTCAACGAAGGCGTGCAAGCGGTCGCTGCCGGCAGGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGc                                                                              <  2:216685/150‑1 (MQ=255)                                                                         aaGGAAGTTGCTTCCAAGGCCAACCAGGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCAct                                                                          >  1:450334/1‑150 (MQ=255)                                                                          aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCActg                                                                         >  2:201903/1‑150 (MQ=255)                                                                          aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCActg                                                                         >  1:1759167/1‑150 (MQ=255)                                                                           ggAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCActgc                                                                        >  1:1879793/1‑150 (MQ=255)                                                                                ttGCTTCCGAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcg                                                                   >  2:1497275/1‑150 (MQ=255)                                                                                ttGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                                                                                                          >  1:1259530/1‑80 (MQ=255)                                                                                ttGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                                                                                                          <  2:1259530/80‑1 (MQ=255)                                                                                   cTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcg                                                                >  1:1355641/1‑150 (MQ=255)                                                                                    ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                           <  2:1209099/107‑1 (MQ=255)                                                                                    ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                           >  1:1209099/1‑107 (MQ=255)                                                                                    ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGcgcg                                                                                           >  1:111618/1‑122 (MQ=255)                                                                                    ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGcgcg                                                                                           <  2:111618/122‑1 (MQ=255)                                                                                      ccAAGGCCAACGACGCTGCCGGTGACGGTACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                             >  2:425068/1‑150 (MQ=255)                                                                                        aaGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                               >  2:954405/1‑146 (MQ=255)                                                                                        aaGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                               <  1:954405/146‑1 (MQ=255)                                                                                            ccAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAAc                                                       <  1:226373/150‑1 (MQ=255)                                                                                                cgacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaa                                                   <  2:1759167/150‑1 (MQ=255)                                                                                                ccacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaa                                                   <  1:791040/148‑1 (MQ=255)                                                                                                 gacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaag                                                  >  2:1506929/1‑150 (MQ=255)                                                                                                      tGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGAc                                                                                                     >  2:517205/1‑94 (MQ=255)                                                                                                      tGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGAc                                                                                                     <  1:517205/94‑1 (MQ=255)                                                                                                      tGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCt                                             >  1:1612366/1‑150 (MQ=255)                                                                                                      tGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCt                                             <  1:52794/150‑1 (MQ=255)                                                                                                       gCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcg                                                                >  2:9691/1‑130 (MQ=255)                                                                                                       gCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcg                                                                <  1:9691/130‑1 (MQ=255)                                                                                                            tGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCCTCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCCGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                       <  1:1842690/150‑1 (MQ=255)                                                                                                              aCGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGc                                     <  1:1506929/150‑1 (MQ=255)                                                                                                                        caccGCTACCGTTCTGGCTCAGGCCATCGTCAACGACGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTGCCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCACCGa                           >  2:359859/1‑150 (MQ=255)                                                                                                                        caccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGa                           >  2:706339/1‑150 (MQ=255)                                                                                                                             cTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCa                      >  2:1134933/1‑150 (MQ=255)                                                                                                                               aCCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGa                           <  2:145956/143‑1 (MQ=255)                                                                                                                               aCCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGa                           >  1:145956/1‑143 (MQ=255)                                                                                                                                   ttCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCAAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCa                >  1:1617134/1‑150 (MQ=255)                                                                                                                                     cTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATc              <  1:218163/150‑1 (MQ=255)                                                                                                                                     cTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAA‑CCGATCTGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATc              <  1:19940/150‑1 (MQ=255)                                                                                                                                         cTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGccc          >  2:1807423/1‑150 (MQ=255)                                                                                                                                          tCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTGCAAGCCATGCGCCGACTCCAAGGCCATCGCCCa         >  1:10142/1‑150 (MQ=255)                                                                                                                                               ccATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtag    <  1:706339/150‑1 (MQ=255)                                                                                                                                                cATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                             <  2:332210/92‑1 (MQ=255)                                                                                                                                                cATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                             >  1:332210/1‑92 (MQ=255)                                                                                                                                                 aTCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtaggt  >  2:711241/1‑150 (MQ=255)                                                                                                                                                 aTCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTAGCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggaacgg  >  2:921854/1‑145 (MQ=255)                                                                                                                                                 |                                                                                                                                                       AAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAT‑GGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGT  >  NC_002947/1549793‑1550086

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