Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,488,043	C→T	7.9%	R139W (CGG→TGG)	PP_2182 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,488,043	0	C	T	7.9%	179.7 / 4.5	63	R139W (CGG→TGG)	PP_2182	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base T (3/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

ACAGCTCAAGCTGGGCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATT  >  NC_002947/2487895‑2488190                                                                                                                                                     |                                                                                                                                                    aCAGCTCAAGCTGGGCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTGAGGCGCTGCATGCCAAGCAGCAGGCTgcgc                                                                                                                                                     <  1:35346/149‑1 (MQ=255)            tGGGCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCGCGCAGCCTGCTGCTTGGCatgca                                                                                                                                             <  7:363192/146‑5 (MQ=255)            tGGGCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCGCGCAGCCTGCTGCTTGGCatgca                                                                                                                                             >  8:363192/1‑142 (MQ=255)            tGGGCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGc                                                                                                                                          >  2:219117/1‑149 (MQ=255)              ggCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGc                                                                                                                                          >  1:338975/1‑147 (MQ=255)                cATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAg                                                                                                                                       >  5:30075/1‑148 (MQ=255)                  tttCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGAc                                                                                                                                    >  7:306873/1‑149 (MQ=255)                       aCGTGGTGCCGGGGCGCTCGATACTGCGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCCCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACg                                                                                                                               >  7:145616/1‑149 (MQ=255)                          tggtgCCGGGGCGCTCGATAATGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCaa                                                                                                                            >  5:247643/1‑149 (MQ=255)                                cGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGc                                                                                                                      <  1:374177/149‑1 (MQ=255)                                 ggggCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGc                                                                                                                      >  8:291688/1‑148 (MQ=255)                                 ggggCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCg                                                                                                                     <  5:303549/149‑1 (MQ=255)                                          gATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCaaa                                                                                                             <  4:1919/148‑1 (MQ=255)                                          gATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCaa                                                                                                              <  2:338975/147‑1 (MQ=255)                                          gATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAag                                                                                                            >  5:261080/1‑149 (MQ=255)                                               tGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAaggcc                                                                                                         >  6:234378/1‑147 (MQ=255)                                                 ggCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAaggcccac                                                                                                      <  7:317530/148‑2 (MQ=255)                                                    gAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAaggcccaggccc                                                                                                  <  3:372987/149‑1 (MQ=255)                                                     aGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAggcccaggccca                                                                                                 <  4:161449/149‑1 (MQ=255)                                                       cATGGTGGAGAATCCACCGCGCCTTGACCTCAATGGGCAGGCCGCAGGTACTGTGACTGCGTCTCAGGCGCTGCATGCCAAGCACCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGcccaggcccagg                                                                                               <  1:319246/149‑1 (MQ=255)                                                            tggAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCTGGCTGCGTTGCCGCGCAGcctgct                                                                                                                                     >  3:238376/1‑102 (MQ=255)                                                            tggAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCTGGCTGCGTTGCCGCGCAGcctgct                                                                                                                                     <  4:238376/106‑5 (MQ=255)                                                            tggAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCg                                                                                          >  8:169507/1‑149 (MQ=255)                                                                       ccGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCccgccgc                                                                               <  6:289911/149‑1 (MQ=255)                                                                          cgcCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGa                                                                            <  7:291688/149‑1 (MQ=255)                                                                              ttGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGTGAAGTCAAAGGCCCAGGCCCAGTCTCCGGCCCccgccgccg                                                                             <  2:65402/144‑1 (MQ=255)                                                                              ttGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGTGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCAccgccgccg                                                                             >  1:65402/1‑144 (MQ=255)                                                                                 aCCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGCCACCGc                                                                      >  6:167683/1‑148 (MQ=255)                                                                                   cTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAAc                                                                   <  5:167683/149‑1 (MQ=255)                                                                                           gCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGt                                                           >  1:297803/1‑149 (MQ=255)                                                                                                  gCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCACCGGTGCAGTCTGCCGa                                                    >  8:243772/1‑149 (MQ=255)                                                                                                    aGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACt                                                  <  8:306873/149‑1 (MQ=255)                                                                                                     ggTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTg                                                 <  6:247643/149‑1 (MQ=255)                                                                                                      gTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACt                                                  >  6:177823/1‑147 (MQ=255)                                                                                                      gTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGa                                                <  5:234378/149‑1 (MQ=255)                                                                                                          tgtgACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGcg                                             >  7:326933/1‑147 (MQ=255)                                                                                                           gtgACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTa                                                        <  7:309077/136‑1 (MQ=255)                                                                                                           gtgACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTa                                                        >  8:309077/1‑136 (MQ=255)                                                                                                           gtgACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCt                                             >  7:63396/1‑147 (MQ=255)                                                                                                            tgACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAg                                          >  8:44201/1‑149 (MQ=255)                                                                                                             gACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGACTACCGACTGAGAACAGt                                         >  4:38292/1‑149 (MQ=255)                                                                                                              aCCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTc                                        <  6:30075/149‑1 (MQ=255)                                                                                                               ccgccgCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTcg                                       >  1:190065/1‑149 (MQ=255)                                                                                                                  ccgcTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGt                                    <  1:219117/149‑1 (MQ=255)                                                                                                                   cgcTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGtt                                   >  5:374877/1‑149 (MQ=255)                                                                                                                     cTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCt                                 >  3:250654/1‑149 (MQ=255)                                                                                                                     cTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCt                                 >  5:345937/1‑149 (MQ=255)                                                                                                                        aGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGcc                              >  5:181315/1‑149 (MQ=255)                                                                                                                            gcTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCg                           >  8:347835/1‑148 (MQ=255)                                                                                                                            gcTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCg                           >  8:21335/1‑148 (MQ=255)                                                                                                                             cTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTa                         <  1:82026/149‑1 (MQ=255)                                                                                                                             cTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTa                         <  7:44201/149‑1 (MQ=255)                                                                                                                                  tGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTAcc                                                      >  6:200980/1‑115 (MQ=255)                                                                                                                                  tGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTAcc                                                      <  5:200980/115‑1 (MQ=255)                                                                                                                                  tGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCAt                    <  6:345937/149‑1 (MQ=255)                                                                                                                                    ccAAGCAGCAGGCTGCGCGGCAACGCAGCCTGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGc                  >  1:34384/1‑149 (MQ=255)                                                                                                                                     cAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCt                 >  7:23642/1‑149 (MQ=255)                                                                                                                                        gcagcaGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGa              <  2:297803/149‑1 (MQ=255)                                                                                                                                        gcagcaGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGa              <  6:181315/149‑1 (MQ=255)                                                                                                                                         cagcagGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAt             <  4:250654/149‑1 (MQ=255)                                                                                                                                             tggCTGCGTTGCCGCGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCa         >  7:270848/2‑149 (MQ=255)                                                                                                                                             agGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCa         <  8:63396/149‑1 (MQ=255)                                                                                                                                             agGCGGCGCGGCAAGGCAGCCAGGACCGACGCAATCGTTCGAAGTCAAAGGCCAAGGCCCAGTCTCCGTCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCTACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCa         <  7:340900/149‑1 (MQ=255)                                                                                                                                               gCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGc       >  6:158821/1‑149 (MQ=255)                                                                                                                                                    gcgGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCAtt  <  5:177823/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ACAGCTCAAGCTGGGCATTTCCACGTGGTGCCGGGGCGCTCGATACTGGGCGAGCATGGTGGAGAATGCACCGCGCCTTGACCTGAATGGGCAGGCCGCAGGTACTGTGACCGCCGCTCAGGCGCTGCATGCCAAGCAGCAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATT  >  NC_002947/2487895‑2488190

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