Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,621,330	G→C	8.5%	D49H (GAT→CAT)	PP_2294 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,621,330	0	G	C	8.5%	211.4 / 7.0	69	D49H (GAT→CAT)	PP_2294	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base C (4/2);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGA  >  NC_002947/2621184‑2621474                                                                                                                                                   |                                                                                                                                                  aGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGa                                                                                                                                                  <  8:146233/148‑1 (MQ=255)  gatgatGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATg                                                                                                                                                <  7:297843/149‑1 (MQ=255)   atgatgTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGc                                                                                                                                               >  5:124749/1‑149 (MQ=255)      atgTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTg                                                                                                                                             <  4:380894/148‑1 (MQ=255)            agaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcatag                                                                                                                                      <  1:172750/149‑6 (MQ=255)                   cGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGc                                                                                                                                  >  2:119245/1‑146 (MQ=255)                    gATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAaga                                                                                                                              <  3:21621/149‑1 (MQ=255)                    gATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAaga                                                                                                                              <  1:328229/149‑1 (MQ=255)                         ggATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCaaa                                                                                                                         >  6:310660/1‑149 (MQ=255)                         ggATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCaaa                                                                                                                         <  6:35095/149‑1 (MQ=255)                            tCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATa                                                                                                                       >  1:144077/1‑148 (MQ=255)                              gtgtTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACa                                                                                                                     <  7:234308/148‑1 (MQ=255)                                 ttGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTg                                                                                                                  <  6:124749/148‑1 (MQ=255)                                    gCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGAt                                                                                                                >  1:255147/1‑147 (MQ=255)                                       aTTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCag‑aa                                                                                                          <  3:175742/149‑1 (MQ=255)                                       aTTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCag                                                                                                             >  3:286977/1‑147 (MQ=255)                                             gctgctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGa                                                                                                    >  2:29111/1‑149 (MQ=255)                                             gctgctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGa                                                                                                    >  2:67452/1‑149 (MQ=255)                                               tgctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGAt                                                                                                   >  7:345645/1‑148 (MQ=255)                                                gctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGa                                                                                                 >  6:160850/1‑149 (MQ=255)                                                 ctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGAt                                                                                                   <  5:284196/146‑1 (MQ=255)                                                 ctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGa                                                                                                 >  2:388047/1‑148 (MQ=255)                                                       ttgtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTc                                                                                          >  6:4613/1‑149 (MQ=255)                                                       ttgtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTc                                                                                          <  8:174769/149‑1 (MQ=255)                                                          tCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGa                                                                                        <  5:4613/148‑1 (MQ=255)                                                            cTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTc                                                                                     <  8:345645/149‑1 (MQ=255)                                                             tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTcc                                                                                    >  6:247283/1‑149 (MQ=255)                                                             tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTcc                                                                                    >  4:84376/1‑149 (MQ=255)                                                             tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTcc                                                                                    >  6:304257/1‑149 (MQ=255)                                                              cAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCct                                                                                   >  2:251318/1‑149 (MQ=255)                                                               aaaGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCct                                                                                   <  5:336929/148‑1 (MQ=255)                                                                      tAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACaa                                                                            >  1:207039/1‑148 (MQ=255)                                                                        gCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTTACCTATGCTCAGCTTGATGCTGAGGTTGGCTTCAAGAGCAAATACAGTGATCAGAAAGAAGGATGACATTTTCGATTCCTCTTACaaa                                                                           <  1:251318/148‑1 (MQ=255)                                                                            aCGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCAc                                                                     >  4:180534/1‑149 (MQ=255)                                                                                  cTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATg                                                               <  4:286977/149‑1 (MQ=255)                                                                                     gCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGAc                                                             >  7:9409/1‑148 (MQ=255)                                                                                        cttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                          >  7:358065/1‑148 (MQ=255)                                                                                        cttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                          >  1:151685/1‑148 (MQ=255)                                                                                            ctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGaaa                                                      <  2:326717/148‑1 (MQ=255)                                                                                            ctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGaaa                                                      <  8:358065/148‑1 (MQ=255)                                                                                              cttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGc                                                   >  8:272064/1‑149 (MQ=255)                                                                                                 cAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCt                                                  <  1:119245/147‑1 (MQ=255)                                                                                                    gTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTg                                                                   <  5:398725/127‑1 (MQ=255)                                                                                                    gTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTg                                                                   >  6:398725/1‑127 (MQ=255)                                                                                                       cTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGa                                           <  3:84376/148‑1 (MQ=255)                                                                                                       cTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAc                                          <  5:304257/149‑1 (MQ=255)                                                                                                        tCCAGAGGTTTATCCCGAAAGCGCAGTGCCCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑ATGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAcc                                         >  7:232163/1‑149 (MQ=255)                                                                                                          catagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                       <  5:310660/146‑1 (MQ=255)                                                                                                          cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                       >  8:210170/1‑149 (MQ=255)                                                                                                            gagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                     <  1:388047/149‑1 (MQ=255)                                                                                                            cagGTTTATCCCGAAACCGCAGTGACCAACGCTCCGCTTCATGCCGAGCTAGGCTGCAAGACCAAAAACACTGAGCAC‑AACAAGGATCCCATTTTCGATTCCTCTTACAAAAAACACTGGACCACCACCAAAGCCTCAATAGACCTtgc                                     <  2:144077/148‑2 (MQ=255)                                                                                                             agGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgtg                                    >  1:383424/1‑149 (MQ=255)                                                                                                               gTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGa                                  >  1:38760/1‑149 (MQ=255)                                                                                                                tttATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAt                                 <  2:255147/149‑1 (MQ=255)                                                                                                                      ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGTCTCAATTGACCTTGTGGATTCTGa                            >  4:141580/1‑148 (MQ=255)                                                                                                                       cGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGAc                           >  6:170885/1‑148 (MQ=255)                                                                                                                        gAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACaa                         <  7:272064/149‑1 (MQ=255)                                                                                                                            gcgcAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTa                     >  4:217886/1‑149 (MQ=255)                                                                                                                                 gTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCt                >  1:87058/1‑149 (MQ=255)                                                                                                                                    agccaaccTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTTACCTTGTGGATTCAGACAACGTATCGCTCaa             >  4:46477/8‑149 (MQ=255)                                                                                                                                    aCCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCa              <  3:217886/148‑1 (MQ=255)                                                                                                                                        aaccTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc           >  1:245293/4‑147 (MQ=255)                                                                                                                                        aaccTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTACTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc           >  6:290426/4‑147 (MQ=255)                                                                                                                                         tGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCt          <  5:160850/147‑1 (MQ=255)                                                                                                                                          gCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGa       >  8:360968/1‑149 (MQ=255)                                                                                                                                           cTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGa       <  2:87058/148‑1 (MQ=255)                                                                                                                                           cTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGAc      >  8:262212/1‑149 (MQ=255)                                                                                                                                           cTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTTCAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGCCAACGTATCGCTCAACCTAGAc      >  4:273283/1‑149 (MQ=255)                                                                                                                                                 ttGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGa  >  1:370599/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAG‑AAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGA  >  NC_002947/2621184‑2621474

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