Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	668,263	G→T	7.8%	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	7.8%	199.8 / 6.6	64	F87L (TTC→TTA) ‡	PP_0571	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base T (3/2);  total (32/32)
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.18e-01

CGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAAT  >  NC_002947/668116‑668410                                                                                                                                                    |                                                                                                                                                    cgAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATg                                                                                                                                                     >  2:119366/1‑148 (MQ=255)  gAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGaa                                                                                                                                                   >  2:60130/1‑149 (MQ=255)     gcgtaaaGCCTGGGGCTGTTGTCCAGCCCTGTCGCGACTCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATgc                                                                                                                                                <  5:188247/145‑1 (MQ=255)       tGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGACGATGATCGAGGCGCGCATGAATgcgc                                                                                                                                              >  8:186885/1‑149 (MQ=255)         aaaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGcgcg                                                                                                                                             >  2:142550/1‑148 (MQ=255)               tGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                                      <  2:240501/149‑1 (MQ=255)                agggCTGCACAGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCCCAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                                      >  3:100059/2‑148 (MQ=255)                 gggCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa                                                                                                                                    <  3:233332/149‑1 (MQ=255)                   gctgctGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGt                                                                                                                                  <  1:14965/149‑1 (MQ=255)                       ctgTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAca                                                                                                                               <  8:208538/148‑1 (MQ=255)                            cAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAgcagca                                                                                                                         <  1:142550/149‑1 (MQ=255)                                    tCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGc                                                                                                                   <  7:62772/147‑1 (MQ=255)                                      gcgACGCAAGGCCGCTCCCAGAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgat                                                                                                               >  1:31957/1‑149 (MQ=255)                                       cgACGCAAGGCCGCTCCCAGAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatg                                                                                                              <  2:31957/149‑1 (MQ=255)                                         aCGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatgat                                                                                                            >  2:190794/1‑149 (MQ=255)                                          cGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatgat                                                                                                            >  6:118041/1‑148 (MQ=255)                                           gCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcct                                                                                                                                          <  7:235671/117‑4 (MQ=255)                                           gCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcct                                                                                                                                          >  8:235671/1‑114 (MQ=255)                                             aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGGCATATGATACCGAGGATGCTCGAGGCGCGCATGAATCCCCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAgt                                                                                                        >  8:71791/1‑148 (MQ=255)                                             aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                                      >  2:163163/1‑119 (MQ=255)                                             aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                                      <  1:163163/119‑1 (MQ=255)                                             aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAgg                                                                                                        <  1:190794/149‑1 (MQ=255)                                                 ccGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGc                                                                                                      >  2:203609/1‑147 (MQ=255)                                                   gCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAgg                                                                                                  <  5:118041/149‑1 (MQ=255)                                                    cTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGa                                                                                                 >  3:247599/1‑149 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcagc                                                                                                                             <  3:15654/118‑1 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcagc                                                                                                                             >  4:15654/1‑118 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgat                                                                                                               >  7:259964/1‑132 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgat                                                                                                               <  8:259964/132‑1 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGc                                                                                              >  8:254289/1‑149 (MQ=255)                                                         acaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCACGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCt                                                                                             <  1:119366/148‑1 (MQ=255)                                                          caAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTca                                                                                           <  1:203609/149‑1 (MQ=255)                                                           aaaGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTcag                                                                                          >  5:224249/1‑149 (MQ=255)                                                                  gcgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCg                                                                                     >  1:124048/1‑147 (MQ=255)                                                                  gcgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                                   <  1:60130/149‑1 (MQ=255)                                                                         tCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCa                                                                            <  6:224249/149‑1 (MQ=255)                                                                            gTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTc                                                                              <  4:304639/144‑1 (MQ=255)                                                                            gTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTc                                                                              <  4:244091/144‑1 (MQ=255)                                                                            gTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTc                                                                              >  3:244091/1‑144 (MQ=255)                                                                            gTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTc                                                                              >  3:304639/1‑144 (MQ=255)                                                                               gATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGccc                                                                       >  4:92425/1‑148 (MQ=255)                                                                                 tCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAg                                                                    <  3:92425/149‑1 (MQ=255)                                                                                      cAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTgcg                                                               >  1:208252/1‑149 (MQ=255)                                                                                          tcatcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGa                                                            >  3:136215/1‑148 (MQ=255)                                                                                            ctcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCTGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAg                                                         <  8:229724/148‑1 (MQ=255)                                                                                               aCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGcc                                                      <  4:19781/149‑1 (MQ=255)                                                                                               aCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGc                                                       >  6:169074/1‑148 (MQ=255)                                                                                                   ataGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATc                                                   >  8:221579/3‑148 (MQ=255)                                                                                                     aGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGAt                                                >  7:85311/1‑149 (MQ=255)                                                                                                           aGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGtt                                          <  4:100059/149‑1 (MQ=255)                                                                                                             cACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGc                                        >  5:31275/1‑149 (MQ=255)                                                                                                                cagccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAg                                            <  5:66150/142‑1 (MQ=255)                                                                                                                cagccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAg                                            >  6:66150/1‑142 (MQ=255)                                                                                                                         gatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTgg                            >  7:11268/1‑149 (MQ=255)                                                                                                                         gatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTgg                            >  1:25243/1‑149 (MQ=255)                                                                                                                          atTCCGAGGATGATCGAGGCGCGCATGATTGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATTATCAGGTCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTgg                            <  7:71791/148‑1 (MQ=255)                                                                                                                           ttCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGcc                          <  2:25243/149‑1 (MQ=255)                                                                                                                              cGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCcagc                       <  7:254289/149‑1 (MQ=255)                                                                                                                                       tCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaa               <  6:31275/148‑1 (MQ=255)                                                                                                                                           tccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCAtt                              <  3:40075/127‑1 (MQ=255)                                                                                                                                           tccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCAtt                              >  4:40075/3‑129 (MQ=255)                                                                                                                                            gcgcgcATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAcccc         >  7:255110/1‑149 (MQ=255)                                                                                                                                            ccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAAccc         >  1:61250/2‑149 (MQ=255)                                                                                                                                                   tGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAAt  <  4:247599/149‑1 (MQ=255)                                                                                                                                                   tGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAAt  <  5:169074/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAAT  >  NC_002947/668116‑668410

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