Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,602,236	A→T	15.1%	Q249H (CAA→CAT)	lsrA →	autoinducer 2 import ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,602,236	0	A	T	15.1%	233.0 / 34.3	93	Q249H (CAA→CAT)	lsrA	autoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base T (7/7);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

ACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGT  >  NC_000913/1602138‑1602334                                                                                                   |                                                                                                   aCCATCGCCTTACGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaa                                                                                                   >  1:730225/1‑100 (MQ=255)  ccATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaa                                                                                                   >  1:768977/1‑99 (MQ=255)   cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                 >  2:2721713/1‑100 (MQ=255)       gctttAAGCAGCAAAACCAGCGACCTGTCTACCGACGACATTATTCAGGCCATCACCCCAACAGAACGGAAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                             <  1:1541455/97‑1 (MQ=255)        ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                            <  2:94200/100‑1 (MQ=255)        ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                            <  2:1244170/100‑1 (MQ=255)          ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                                                                          >  2:2279556/1‑100 (MQ=255)             aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCTCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                       <  2:2429111/100‑1 (MQ=255)             aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                       <  2:2102640/100‑1 (MQ=255)             aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGGACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                       <  2:861797/100‑1 (MQ=255)             aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAAAAGTAAGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTgg                                                                                       <  1:88759/100‑1 (MQ=255)              gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                      <  1:2453489/100‑2 (MQ=255)              gCGGCAAAACCAGCGAACTGTCTACAGACGACATAATAAACGCCGTAACCCCAGCGGTACGGGAAAACTCTCACACTGCAAGCCAAAAATTACGGc                                                                                          >  1:1788310/1‑96 (MQ=255)                ggCAAAACCAGCGAACTGTCTACCGACGAGATTATTCAGGCCATCTCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                      <  1:1965651/98‑2 (MQ=255)                ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                    <  1:2289539/100‑4 (MQ=255)                ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                    <  2:879224/100‑4 (MQ=255)                ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                     <  2:2719935/99‑3 (MQ=255)                  cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                  <  1:723508/100‑1 (MQ=255)                  cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                  >  2:860106/1‑100 (MQ=255)                    aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcagagag                                                                                <  2:1822940/100‑8 (MQ=38)                     aaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                <  2:1454783/99‑1 (MQ=255)                           cGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAc                                                                         <  2:2869399/100‑1 (MQ=255)                             aaCTGACTACCGACGACATTATTCAGGCCATCACCCCAGCGGGACGGGAGAAATGGCTCTCTGCCGGCCAAAAATTATGGCTGGAGGTAGCTGGtgagcg                                                                       >  1:2876958/1‑95 (MQ=255)                               cTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGcc                                                                     <  2:646007/100‑1 (MQ=255)                                 gacTACCGACGACATTATTCAGGCCATCACCACAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCACAAATTATGGCTGGAGTTACCTGGTAACCGccta                                                                   >  2:1959867/3‑98 (MQ=255)                                 gTCTACCGGCGACATTATTCCGGCCCTCACCCCAGCGGTAGGGGAGAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                                                                   <  1:400409/100‑1 (MQ=255)                                 gTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                                                                   <  2:1089183/100‑1 (MQ=255)                                 gTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGATAAAAAAAAACCAAAAATTATGGCTGGAGTTACCCGGTAACCGccc                                                                    <  1:2518588/99‑1 (MQ=255)                                  tCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcac                                                                  <  2:2478717/100‑1 (MQ=255)                                   cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcaca                                                                 >  1:2333208/1‑100 (MQ=255)                                   cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcaca                                                                 >  2:201660/1‑100 (MQ=255)                                     aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                               >  2:1742061/1‑100 (MQ=255)                                      ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                              >  1:164135/1‑100 (MQ=255)                                       cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACAt                                                             <  2:2206620/100‑1 (MQ=255)                                       cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACAt                                                             >  1:1117678/1‑100 (MQ=255)                                        nacgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg                                                            >  1:2505102/2‑100 (MQ=255)                                        gacgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg                                                            <  2:1508303/100‑1 (MQ=255)                                          cgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgcc                                                          <  2:1047814/100‑1 (MQ=255)                                             cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                       >  1:2776175/1‑100 (MQ=255)                                             cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                       >  1:2142371/1‑100 (MQ=255)                                               ttattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCgg                                                     >  2:551224/1‑100 (MQ=255)                                               ttattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGCCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCgg                                                     <  1:1791719/100‑1 (MQ=255)                                                tattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGa                                                    <  1:2237517/100‑1 (MQ=255)                                                  ttCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAAc                                                  >  2:182544/1‑100 (MQ=255)                                                    cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACg                                                 <  1:551223/99‑1 (MQ=255)                                                    cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGc                                                >  1:2897780/1‑100 (MQ=255)                                                      ggCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGcc                                               >  1:2523600/1‑99 (MQ=255)                                                      ggCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCg                                              >  1:2521479/1‑100 (MQ=255)                                                        ccATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGt                                            <  1:2079171/100‑1 (MQ=255)                                                         cATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                           <  2:2627226/100‑1 (MQ=255)                                                          aTCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                          <  1:1007908/100‑1 (MQ=255)                                                           tCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCt                                         <  1:830919/100‑1 (MQ=255)                                                           tCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCAGGAACGCCGGTGc                                          >  1:1985141/1‑99 (MQ=255)                                                             aCCCCAGCGGTACGGGAAAAATCGCTCTCTGACAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGa                                       <  1:760642/100‑1 (MQ=255)                                                              ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                      >  2:992200/1‑100 (MQ=255)                                                              ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                      >  1:62086/1‑100 (MQ=255)                                                              ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                      <  1:2728829/100‑1 (MQ=255)                                                               cccAGCAGTACGGGAAAAATCGCTCTCTGACAGCCAAAAATCATGGCTGGAGTTACCTGGTAAACGCCCACAACATGCCGCCGGAACGCCGGTGc                                          <  1:51867/95‑1 (MQ=255)                                                                 cAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACt                                   >  2:1692266/1‑100 (MQ=255)                                                                  aGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                  <  1:975872/100‑1 (MQ=255)                                                                     ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                               <  1:2117405/100‑1 (MQ=255)                                                                      gTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaaa                              >  2:1234366/1‑100 (MQ=255)                                                                        aagggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAAt                            >  2:1960709/3‑100 (MQ=255)                                                                        aCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAAt                            >  2:1692421/1‑100 (MQ=255)                                                                         cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaaaa                             >  1:2320433/1‑98 (MQ=255)                                                                         cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                           >  1:692185/1‑100 (MQ=255)                                                                         cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                           <  1:2896827/100‑1 (MQ=255)                                                                         cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                           >  1:2153350/1‑100 (MQ=255)                                                                          gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                          >  1:824474/1‑100 (MQ=255)                                                                           ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                         >  1:1016012/1‑100 (MQ=255)                                                                                  ccaggtaactccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCg                  >  2:2870730/11‑100 (MQ=37)                                                                                     ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAg               >  1:742245/1‑100 (MQ=255)                                                                                      tctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAgg              >  1:1987610/1‑100 (MQ=255)                                                                                       ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGGAAATCTGACCGGCGAAgg              >  1:623638/1‑99 (MQ=255)                                                                                       ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGt             <  2:1466528/100‑1 (MQ=255)                                                                                       ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGt             <  2:576159/100‑1 (MQ=255)                                                                                         ctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt           >  2:407337/1‑100 (MQ=255)                                                                                          tGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc          >  2:2365811/1‑100 (MQ=255)                                                                                          tGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc          >  2:2365855/1‑100 (MQ=255)                                                                                          ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc          >  2:1034593/3‑100 (MQ=255)                                                                                          ctccagccaTAATTTTTGACTGGAGTTACCTGGTAACAACACACAACATGAGGCCGGAACGCCGGTTCTGACAGTAGAAAATCTGACCGCGGACGGttt           >  1:964747/3‑99 (MQ=255)                                                                                           gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa         <  1:2721703/100‑1 (MQ=255)                                                                                            ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  2:1146807/1‑100 (MQ=255)                                                                                            ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  2:1315628/1‑100 (MQ=255)                                                                                            ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  1:2663287/1‑100 (MQ=255)                                                                                            ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  2:1617231/1‑100 (MQ=255)                                                                                            ccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  1:2089748/1‑100 (MQ=255)                                                                                            ccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        <  2:124526/100‑1 (MQ=255)                                                                                              agccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGa      <  1:592167/100‑1 (MQ=255)                                                                                              agccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGa      >  2:1391362/1‑100 (MQ=255)                                                                                               gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa     <  1:639070/100‑1 (MQ=255)                                                                                               gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa     >  1:1829742/1‑100 (MQ=255)                                                                                               gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGAGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa     >  1:1030155/1‑100 (MQ=255)                                                                                                ccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGACAATCTGACCGGCGAAGGTTTCAGGaa     >  1:393711/1‑99 (MQ=255)                                                                                                 caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg   >  2:1945967/1‑100 (MQ=255)                                                                                                 caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg   >  1:860424/1‑100 (MQ=255)                                                                                                  aaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGt  >  2:2297694/1‑100 (MQ=255)                                                                                                  aaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGt  <  1:1939357/100‑1 (MQ=255)                                                                                                  aaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGt  >  1:1603094/1‑100 (MQ=255)                                                                                                   |                                                                                                   ACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGT  >  NC_000913/1602138‑1602334

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑