Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,621,328	T→A	6.3%	L48H (CTT→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,328	0	T	A	6.3%	193.4 / 2.9	63	L48H (CTT→CAT)	PP_2294	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

GATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGA  >  NC_002947/2621185‑2621474                                                                                                                                                |                                                                                                                                                    gatgatGTCGAGACGGGCGGTTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATg                                                                                                                                                >  3:185300/1‑149 (MQ=255) gatgatGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATg                                                                                                                                                >  3:188788/1‑149 (MQ=255)      tgTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑g                                                                                                                                          <  2:35762/149‑1 (MQ=255)        tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGt                                                                                                                                        >  5:28229/1‑149 (MQ=255)          gagaCGGGCGATTGGTATCGTGGTGGCCATTTGAGCTGCTTATGTTGTACTCAAAGGAGTAGCTAACGCGCCTAGCTCTTGCTCTTGAAGTGCTCGAGTGGTTTTTCCCGTAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑TGtt                                                                                                                                       <  7:87451/148‑1 (MQ=255)          gagaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGA‑gcat                                                                                                                                       <  5:106888/148‑4 (MQ=255)          gagaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGA‑gcat                                                                                                                                       <  6:160762/148‑4 (MQ=255)                 gCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCaa                                                                                                                               <  5:128835/149‑1 (MQ=255)                 gCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCaa                                                                                                                               <  8:331759/149‑1 (MQ=255)                          aTCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATa                                                                                                                      <  6:374789/149‑1 (MQ=255)                              tgtTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGt                                                                                                                  <  8:103076/149‑1 (MQ=255)                                   gCCATTTGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCa                                                                                                             <  8:77945/149‑1 (MQ=255)                                     cATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCag                                                                                                            >  3:91978/1‑148 (MQ=255)                                          gAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGaagaag                                                                                                      >  6:152550/1‑149 (MQ=255)                                             ctgctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                   <  4:269098/149‑1 (MQ=255)                                                   aTGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACAtt                                                                                             >  5:38764/1‑149 (MQ=255)                                                   aTGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACAt                                                                                              <  3:14296/148‑1 (MQ=255)                                                         tCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGAt                                                                                       >  5:239830/1‑149 (MQ=255)                                                           cTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGAtt                                                                                      >  3:194533/1‑148 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                   >  1:327129/1‑134 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                   <  2:327129/134‑1 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                    >  2:288540/1‑149 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                    >  1:201243/1‑149 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                    >  6:241968/1‑149 (MQ=255)                                                             cAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAAt                                                                                                                       <  6:93750/113‑1 (MQ=255)                                                             cAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAAt                                                                                                                       >  5:93750/1‑113 (MQ=255)                                                                 ggAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTa                                                                               <  4:185300/149‑1 (MQ=255)                                                                  gAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTTGCTTCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTAc                                                                              <  8:173146/149‑1 (MQ=255)                                                                       gCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACaaaaa                                                                         >  6:282684/1‑149 (MQ=255)                                                                          aaCGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCa                                                                      >  5:36691/1‑149 (MQ=255)                                                                            cgcgCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACt                                                                    <  6:92941/149‑1 (MQ=255)                                                                               gcCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGCATGACATTTTCGATTCCTCTTACAAAAATCACTGGa                                                                 >  6:364701/1‑149 (MQ=255)                                                                                      tcttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                          >  5:37502/1‑149 (MQ=255)                                                                                        ttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GATTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGAt                                                                                       >  1:15226/1‑118 (MQ=255)                                                                                        ttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GATTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGAt                                                                                       <  2:15226/118‑1 (MQ=255)                                                                                            tcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTgg                                                                  >  4:205515/1‑135 (MQ=255)                                                                                            tcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTgg                                                                  <  3:205515/135‑1 (MQ=255)                                                                                            tcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGaaac                                                     >  1:136750/1‑147 (MQ=255)                                                                                              ttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAgg                                                    <  3:297447/147‑1 (MQ=255)                                                                                                cAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGTATGACATTTTCGATTCGTCTTACAAAAATCACTGGATGACCATGAAAGGCTc                                                 >  5:371070/1‑148 (MQ=255)                                                                                                 aaGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAATTACAGTGAt                                                                                                               <  2:316072/85‑1 (MQ=255)                                                                                                   gTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAAtt                                             <  6:36691/149‑1 (MQ=255)                                                                                                    tGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTg                                            <  6:239830/149‑1 (MQ=255)                                                                                                      cTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGa                                           <  4:194533/148‑1 (MQ=255)                                                                                                         cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAcc                                         <  4:297392/147‑1 (MQ=255)                                                                                                         cAGAGGTTTATCCCGAAAGCGCAGTGACCAATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                       >  1:219900/1‑149 (MQ=255)                                                                                                          agagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                       >  1:283836/1‑149 (MQ=255)                                                                                                           gagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                     <  6:142861/149‑1 (MQ=255)                                                                                                                  aTCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAtt                                >  4:34365/1‑147 (MQ=255)                                                                                                                  aTCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAtt                                >  3:24880/1‑147 (MQ=255)                                                                                                                   tCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGaaa                                                      >  3:304650/1‑124 (MQ=255)                                                                                                                   tCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGaaa                                                      <  4:304650/124‑1 (MQ=255)                                                                                                                     ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATTACATTTTCGATTCCTCTTACAAAAATCACTGTAAGACCATGAAAGGCTCAATTGACATTGTGGATTCata                            >  4:26114/1‑146 (MQ=255)                                                                                                                     ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGa                            >  6:194048/1‑148 (MQ=255)                                                                                                                        aaaGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAAc                        <  6:28229/149‑1 (MQ=255)                                                                                                                                     ccaaccTCAGCATCAAGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAc            >  2:378225/6‑148 (MQ=255)                                                                                                                                      tATGCTCATCTTGATGCTGA‑GGTTTGCTGCAATAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCt          <  3:26114/149‑1 (MQ=255)                                                                                                                                      tATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGCATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc           >  2:135219/1‑148 (MQ=255)                                                                                                                                       aaccTCAGCATCAAGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc           >  7:26889/4‑147 (MQ=255)                                                                                                                                       aTGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc           >  6:231947/1‑147 (MQ=255)                                                                                                                                             aGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACgg    <  6:37502/148‑1 (MQ=255)                                                                                                                                             aGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACg     <  5:231947/147‑1 (MQ=255)                                                                                                                                               cTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGa  >  3:207786/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    GATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGA  >  NC_002947/2621185‑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: ‑