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%	200.5 / 6.7	68	D49H (GAT→CAT)	PP_2294	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (2/4);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGG  >  NC_002947/2621184‑2621472                                                                                                                                                    |                                                                                                                                               aGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGa                                                                                                                                               <  6:129622/148‑1 (MQ=255) aGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGAt                                                                                                                                              >  5:345869/1‑149 (MQ=255)        gTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAgg                                                                                                                                       >  5:313891/1‑149 (MQ=255)         tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgca                                                                                                                                      <  2:27429/149‑3 (MQ=255)         tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgc                                                                                                                                       <  5:404851/148‑2 (MQ=255)         tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTACTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAATTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGt                                                                                                                                      >  3:273626/1‑149 (MQ=255)           gagaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGtt                                                                                                                                     <  3:260403/148‑1 (MQ=255)           gagaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcat                                                                                                                                     <  2:281655/148‑4 (MQ=255)            agaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcatag                                                                                                                                   <  2:122107/149‑6 (MQ=255)                         ggATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCaaa                                                                                                                      <  5:326358/149‑1 (MQ=255)                          gATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCaaa                                                                                                                      >  6:225982/1‑148 (MQ=255)                                 ttGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGa                                                                                                              >  2:165876/1‑149 (MQ=255)                                   ggCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAg                                                                                                                                        <  8:10954/121‑1 (MQ=255)                                   ggCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAg                                                                                                                                        >  7:10954/1‑121 (MQ=255)                                     ccATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCag                                                                                                          <  7:394304/149‑1 (MQ=255)                                     ccATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCag                                                                                                          <  6:345869/149‑1 (MQ=255)                                        ttGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCag                                                                                                          >  7:255756/1‑146 (MQ=255)                                           gAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAgaagaa                                                                                                     <  5:225982/148‑1 (MQ=255)                                               tgctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCATTAAAGTGATCCAGAGGTTTATCCCGAAAGCGCAGTGCCCTATGCTAAGATTGATGCTGAGGTTGGCTGCAAGAGAAAATACAGTGATCAGAAGAAGGa                                                                                                  >  1:103037/1‑148 (MQ=255)                                                  ttATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACa                                                                                             >  6:151467/1‑149 (MQ=255)                                                    aTGTTTTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAATGCTCCAGAGGTTTATCACGAAAGCGCACATACCTATGCTCAGCTTCATGCTGAGGTTTGCTGCAAGATCAAATACAGTCATCAGAAGAAGTATGACAtt                                                                                           >  1:101296/1‑149 (MQ=255)                                                          tCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGa                                                                                      <  6:313891/148‑1 (MQ=255)                                                           ccTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGAtt                                                                                    >  7:95079/1‑149 (MQ=255)                                                            cTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAag                                                                                                                            >  6:59522/1‑108 (MQ=255)                                                            cTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAag                                                                                                                            <  5:59522/108‑1 (MQ=255)                                                            cTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGCGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTc                                                                                   >  1:292089/1‑149 (MQ=255)                                                             tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                  >  7:194131/1‑149 (MQ=255)                                                             tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                  >  2:430206/1‑149 (MQ=255)                                                             tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGAAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                  >  5:444635/1‑149 (MQ=255)                                                               aaaGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCct                                                                                 <  8:136683/148‑1 (MQ=255)                                                                        gCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACaaaa                                                                        <  3:343969/148‑1 (MQ=255)                                                                          tAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAAt                                                                      >  3:41244/1‑148 (MQ=255)                                                                            aCGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATc                                                                     <  1:430206/147‑1 (MQ=255)                                                                              gcgcCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTg                                                                 >  3:171568/1‑149 (MQ=255)                                                                               cgcCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTg                                                                 >  1:192351/1‑148 (MQ=255)                                                                               cgcCTAGCTCTTCCTC‑TTCAAGCGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTgg                                                                <  2:292089/149‑1 (MQ=255)                                                                                 ccTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGAt                                                              <  5:68265/149‑1 (MQ=255)                                                                                  cTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTg                                                                 >  6:419268/1‑145 (MQ=255)                                                                                  cTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTg                                                                 <  5:419268/145‑1 (MQ=255)                                                                                  cTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATg                                                             <  2:192351/149‑1 (MQ=255)                                                                                     gCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGa                                                            <  2:50603/147‑1 (MQ=255)                                                                                       tcttcctc‑ttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                        >  4:13806/1‑149 (MQ=255)                                                                                                ttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAgg                                                  <  5:19382/147‑1 (MQ=255)                                                                                                  cAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTc                                               <  5:151467/148‑1 (MQ=255)                                                                                                     gTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAAtt                                           >  6:290033/1‑149 (MQ=255)                                                                                                        cTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAc                                        >  8:223323/1‑149 (MQ=255)                                                                                                         gtcAGAGGTTTATCCCGAAAGCGCATTGACCTTTGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAATGCTCAATTGAc                                        <  8:143543/146‑1 (MQ=255)                                                                                                          ccAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCt                                      >  7:29376/1‑149 (MQ=255)                                                                                                           cAGAGGTTTATCCAGAAATCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTATCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                     <  4:273626/149‑1 (MQ=255)                                                                                                            agagGTTTATCCCGAAAGCGCAGTGACCTATTCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCaaattac                                        >  7:143543/1‑140 (MQ=255)                                                                                                             gagGTTTTTCCCGAAAGCGCAGTGACCAATGCTGAGCTTGATGCTGATGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGTTGACATTTTCGATTCCTCTTTCAAAAATCACTGTATGACGTTGAAAGGCTCAATTTACCtttt                                   <  2:101296/149‑3 (MQ=255)                                                                                                                 tttATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGa                                <  7:223323/148‑1 (MQ=255)                                                                                                                   tATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAtt                              <  8:95079/148‑1 (MQ=255)                                                                                                                     tCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCa                            <  4:171568/148‑1 (MQ=255)                                                                                                                     tCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAg                           <  5:290033/149‑1 (MQ=255)                                                                                                                       ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGa                          >  6:239704/1‑148 (MQ=255)                                                                                                                       ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGa                          >  5:254787/1‑148 (MQ=255)                                                                                                                             gcgcAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTa                   <  4:41244/149‑1 (MQ=255)                                                                                                                              cgcAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTAt                  <  1:165876/149‑1 (MQ=255)                                                                                                                                 aGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGTCTCAATTGACCTTGTGGATTCAGACAACGTATCGc               <  6:444635/149‑1 (MQ=255)                                                                                                                                 aGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAt                               >  4:255905/1‑133 (MQ=255)                                                                                                                                 aGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAt                               <  3:255905/133‑1 (MQ=255)                                                                                                                                      ccTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCACTCaa           <  8:255756/148‑1 (MQ=255)                                                                                                                                       cTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc         >  5:176187/1‑149 (MQ=255)                                                                                                                                          accTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAAcc         >  3:407921/3‑146 (MQ=255)                                                                                                                                            cTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGa     <  8:194131/148‑1 (MQ=255)                                                                                                                                             tCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACg   >  4:193722/1‑149 (MQ=255)                                                                                                                                              cAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACgg  >  2:313613/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               AGATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTC‑TTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGG  >  NC_002947/2621184‑2621472

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑