Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,802	A→G	100%	V222A (GTA→GCA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,802	0	A	G	98.0%	327.4 / ‑4.8	96	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (42/51);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGC  >  NC_000913/953705‑953893                                                                                                  |                                                                                            gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCa                                                                                           <  2:2885918/100‑1 (MQ=255)   tgaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                         >  2:2535496/6‑100 (MQ=255)      ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                      >  1:346787/3‑100 (MQ=255)      ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                       >  2:1876203/3‑99 (MQ=255)      cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGGCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                      >  1:2908488/1‑100 (MQ=255)        gTGATGAAATTCATCACGGTCAGGTGAG‑AAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                   >  2:441488/1‑100 (MQ=255)         tgatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                   >  2:526135/1‑100 (MQ=255)          gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                  >  1:2459968/1‑100 (MQ=255)          gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                  >  2:681585/1‑100 (MQ=255)           atgaAATTCATCACGGTCAGGTGAGAAAAATTTTACGGTGCAAAACCTACTGCGGTCCAAAATTTCGGGGATGCGAAGTCGCGGATAACAATAc                                                                                       >  1:2340735/1‑94 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                <  2:2738246/100‑1 (MQ=255)               aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                             <  1:1536848/100‑1 (MQ=255)               aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                             <  2:934548/100‑1 (MQ=255)               aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                             <  2:1990013/100‑1 (MQ=255)                   catcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                         <  2:1088100/100‑1 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        >  2:650427/1‑100 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        <  2:2676745/100‑1 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        <  2:837652/100‑1 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        >  2:564300/1‑100 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        <  2:610530/100‑1 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATATTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        >  1:1689635/1‑100 (MQ=255)                     tcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaaa                                                                       <  2:65025/100‑1 (MQ=255)                        cGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACAt                                                                    <  1:1699561/100‑1 (MQ=255)                        cGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACAt                                                                    >  2:2763628/1‑100 (MQ=255)                              ggcgagaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                              >  1:2834212/4‑100 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACagtacggt                                                              >  1:719099/1‑93 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                              >  1:1080968/1‑100 (MQ=255)                               gTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                             <  1:1891915/100‑1 (MQ=255)                                  agaAAAATTTTCCGGTGCAGAAACGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                           <  2:1222067/99‑1 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGCTTGCAATACCCATCGGGATCATAACCACGTTTGCGATAc                                                         >  1:2286664/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                         >  1:2477123/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                         >  2:1238398/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                         >  1:517217/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                         >  1:2398428/1‑100 (MQ=255)                                    aaaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACt                                                        <  1:1490851/100‑1 (MQ=255)                                     aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                       <  1:2877299/100‑1 (MQ=255)                                     aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                       >  2:2528474/1‑100 (MQ=255)                                     aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                       >  2:2528653/1‑100 (MQ=255)                                        aTTTTCCGGTGCAGAACCGGCTGCGGGCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                    <  1:1181441/100‑1 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                    <  1:747989/100‑1 (MQ=255)                                         ttttCCGGTGCAGAACCGACTGCAGTCCAAAATTCCGGGGATGCGAAGTCTCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                    <  1:1964293/99‑1 (MQ=255)                                            tCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaa                                                <  1:2029531/100‑1 (MQ=255)                                             ccGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                                               <  1:326524/100‑1 (MQ=255)                                              cGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaaa                                              >  2:121360/1‑100 (MQ=255)                                              cGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaaa                                              >  2:1586161/1‑100 (MQ=255)                                               ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                             >  2:2353322/1‑100 (MQ=255)                                                gTGCAGAACCGACTGTGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAcc                                            <  2:101300/100‑1 (MQ=255)                                                gTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAcc                                            >  2:885976/1‑100 (MQ=255)                                                   cAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                         <  2:2763151/100‑1 (MQ=255)                                                   cAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                         <  1:2381412/100‑1 (MQ=255)                                                   cAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCACGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                         >  1:147011/1‑100 (MQ=255)                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                        >  2:2813153/1‑100 (MQ=255)                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                        <  2:1888777/100‑1 (MQ=255)                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                        >  1:1146007/1‑100 (MQ=255)                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                        <  2:1657493/100‑1 (MQ=255)                                                     gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                       <  1:1680487/100‑1 (MQ=255)                                                     gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                       <  2:819463/100‑1 (MQ=255)                                                     gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAACACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                       <  2:653263/100‑1 (MQ=255)                                                          gACTGCGGTTCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGTGATCCTGTGCTCAAAACCGCTGGcaaca                                  <  1:2544392/100‑1 (MQ=255)                                                          gACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaaca                                  <  2:2188980/100‑1 (MQ=255)                                                             tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                               <  1:2235139/100‑1 (MQ=255)                                                             tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                               <  2:650021/100‑1 (MQ=255)                                                              gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                              >  1:1465580/1‑100 (MQ=255)                                                              gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAAGACCCAGCGGGATCATAAACATTTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                              <  2:277051/100‑1 (MQ=255)                                                                ggTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATc                            >  2:2563939/1‑100 (MQ=255)                                                                  tCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgc                          <  1:364234/100‑1 (MQ=255)                                                                   ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                         <  1:2109782/100‑1 (MQ=255)                                                                   ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAGACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                         <  1:1420618/100‑1 (MQ=255)                                                                   ccAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                         <  2:1155043/100‑1 (MQ=255)                                                                    cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                         <  2:2516439/99‑1 (MQ=255)                                                                    cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                        <  2:916809/100‑1 (MQ=255)                                                                    cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTACTCAAAACCGCTGGCAACAAACATcgcg                         >  1:1271904/1‑99 (MQ=255)                                                                     aaaaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGAc                       >  1:1988372/1‑100 (MQ=255)                                                                      aaaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGAcc                      <  1:2034314/100‑1 (MQ=255)                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                   <  1:1453956/100‑1 (MQ=255)                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                   <  2:2398439/100‑1 (MQ=255)                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                   <  2:105184/100‑1 (MQ=255)                                                                          tCCGGGGATGCGAAGTCGCGGATTGTAATACCCATCCGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                  <  1:1932136/100‑1 (MQ=255)                                                                           ccGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcag                 >  1:2841408/1‑100 (MQ=255)                                                                           ccGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcag                 >  1:548868/1‑100 (MQ=255)                                                                           ccGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgaag                 <  1:123643/100‑1 (MQ=255)                                                                            cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                <  2:2456337/100‑1 (MQ=255)                                                                             ggggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagca               <  2:511404/100‑1 (MQ=255)                                                                             ggggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagca               >  1:1612106/1‑100 (MQ=255)                                                                              gggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCAc              <  2:2374980/100‑1 (MQ=255)                                                                                   gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa         >  1:2823390/1‑100 (MQ=255)                                                                                    cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt        <  1:2763618/100‑1 (MQ=255)                                                                                    cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt        >  2:825783/1‑100 (MQ=255)                                                                                    cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt        >  2:825884/1‑100 (MQ=255)                                                                                     gAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg       <  2:1756949/100‑1 (MQ=255)                                                                                      aaGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGa      <  2:2841418/100‑1 (MQ=255)                                                                                      aaGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGa      <  1:85622/100‑1 (MQ=255)                                                                                        gTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAAc    >  1:49317/1‑100 (MQ=255)                                                                                        gTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACAGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAAc    >  2:1123021/1‑100 (MQ=255)                                                                                         tCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACg   <  2:2687306/100‑1 (MQ=255)                                                                                          cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc  >  2:1947479/1‑100 (MQ=255)                                                                                                  |                                                                                            GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGC  >  NC_000913/953705‑953893

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