Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	668,263	G→T	10.2%	F87L (TTC→TTA) ‡	PP_0571 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	668,263	0	G	T	10.2%	254.8 / 17.2	88	F87L (TTC→TTA) ‡	PP_0571	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base T (5/4);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCT  >  NC_002947/668117‑668406                                                                                                                                                   |                                                                                                                                                gAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGaa                                                                                                                                               >  8:18315/1‑149 (MQ=255)        aaaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGCTCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGCTGATTCCGAGGTTGAGCGAGGCGCGCATGCATGcgcgt                                                                                                                                        >  8:161515/1‑148 (MQ=38)         aaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcc                                                                                                                                       <  7:327342/149‑3 (MQ=255)         aaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcc                                                                                                                                       <  2:108945/149‑3 (MQ=255)         aaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGa                                                                                                                                       >  3:93283/1‑149 (MQ=255)            ccTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcct                                                                                                                                      <  4:264248/147‑4 (MQ=255)            ccTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcct                                                                                                                                      >  3:264248/1‑144 (MQ=255)               agggCTGCACAGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                                  >  4:210161/2‑148 (MQ=255)                gggCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                                  <  2:112487/147‑1 (MQ=255)                  gctgctGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAg                                                                                                                               >  1:315743/1‑148 (MQ=255)                      cacagCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcag                                                                                                                          <  3:210161/145‑1 (MQ=255)                       tgtgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcag                                                                                                                          >  6:50313/1‑148 (MQ=255)                         tgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcagcag                                                                                                                       >  3:162523/1‑149 (MQ=255)                         tgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcagca                                                                                                                        <  4:303862/148‑1 (MQ=255)                            aGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCAt                                                                                                                    >  6:208810/1‑149 (MQ=255)                                   tCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa                                                                                                                                <  4:401308/130‑1 (MQ=255)                                   tCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa                                                                                                                                >  3:401308/1‑130 (MQ=255)                                   tCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCg                                                                                                             <  4:162523/149‑1 (MQ=255)                                    cgcgACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATTATTCCGAGGATGATCGAGGCGCGCATGAATTCGCGGATTGCCGTGTACAGCAGCATCTGGCCga                                                                                                            <  2:98776/149‑1 (MQ=255)                                     gcgACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGTTCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCga                                                                                                            <  1:272616/148‑1 (MQ=255)                                       gACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatga                                                                                                         >  3:349178/1‑149 (MQ=255)                                          gCAAGGCCGCTCCCACAAAGACAGCGCAGATCAGAGGATCACCCAGCTCATCACACACAGCACCAGCCCCATCCAGATCAATCCGAGGATGAACGAGGCGCGCATGAATGAGCGGATTGACGAGTACAGCAGCATCAGCCCgatgaac                                                                                                       >  2:201648/1‑146 (MQ=255)                                            aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa                                                                                                                                <  1:214880/121‑1 (MQ=255)                                            aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa                                                                                                                                >  2:214880/1‑121 (MQ=255)                                            aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCg                                                                                                             <  5:374768/140‑1 (MQ=255)                                            aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCg                                                                                                             >  6:374768/1‑140 (MQ=255)                                            aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGTTCAgg                                                                                                    >  4:183311/1‑149 (MQ=255)                                            aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAgg                                                                                                    >  6:22954/1‑149 (MQ=255)                                                   cTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGa                                                                                             <  5:402903/149‑1 (MQ=255)                                                      ccACAAAGCCCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATAAGGGCCAGGATGc                                                                                          >  8:214465/1‑149 (MQ=255)                                                      ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGAt                                                                                            <  1:27319/147‑1 (MQ=255)                                                      ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGAt                                                                                            <  8:144062/147‑1 (MQ=255)                                                      ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGAt                                                                                            <  2:315743/147‑1 (MQ=255)                                                        acaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCt                                                                                         >  4:357467/1‑148 (MQ=255)                                                        acaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCt                                                                                         >  6:193990/1‑148 (MQ=255)                                                          aaaGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatga                                                                                                         >  7:1511/1‑130 (MQ=255)                                                          aaaGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatga                                                                                                         <  8:1511/130‑1 (MQ=255)                                                          aaaGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTcag                                                                                      <  5:22954/149‑1 (MQ=255)                                                             gACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTcagcag                                                                                   >  1:73933/1‑149 (MQ=255)                                                                cgcgCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCg                                                                                 >  6:77513/1‑148 (MQ=255)                                                                     aGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAgtgt                                                                           <  5:193990/149‑1 (MQ=255)                                                                      gATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTc                                                                          <  4:93283/149‑1 (MQ=255)                                                                         cAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCAt                                                                       >  8:299545/1‑149 (MQ=255)                                                                          aGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATg                                                                      <  2:73933/149‑1 (MQ=255)                                                                                   cccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTg                                                             >  1:312303/1‑149 (MQ=255)                                                                                    ccAGCTCATCACCCACAGCCACAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCACGCCGATGATCAGGGCCAGGATGCTCGGCAGCGAAGTGTCCATGCCCAAGGTgc                                                            >  7:177642/1‑149 (MQ=255)                                                                                         tcatcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGa                                                        >  6:204465/1‑148 (MQ=255)                                                                                           atcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGTTCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACa                                                      <  3:183311/148‑1 (MQ=255)                                                                                              aCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGc                                                   >  5:52339/1‑148 (MQ=255)                                                                                              aCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGc                                                   >  5:275329/1‑148 (MQ=255)                                                                                                ccACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAg                                                                                                     >  7:117560/1‑96 (MQ=255)                                                                                                ccACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAg                                                                                                     <  8:117560/96‑1 (MQ=255)                                                                                                ccACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCAt                                                <  5:50313/149‑1 (MQ=255)                                                                                                  acaGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATc                                               >  3:268419/1‑148 (MQ=255)                                                                                                  acaGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCg                                              <  6:275329/149‑1 (MQ=255)                                                                                                  acaGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGAGCGACAGGACATc                                               >  6:314812/1‑148 (MQ=255)                                                                                                      ccccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACGGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATg                                           >  5:375031/1‑148 (MQ=255)                                                                                                      ccccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATg                                           >  6:145898/1‑148 (MQ=255)                                                                                                      ccccAGCACCAGCCACATGAAAACGAGGATAATCGAGGCCCCCATGAATGCCCGGATTGCCGAGAACAGCAGCATCAGGCCCATGAACAGGGCCAGCATGCTCAGCAGCGAAGAGTCAATGACCAAGCTGCGCGACAGGCCATCGATg                                           >  2:335358/1‑148 (MQ=255)                                                                                                       cacAGCACCAGCCAGATTATTCGGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACATCAGCATCAGGCCGATGATCTGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGaa                                         <  1:201648/147‑1 (MQ=255)                                                                                                        ccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCga                                                                                                            >  7:327150/1‑81 (MQ=255)                                                                                                        ccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCga                                                                                                            <  8:327150/81‑1 (MQ=255)                                                                                                         cAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGt                                       >  3:242687/1‑149 (MQ=255)                                                                                                          aGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGtt                                      <  2:155441/149‑1 (MQ=255)                                                                                                          aGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGtt                                      <  4:349178/149‑1 (MQ=255)                                                                                                              ccagccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCa                                  >  5:92922/1‑149 (MQ=255)                                                                                                              ccagccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCa                                  <  5:204465/149‑1 (MQ=255)                                                                                                                     gatgatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCAt                           <  6:52339/149‑1 (MQ=255)                                                                                                                     gatgatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGCGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCAt                           >  1:207347/1‑149 (MQ=255)                                                                                                                      atgatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCAtt                          <  3:357467/149‑1 (MQ=255)                                                                                                                      atgatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCAtt                          >  4:164409/1‑149 (MQ=255)                                                                                                                       tgatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTg                         <  2:312303/149‑1 (MQ=255)                                                                                                                          ttCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGAGGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGcc                      <  7:161515/149‑1 (MQ=255)                                                                                                                            ccGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCca                     <  4:242687/148‑1 (MQ=255)                                                                                                                                 gatgatCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGt               >  4:37514/1‑149 (MQ=255)                                                                                                                                   tgatCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTg             <  5:77513/149‑1 (MQ=255)                                                                                                                                     atCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGTGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGa            <  4:269716/148‑1 (MQ=255)                                                                                                                                      tCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaac          <  5:145898/149‑1 (MQ=255)                                                                                                                                       caatccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGcc                      <  5:199491/131‑1 (MQ=255)                                                                                                                                       caatccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGcc                      >  6:199491/6‑136 (MQ=255)                                                                                                                                        gAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaacaa        >  4:200415/1‑149 (MQ=255)                                                                                                                                        aatccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaacaa        >  7:250174/5‑149 (MQ=255)                                                                                                                                         aGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaacaac       <  4:268419/149‑1 (MQ=255)                                                                                                                                          tccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAAcc      >  3:341790/3‑149 (MQ=255)                                                                                                                                          tccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAAcc      >  5:305793/3‑149 (MQ=255)                                                                                                                                             gcgcATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGc   >  8:190343/1‑149 (MQ=255)                                                                                                                                             gcgcATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTCAACAACCCGc   >  3:245660/1‑149 (MQ=255)                                                                                                                                              cgcATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCt  <  3:37514/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCT  >  NC_002947/668117‑668406

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