Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,621,332	T→A	6.2%	D49E (GAT→GAA) ‡	PP_2294 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,621,332	0	T	A	6.2%	195.8 / 2.9	64	D49E (GAT→GAA) ‡	PP_2294	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (33/31)
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.03e-01

GATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAGGCA  >  NC_002947/2621185‑2621480                                                                                                                                                    |                                                                                                                                                      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                                                                                                         <  2:327129/134‑1 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                         >  1:327129/1‑134 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                          >  6:241968/1‑149 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                          >  1:201243/1‑149 (MQ=255)                                                            tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                          >  2:288540/1‑149 (MQ=255)                                                             cAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAAt                                                                                                                             >  5:93750/1‑113 (MQ=255)                                                             cAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAAt                                                                                                                             <  6:93750/113‑1 (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                                                                        <  3:205515/135‑1 (MQ=255)                                                                                            tcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTgg                                                                        >  4:205515/1‑135 (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)                                                                                                                                                  gATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAg     <  2:201243/148‑1 (MQ=255)                                                                                                                                                    tGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAGGCa  <  5:241968/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGA‑GGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAGGCA  >  NC_002947/2621185‑2621480

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