Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,013,287	G→C	5.6%	L36L (CTG→CTC)	PP_5678 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,013,287	0	G	C	5.6%	296.6 / 6.0	90	L36L (CTG→CTC)	PP_5678	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base C (3/2);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

CCCGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTC  >  NC_002947/5013140‑5013434                                                                                                                                                    |                                                                                                                                                     cccGAGCAGTGCAACAAGCCCCAGACTTAGTGGACAGCCAATGATGATCAATGAGTCGAGCCTCGCAGCCGTTGCGGTGAGCGTAGTTTGCGACCGTTGCTGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCGGTGCTGa                                                                                                                                                     <  6:214295/149‑1 (MQ=255)   cGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAgc                                                                                                                                                   >  1:38733/1‑149 (MQ=255)      gCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTATGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGcgcg                                                                                                                                                <  1:121483/149‑1 (MQ=255)          tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGt                                                                                                                                            <  6:64593/149‑1 (MQ=255)          tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGt                                                                                                                                            >  6:235694/1‑149 (MQ=255)          tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGt                                                                                                                                            >  7:384888/1‑149 (MQ=255)          tGCAACAAGCCTAAGAATTAGGCGAGAGACAATGATGATCAAACAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGAAGGCCAACCACAGCACGCCGTCATGAGCGCGGCGt                                                                                                                                            >  2:24684/1‑149 (MQ=255)             aacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgct                                                                                                                                         <  4:386826/149‑3 (MQ=255)              acaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGaa                                                                                                                                        <  5:235694/149‑1 (MQ=255)                     tCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                  >  6:243913/1‑148 (MQ=255)                     tCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGc                                                                                                                                 >  4:162868/1‑149 (MQ=255)                       aGACTTAGGGGATAGCCAATGATGATAAAATAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCt                                                                                                                                <  1:128738/148‑1 (MQ=255)                               gggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaa                                                                                                                        >  2:2259/1‑148 (MQ=255)                                ggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaa                                                                                                                        >  8:310895/1‑147 (MQ=255)                                ggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaa                                                                                                                        <  7:310895/147‑1 (MQ=255)                                 gagaGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAAc                                                                                                                     <  3:162868/149‑1 (MQ=255)                                  agagCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCTGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGTCCATCCACAGCACTCCGTGCTGAGCGCTGCGTGGAAATCAGGCTCGCAGCAAAAcc                                                                                                                    <  1:24684/149‑1 (MQ=255)                                      ccAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCg                                                                                                                 >  5:337276/1‑148 (MQ=255)                                        aaTGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAg                                                                                                              <  1:25220/149‑1 (MQ=255)                                            atgatCAAAGAGTCGCCCCTCGCAGCCGGTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTGCGTCTGATGAGGGCCCTCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCCGGAAGCGCCCGCAGC‑Ctc                                                                                                         >  8:298167/1‑148 (MQ=255)                                            atgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCATAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑Ct                                                                                                          <  5:148763/148‑1 (MQ=255)                                             tgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTAt                                                                                                        <  3:192918/149‑1 (MQ=255)                                             tgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTAt                                                                                                        <  2:38733/149‑1 (MQ=255)                                                  aaaGAGTCGACCCTCGCAGCCGTTGCTGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACt                                                                                                   <  1:2259/149‑1 (MQ=255)                                                  aaaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAgc                                                                                                                                                   >  1:60951/1‑102 (MQ=255)                                                  aaaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAgc                                                                                                                                                   <  2:60951/102‑1 (MQ=255)                                                   aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTc                                                                                                  >  6:116364/1‑149 (MQ=255)                                                    agagTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTcc                                                                                                 <  5:116364/149‑1 (MQ=255)                                                    agagTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTc                                                                                                  >  1:266206/1‑148 (MQ=255)                                                        tCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCt                                                                                             >  6:196722/1‑149 (MQ=255)                                                             ccTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGc                                                                                         >  5:21979/1‑148 (MQ=255)                                                              cTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                  >  3:314293/1‑107 (MQ=255)                                                              cTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                  <  4:314293/107‑1 (MQ=255)                                                               tCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGTCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGc                                                                                         >  5:58450/1‑146 (MQ=255)                                                                 gCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCgaga                                                                                     <  3:87570/148‑1 (MQ=255)                                                                    gCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaaa                                                                                                                       <  2:100979/112‑1 (MQ=255)                                                                    gCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaaa                                                                                                                       >  1:100979/1‑112 (MQ=255)                                                                    gCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAgctg                                                                                 <  5:43933/149‑1 (MQ=255)                                                                           cGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCttt                                                                          >  4:330424/1‑149 (MQ=255)                                                                            ggTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTg                                                                         <  2:266206/149‑1 (MQ=255)                                                                                       tttGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAAt                                                                                                                                      >  7:232135/1‑78 (MQ=255)                                                                                       tttGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAAt                                                                                                                                      <  8:232135/78‑1 (MQ=255)                                                                                       tttGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCgaga                                                                                     >  7:350081/1‑126 (MQ=255)                                                                                       tttGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCgaga                                                                                     <  8:350081/126‑1 (MQ=255)                                                                                          gCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCaa                                                             <  7:95802/147‑1 (MQ=255)                                                                                          gCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCaa                                                             >  8:95802/1‑147 (MQ=255)                                                                                          gCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAAc                                                            <  6:337276/148‑1 (MQ=255)                                                                                           cGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGa                                                          <  6:21979/149‑1 (MQ=255)                                                                                             aCCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGAt                                                         >  3:197656/1‑148 (MQ=255)                                                                                             aCCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGAt                                                         >  3:204022/1‑148 (MQ=255)                                                                                                gTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAg                                                     >  2:88348/1‑149 (MQ=255)                                                                                                      gCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGaa                                               <  3:58219/149‑1 (MQ=255)                                                                                                        cAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACg                                             >  2:313211/1‑149 (MQ=255)                                                                                                         aGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAACCCCCGAAGC‑CTATGAACTCCATCttt                                                                                           >  7:379614/1‑102 (MQ=255)                                                                                                         aGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGc                                            >  5:279200/1‑149 (MQ=255)                                                                                                          gACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCt                                           <  8:384888/149‑1 (MQ=255)                                                                                                            cTTGTCCGTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTgg                                         <  5:83743/149‑1 (MQ=255)                                                                                                             ttGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTggg                                        >  3:191289/1‑149 (MQ=255)                                                                                                             ttGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTggg                                        >  6:134447/1‑149 (MQ=255)                                                                                                                  cTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCaa                                   <  6:58450/149‑1 (MQ=255)                                                                                                                   ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAt                                  >  1:354512/1‑149 (MQ=255)                                                                                                                    tCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAtt                                 >  1:282479/1‑149 (MQ=255)                                                                                                                        atgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGAt                              <  3:330424/148‑1 (MQ=255)                                                                                                                        atgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGAt                              <  5:243913/148‑1 (MQ=255)                                                                                                                        atgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATg                             >  3:9710/1‑149 (MQ=255)                                                                                                                             ggCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGAt                                                         >  7:152282/1‑116 (MQ=255)                                                                                                                             ggCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGAt                                                         <  8:152282/116‑1 (MQ=255)                                                                                                                             ggCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTtctc                        >  7:146369/1‑149 (MQ=255)                                                                                                                                  tCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAg                    >  2:188644/1‑148 (MQ=255)                                                                                                                                   ccACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAgg                   >  4:354099/1‑148 (MQ=255)                                                                                                                                    atttccACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAgg                   >  1:75152/6‑147 (MQ=255)                                                                                                                                      caGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGc               <  2:354512/149‑1 (MQ=255)                                                                                                                                      caGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGc               >  2:126549/1‑149 (MQ=255)                                                                                                                                       aGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCt              >  7:298865/1‑149 (MQ=255)                                                                                                                                         cACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑ATATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTCGGCCCAATTGATGTTCTCTGAGGAAGCTGa            >  3:198710/1‑149 (MQ=255)                                                                                                                                          aCGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGAc           >  6:103973/1‑149 (MQ=255)                                                                                                                                            gCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATATATCGTGAACGCTGGGACCAATTGATGTTCTCTGAGGAAGCTGAca          >  4:54530/1‑147 (MQ=255)                                                                                                                                            gCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTc         >  6:1995/1‑149 (MQ=255)                                                                                                                                             ccGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTcc        <  3:354099/149‑1 (MQ=255)                                                                                                                                              cGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCt       <  2:282479/149‑1 (MQ=255)                                                                                                                                               gTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCtt      <  4:191289/149‑1 (MQ=255)                                                                                                                                               gTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCtt      <  4:197656/149‑1 (MQ=255)                                                                                                                                               gTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCtt      <  4:204022/149‑1 (MQ=255)                                                                                                                                                tGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTa     <  1:188644/149‑1 (MQ=255)                                                                                                                                                  cTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCTCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTAc    >  1:84447/1‑149 (MQ=255)                                                                                                                                                  cTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGc                                                                        >  3:124870/1‑80 (MQ=255)                                                                                                                                                  cTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGc                                                                        <  4:124870/80‑1 (MQ=255)                                                                                                                                                   tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGAc           >  2:100308/1‑140 (MQ=255)                                                                                                                                                   tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGAc           <  1:100308/140‑1 (MQ=255)                                                                                                                                                   tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTActc  >  6:149910/1‑149 (MQ=255)                                                                                                                                                    gAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTActc  >  4:253706/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCCGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGC‑CTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTC  >  NC_002947/5013140‑5013434

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