Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,540,881	G→T	6.2%	C35* (TGC→TGA)	PP_4872 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,540,881	0	G	T	6.2%	250.6 / 3.6	78	C35* (TGC→TGA)	PP_4872	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base T (2/3);  total (35/43)
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

GCAGACGGCCAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTGGTGAG  >  NC_002947/5540742‑5541025                                                                                                                                            |                                                                                                                                                 gcagACGGCCAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACAGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGcc                                                                                                                                         >  4:238205/1‑149 (MQ=255)   agACGGCCAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGcc                                                                                                                                         >  2:282247/1‑147 (MQ=255)    gACGGCCAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAg                                                                                                                                      <  4:27030/149‑1 (MQ=255)       ggCCAATGATGGCCATGCCGAACGACGGATGGTAGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTAc                                                                                                                                   >  1:50531/1‑149 (MQ=255)         ccAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTACTGGGCCGAATGCTTCGGACCAGCAcct                                                                                                                                 <  6:155294/149‑3 (MQ=255)         ccAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTACCGGGCCGAATGCTTCGGACCAGCAcct                                                                                                                                 >  6:229906/1‑147 (MQ=255)         ccAAGGCTGGCCAGGCCGCACGACGGAGGGTCGACGTAGGCATCACCGGGCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACtt                                                                                                                                 <  4:271407/149‑1 (MQ=255)          cAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGGCCGAATGCTTCGGACCAGCAcctt                                                                                                                                <  5:26904/149‑4 (MQ=255)          cAATGATGGCCATGCCGAACGACGAATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTg                                                                                                                                <  4:210698/149‑1 (MQ=255)           aaTGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTgg                                                                                                                               <  1:51501/149‑1 (MQ=255)           aaTGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTgg                                                                                                                               <  5:385887/149‑1 (MQ=255)           aaTGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTgg                                                                                                                               >  4:13276/1‑149 (MQ=255)             tgatgGCCATGCAGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGa                                                                                                                             <  5:320684/149‑1 (MQ=255)               atgGCCATGCCGAACGACGGATGCTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTACTGGGCCGAATGCTTCGGACCAGCAcctttcct                                                                                                                            <  5:65715/148‑8 (MQ=38)                tgGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAg                                                                                                                          <  5:193084/149‑1 (MQ=255)                     aTGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTa                                                                                                                     <  3:13276/149‑1 (MQ=255)                     aTGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTa                                                                                                                     <  3:175748/149‑1 (MQ=255)                      tGCCGAACGACGTATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTag                                                                                                                    <  6:293073/149‑1 (MQ=255)                      tGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATACATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCg                                                                                                                       >  8:115250/1‑146 (MQ=255)                       gCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTa                                                                                                                     >  4:68705/1‑147 (MQ=255)                        ccGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATTGTCGTagag                                                                                                                  >  5:273011/1‑149 (MQ=255)                         cGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGt                                                                                                                 >  1:366710/1‑149 (MQ=255)                          gAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTg                                                                                                                >  2:154355/1‑149 (MQ=255)                          gAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTg                                                                                                                >  7:194149/1‑149 (MQ=255)                           aaCGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTgg                                                                                                               >  3:296234/1‑149 (MQ=255)                              gacgGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGt                                                                                                                                     <  7:18465/124‑1 (MQ=255)                              gacgGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGt                                                                                                                                     >  8:18465/1‑124 (MQ=255)                              gacgGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGttt                                                                                                            <  1:292687/149‑1 (MQ=255)                                cgGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTgg                                                                                                          >  3:38598/1‑149 (MQ=255)                                 ggATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGgc                                                                                                         <  6:384109/149‑1 (MQ=255)                                              aGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGt                                                                                              >  8:240092/1‑147 (MQ=255)                                               ggCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAg                                                                                           <  3:68705/149‑1 (MQ=255)                                                   tCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTTGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGAc                                                                                       <  3:238205/149‑1 (MQ=255)                                                   tCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGAc                                                                                       >  8:117738/1‑149 (MQ=255)                                                     aCCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCg                                                                                     <  7:240092/149‑1 (MQ=255)                                                        ggTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGtt                                                                                  <  8:208021/149‑1 (MQ=255)                                                            aCGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTg                                                                               <  4:38598/148‑1 (MQ=255)                                                              gatgatgatGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGgtg                                                                            <  4:256171/149‑1 (MQ=255)                                                                     atgTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAggg                                                                     <  2:50531/149‑1 (MQ=255)                                                                          gCAGGAATCCCCGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATc                                                                >  1:255000/1‑149 (MQ=255)                                                                           cAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTg                                                                               <  2:364944/133‑1 (MQ=255)                                                                           cAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTg                                                                               >  1:364944/1‑133 (MQ=255)                                                                            aGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTTCTGGTCCGAAGGATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGTTGTGCAGGGAAATcgc                                                              <  3:66730/149‑1 (MQ=255)                                                                            aGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATcgc                                                              <  7:117738/149‑1 (MQ=255)                                                                             ggAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATcgcg                                                             <  7:191105/149‑1 (MQ=255)                                                                             ggAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATcgcg                                                             >  8:191105/1‑149 (MQ=255)                                                                              gAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGa                                                                                                                             <  4:285484/84‑1 (MQ=255)                                                                              gAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGa                                                                                                                             >  3:285484/1‑84 (MQ=255)                                                                              gAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCgg                                                            <  3:347011/149‑1 (MQ=255)                                                                                  cccAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGgcgc                                                        >  5:10710/1‑149 (MQ=255)                                                                                   acAGCCGAGCAGATCAATCTCATCCCTGCTCATTGTCATGAAACGTGCTGGTCCGAATAATTAGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGAGTGCATGGAAATCGAgggagag                                                       >  2:18603/2‑145 (MQ=255)                                                                                     aGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTg                                                                                                                >  4:42484/1‑90 (MQ=255)                                                                                     aGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTg                                                                                                                <  3:42484/90‑1 (MQ=255)                                                                                      gCCGAGCAGATCCATCTCCTGCCTGCTCATTGGCAGGTAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGaa                                                    <  7:313310/149‑1 (MQ=255)                                                                                       ccGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAat                                                   <  2:366710/149‑1 (MQ=255)                                                                                         gAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGACCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGAATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAatat                                                 >  4:112791/1‑149 (MQ=255)                                                                                              gATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACATGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTGCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGca                                             >  5:68760/1‑148 (MQ=255)                                                                                               aTCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAACCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACAGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGcag                                            >  5:138097/1‑147 (MQ=255)                                                                                                  cATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACaaa                                         <  1:282247/148‑1 (MQ=255)                                                                                                   aTCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAAtt                                       <  2:255000/149‑1 (MQ=255)                                                                                                        ctccCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTg                                    <  6:149756/147‑1 (MQ=255)                                                                                                        ctccCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTg                                    >  5:149756/1‑147 (MQ=255)                                                                                                        ctccCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGAc                                  >  8:35279/1‑149 (MQ=255)                                                                                                        ctccCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGAc                                  <  4:296234/149‑1 (MQ=255)                                                                                                               ccaatgagCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAAt                             <  4:155724/140‑1 (MQ=255)                                                                                                               cTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATc                            <  7:34670/148‑1 (MQ=255)                                                                                                                tCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCg                          >  8:79659/1‑149 (MQ=255)                                                                                                                  aTTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGGCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGa                         <  6:110706/148‑1 (MQ=255)                                                                                                                     ggCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTg                      <  8:194149/148‑1 (MQ=255)                                                                                                                           aaaGGTGCTGGTGCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGt               <  1:18603/149‑1 (MQ=255)                                                                                                                           aaaGGTGCTGGTCCGAAGCATTAGGCCCAGTACTTGGGATATTCGTAGAGTGGATTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATAGCGGGCGCGGAATATAGAACAAATTATGACAAAATCCGa                         >  7:70011/1‑139 (MQ=255)                                                                                                                            aaGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTg              <  7:115250/149‑1 (MQ=255)                                                                                                                             aGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTgg             <  6:10710/149‑1 (MQ=255)                                                                                                                                    ggTCCGAAGCATTCGGCCCAGTACTTGGAATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTgg      >  3:303477/1‑149 (MQ=255)                                                                                                                                     gTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTGGt     >  3:248518/1‑149 (MQ=255)                                                                                                                                     gTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAACACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTGGt     >  5:213379/1‑149 (MQ=255)                                                                                                                                         gAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTGGTgag  >  1:346395/1‑148 (MQ=255)                                                                                                                                         gAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGACAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTGGTgag  >  6:292779/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                                 GCAGACGGCCAATGATGGCCATGCCGAACGACGGATGGTCGACGTAGGCATCACCGGTCACGATGATGATGTCGCAGGAATCCCAGCCGAGCAGATCCATCTCCTCCCTGCTCATTGGCAGGAAAGGTGCTGGTCCGAAGCATTCGGCCCAGTACTTGGGATAGTCGTAGAGTGGTTTGGCTGCTTGCATGTCAGTGACCGGTTCTGGTGTGCAGGGAAATCGCGGGCGCGGAATATAGCACAAATTTTGACCAAATCCGACTGGATTCGTGGGGATTGGTGAG  >  NC_002947/5540742‑5541025

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