Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,808,839	A→G	7.0%	intergenic (+44/‑111)	kdsA‑I → / → eno	2‑dehydro‑3‑deoxyphosphooctonate aldolase/enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,808,839	0	A	G	7.0%	249.6 / 6.7	86	intergenic (+44/‑111)	kdsA‑I/eno	2‑dehydro‑3‑deoxyphosphooctonate aldolase/enolase
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base G (2/4);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CCATGCGCCCTGCGCCTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGTGAA  >  NC_002947/1808697‑1808982                                                                                                                                                    |                                                                                                                                                ccATGCGCCCTGCGCCTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGt                                                                                                                                           >  4:3114/1‑149 (MQ=255)       gcCCTGCGCCTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCGACTGGTCAgggg                                                                                                                                     <  8:73940/149‑1 (MQ=255)       gcCCTGCGCCTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCGACTGGTCAgggg                                                                                                                                     <  4:155699/149‑1 (MQ=255)        cccTGCGCCTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGc                                                                                                                                    <  6:42523/149‑1 (MQ=255)               ccTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAgg                                                                                                                                       <  5:275012/139‑1 (MQ=255)               ccTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAgg                                                                                                                                       >  6:275012/1‑139 (MQ=255)                cTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTAGCTCAGCGGATCCCTtc                                                                                                                             >  8:223755/1‑148 (MQ=255)                      aaaCTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCG‑CCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAg                                                                                                                                        >  5:64618/1‑130 (MQ=255)                      aaaCTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCG‑CCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAg                                                                                                                                        <  6:64618/130‑1 (MQ=255)                       aaCTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGc                                                                                                                      <  2:407672/148‑1 (MQ=255)                        aCTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCt                                                                                                                    >  1:45948/1‑149 (MQ=255)                         cTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTg                                                                                                                   >  4:243863/1‑149 (MQ=255)                          tGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGc                                                                                                                  >  1:45825/1‑149 (MQ=255)                              gCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCa                                                                                                               >  4:261261/1‑148 (MQ=255)                               ccGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCa                                                                                                               >  2:72143/1‑147 (MQ=255)                                 gTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCGACTGGTCAggggt                                                                                                                                    <  4:11580/124‑2 (MQ=255)                                 gTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCGACTGGTCAggggt                                                                                                                                    >  3:11580/1‑123 (MQ=255)                                     cTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGt                                                                                                       <  6:302344/149‑1 (MQ=255)                                         tCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGc                                                                                                   >  7:144342/1‑149 (MQ=255)                                          ccAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGcc                                                                                                  >  2:372738/1‑149 (MQ=255)                                           cAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGcc                                                                                                  <  2:45825/148‑1 (MQ=255)                                               tCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGt                                                                                              >  7:184632/1‑148 (MQ=255)                                                cAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGt                                                                                              >  2:410820/1‑147 (MQ=255)                                                cAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGt                                                                                              <  1:410820/147‑1 (MQ=255)                                                 aaGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGaa                                                                                           >  8:268657/1‑149 (MQ=255)                                                     aaTTTGACGACCTGTTGAAAAGTTTTCCGACGTTAGAAACCG‑CGTAAATCTCATTTCTCGGGTAAAGTACCGACAGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGc                                                                                       <  7:223755/149‑1 (MQ=255)                                                      aTTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGc                                                                                       >  4:123589/1‑148 (MQ=255)                                                          gacgacCTGGTGAAAAGTTTTCCGACGGTAGAAACCGCCGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAg                                                                                    >  1:72612/1‑148 (MQ=255)                                                          gacgacCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAAc                                                                                          >  4:129523/1‑141 (MQ=255)                                                          gacgacCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAAc                                                                                          <  3:129523/141‑1 (MQ=255)                                                                cTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCct                                                                            <  8:136077/149‑1 (MQ=255)                                                                   gTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCctaag                                                                         >  4:372815/1‑149 (MQ=255)                                                                   gTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCcta                                                                           <  1:282337/147‑1 (MQ=255)                                                                     gAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTc                                                                                                     <  3:26736/119‑1 (MQ=255)                                                                     gAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTc                                                                                                     >  4:26736/1‑119 (MQ=255)                                                                      aaaaGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCctaagct                                                                       <  3:372815/148‑1 (MQ=255)                                                                         aGTTTTCCGACGGTAGAAACCG‑CGTAACGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAgt                                                                   <  3:123589/149‑1 (MQ=255)                                                                          gTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAgtt                                                                  >  2:377596/1‑149 (MQ=255)                                                                              tCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGttt                                                              >  3:274753/1‑149 (MQ=255)                                                                                 gACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCg                                                           >  4:233384/1‑149 (MQ=255)                                                                                  aCGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGt                                                          <  8:243346/149‑1 (MQ=255)                                                                                   cGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTc                                                         <  8:160812/149‑1 (MQ=255)                                                                                      tAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCaaa                                                      >  6:113035/1‑149 (MQ=255)                                                                                      tAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCAACAGGCCTGCCCAAAGCAAGTTCGCCCgg                                                                                               <  1:231362/108‑1 (MQ=255)                                                                                      tAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCAACAGGCCTGCCCAAAGCAAGTTCGCCCgg                                                                                               >  2:231362/1‑108 (MQ=255)                                                                                         aaaCCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTc                                                   >  3:110490/1‑149 (MQ=255)                                                                                          aaCCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCCCTGCCCATCTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCaaa                                                      <  1:136455/149‑1 (MQ=255)                                                                                          aaCCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCt                                                  <  8:184632/149‑1 (MQ=255)                                                                                             cg‑cgTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGa                                               <  5:76669/149‑1 (MQ=255)                                                                                              g‑cgTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCACGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAg                                              <  3:3114/149‑1 (MQ=255)                                                                                                   aaaGCTCATTTCTCGGGTAAAGTAAAGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCAACCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAAAGGCAAGAATTTACTAAGCTAATTTGTTTTCGTCAAATTCTGGAGTCCt                                          >  4:354392/1‑149 (MQ=255)                                                                                                       cTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTac                                       >  6:357300/1‑148 (MQ=255)                                                                                                       cTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTac                                       >  1:33080/1‑148 (MQ=255)                                                                                                          aTTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTacaac                                    >  4:160047/1‑148 (MQ=255)                                                                                                          aTTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTacaac                                    >  6:27885/1‑148 (MQ=255)                                                                                                          aTTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTacaac                                    <  2:187109/148‑1 (MQ=255)                                                                                                              ctcGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTACGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATg                                <  3:261261/148‑1 (MQ=255)                                                                                                               tcGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACACCAATgg                               >  2:287570/1‑148 (MQ=255)                                                                                                               tcGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATgg                               >  8:14468/1‑148 (MQ=255)                                                                                                                 gggTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGa                                               <  1:61622/130‑1 (MQ=255)                                                                                                                 gggTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGa                                               >  2:61622/1‑130 (MQ=255)                                                                                                                 gggTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGc                              <  2:45948/147‑1 (MQ=255)                                                                                                                  ggTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGCGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATTGCaaa                           <  5:357300/149‑1 (MQ=255)                                                                                                                   gTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCaaaa                          <  7:268657/149‑1 (MQ=255)                                                                                                                    tAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCaaaaa                         <  5:27885/149‑1 (MQ=255)                                                                                                                      aaGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtc                       <  5:113035/149‑1 (MQ=255)                                                                                                                         tACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtcgtc                    <  4:110490/149‑1 (MQ=255)                                                                                                                           ccGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtcgtcg                   <  4:274753/148‑1 (MQ=255)                                                                                                                           ccGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGa                  <  1:372738/149‑1 (MQ=255)                                                                                                                            cGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTc                                                         >  8:118786/1‑109 (MQ=255)                                                                                                                              cccGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtcgtc                    <  6:247125/144‑1 (MQ=255)                                                                                                                              cccGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtcgtc                    >  5:247125/1‑144 (MQ=255)                                                                                                                              cccGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACAt               >  6:126852/1‑149 (MQ=255)                                                                                                                              cccGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACAt               <  2:72612/149‑1 (MQ=255)                                                                                                                              cccGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAACTTCGCACGGTGAACGGCAAGAAATTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGAATACAACAATCGCAAAAATCGTCGACAt               >  2:279311/1‑149 (MQ=255)                                                                                                                                   tCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAg          <  7:14468/149‑1 (MQ=255)                                                                                                                                    gggagccc‑‑cTGACCGACTGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtcgtc                    <  2:58138/135‑1 (MQ=255)                                                                                                                                    gggagccc‑‑cTGACCGACTGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAAtcgtc                    >  1:58138/6‑140 (MQ=255)                                                                                                                                     c‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGt        >  6:247755/1‑149 (MQ=255)                                                                                                                                          cc‑‑cTGCCCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGa                                               <  4:58351/107‑1 (MQ=255)                                                                                                                                          cc‑‑cTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGa                                               >  3:58351/1‑107 (MQ=255)                                                                                                                                          cc‑‑cTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGt     >  1:283287/1‑149 (MQ=255)                                                                                                                                              cTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGTGa   <  1:377596/149‑1 (MQ=255)                                                                                                                                              cTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGTGa   >  3:307195/1‑149 (MQ=255)                                                                                                                                                 aCCAGTCGGTCAGGGGCTCCCTTCCTCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGTGaa  <  3:243863/147‑1 (MQ=255)                                                                                                                                                 aCCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGTGaa  <  2:283287/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CCATGCGCCCTGCGCCTGGATAAACTGGAGCCGTTCCTGGTCCAGCTCAAGCAATTGGACGACCTGGTGAAAAGTTTTCCGACGGTAGAAACCG‑CGTAAAGCTCATTTCTCGGGTAAAGTACCGCCCGATCC‑‑ACCC‑‑CTGACCAGTCGGTCAGGGGCTCCCTTCATCAGGCCTGCCCAAAGCAAGTTCGCCCGGTGAACGGCAAGAATTTCCTAAGCTAAGTTGTTTTCGTCAAATTCTGGAGTGCTTACAACAATGGCAAAAATCGTCGACATCAAAGGTCGTGAA  >  NC_002947/1808697‑1808982

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