Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,370,614	T→A	intergenic (‑202/‑74)	yjeH ← / → groS	putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,370,614	0	T	A	96.0%	324.9 / 0.2	100	intergenic (‑202/‑74)	yjeH/groS	putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (45/51);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAT  >  NC_000913/4370519‑4370711                                                                                                 |                                                                                                  gAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACAcc                                                                                               <  2:430448/100‑1 (MQ=255) gAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACAcc                                                                                               <  1:387761/100‑1 (MQ=255) gAGGCGAATGAGGGTTTTCCCCCGATTTTGGGCTGATCAGAATTTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATGTCTCTGGACa                                                                                                 <  1:359754/99‑1 (MQ=255)   ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTTTTTTTTTTCCTTGAAGGGGCGAAGCCTCATCTCCATTTCTCTGGAcagcag                                                                                             >  2:281838/1‑100 (MQ=255)   ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAg                                                                                             >  1:888796/1‑100 (MQ=255)   ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAg                                                                                             >  1:8047/1‑100 (MQ=255)   ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATAAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAg                                                                                             >  1:502138/1‑100 (MQ=255)      gAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCg                                                                                          >  1:779045/1‑100 (MQ=255)      gAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCg                                                                                          >  1:1830938/1‑100 (MQ=255)       aaTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCgg                                                                                         <  1:424905/100‑1 (MQ=255)        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCggg                                                                                        >  2:1506416/1‑100 (MQ=255)        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCggg                                                                                        >  1:2110614/1‑100 (MQ=255)             tgTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAAcc                                                                                   >  1:221773/2‑100 (MQ=255)              gTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCa                                                                                  <  1:1692669/100‑1 (MQ=255)              gTTTTCACCCGATTCTGTGCTGATCAGAA‑TTTTTTTGCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATT‑‑‑‑TGGTCACCAGCCGGGAAACCACGTa                                                                              <  1:1195458/100‑1 (MQ=255)                tttCACCCGATTTTGTGCTGATCAGAATTTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCAc                                                                                 <  2:182854/100‑1 (MQ=255)                 ttCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGt                                                                               <  1:2581265/100‑1 (MQ=255)                  tCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTa                                                                              <  1:297470/100‑1 (MQ=255)                    aCCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAg                                                                            <  2:826491/100‑1 (MQ=255)                     cccGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGc                                                                           >  1:540782/1‑100 (MQ=255)                      ccGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCt                                                                          <  1:1570357/100‑1 (MQ=255)                      ccGATTTAGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCt                                                                          >  1:1419563/1‑100 (MQ=255)                       cGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTc                                                                         >  2:234623/1‑100 (MQ=255)                         aTTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCg                                                                       >  1:2427776/1‑100 (MQ=255)                            ttGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                    <  1:1497475/100‑1 (MQ=255)                             tgtgCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGgtgt                                                                   >  2:194593/1‑100 (MQ=255)                              gtgCTGATCAGAA‑TTTTTTTTCTTTTTCCTCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTc                                                                  >  1:910886/1‑100 (MQ=255)                              gtgCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTc                                                                  >  1:971787/1‑100 (MQ=255)                              gtgCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTc                                                                  <  2:842475/100‑1 (MQ=255)                               tgCTGATCAGAATTTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACCCCACCCGGGAAACCCCGTAAGCCCCg                                                                       >  1:2494923/1‑95 (MQ=255)                                gCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTTGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAc                                                                >  2:1827494/1‑100 (MQ=255)                                 cTGATCAGAA‑TTTTTTTTTTTTTTTTTTCTTGAAGAAGCGAAGCCTCATCCTCATTTTTTTGCTCACCAGCCGGGAAACCACTTAAGCTCCTGCGTCAcc                                                               >  2:820013/1‑100 (MQ=255)                                 cTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAcc                                                               >  2:819234/1‑100 (MQ=255)                                 cTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAcc                                                               <  2:1450663/100‑1 (MQ=255)                                 cTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAcc                                                               <  2:1679751/100‑1 (MQ=255)                                  tGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAccc                                                              <  2:1361458/100‑1 (MQ=255)                                     tCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATa                                                           <  1:2471124/100‑1 (MQ=255)                                     tCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATa                                                           >  2:1623714/1‑100 (MQ=255)                                       aGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAAc                                                         <  1:1238487/100‑1 (MQ=255)                                       aGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAAc                                                         <  1:1191363/100‑1 (MQ=255)                                       aGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAAc                                                         <  1:1928960/100‑1 (MQ=255)                                       aGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAAc                                                         <  1:1930443/100‑1 (MQ=255)                                        gAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                        <  2:1119205/100‑1 (MQ=255)                                        aaa‑TTTTTTTTCTTGTTCCCCCTTGAAGGGGCGCAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAccc                                                              <  1:995674/93‑1 (MQ=255)                                          a‑TTTTTTTTCTTTTTCCCCCATGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTNCGGCGTCACCCATAACAGa                                                      <  1:1472799/100‑1 (MQ=255)                                            ttttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCTATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGa                                                      >  1:286913/1‑99 (MQ=255)                                            ttttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGAt                                                     <  1:1134807/100‑1 (MQ=255)                                             tttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATa                                                    <  1:2430422/100‑1 (MQ=255)                                              ttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATAc                                                   >  2:1303344/1‑100 (MQ=255)                                                    ctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACtt                                             >  2:415348/1‑100 (MQ=255)                                                     tttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACttt                                            >  1:852082/1‑100 (MQ=255)                                                      ttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                           <  2:123991/100‑1 (MQ=255)                                                       tttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtct                                          <  1:1370544/100‑1 (MQ=255)                                                       ttccccccTTGAAGGGGCGAGGCCTCATCCCCATTCCTCTGGACACCAGCCGGGAATCCACGTAAGCTCCGGCGTCACCCATACCAGATACGGACTTtc                                           <  2:1867711/96‑1 (MQ=255)                                                         tcCCCCTAGAAAGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCa                                        >  2:1412404/1‑100 (MQ=255)                                                          cccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                       <  2:2494898/100‑1 (MQ=255)                                                          cccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                       <  2:1397882/100‑1 (MQ=255)                                                           ccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaaa                                      >  1:997463/1‑100 (MQ=255)                                                            cccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                     <  1:1756114/100‑1 (MQ=255)                                                              cTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGga                                   >  2:1651028/1‑100 (MQ=255)                                                               ttGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgag                                  <  1:234625/100‑1 (MQ=255)                                                               ttGAAGGGGCGAAGCCTCAGCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCGATAACAGATACGGACTTTCTCAAAGgag                                  <  2:751428/100‑1 (MQ=255)                                                                ttaaGGGGCGAAGGCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                 >  1:1824984/3‑100 (MQ=255)                                                                tGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                 >  1:177000/1‑100 (MQ=255)                                                                tGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                 <  1:214831/100‑1 (MQ=255)                                                                 gAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                <  2:971779/100‑1 (MQ=255)                                                                 gAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                <  2:363427/100‑1 (MQ=255)                                                                 gAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                <  1:1056352/100‑1 (MQ=255)                                                                 gAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                <  1:2090403/100‑1 (MQ=255)                                                                  aaGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGt                               <  2:779038/100‑1 (MQ=255)                                                                      ggCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           >  1:695187/1‑100 (MQ=255)                                                                      ggCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           <  1:1713430/100‑1 (MQ=255)                                                                      ggCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           >  1:65536/1‑100 (MQ=255)                                                                       gCGAAGCCTCATCCCCCTTTCTCTGGACCCCAGACGCAACACCACCTAAGCTCCGGCGTCACCAATAACAGATACGGACTTTCTCAAAGGAGATTTAtn                           <  2:221772/99‑2 (MQ=255)                                                                       gCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCa                          <  1:421662/100‑1 (MQ=255)                                                                          aaGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATg                       <  2:502131/100‑1 (MQ=255)                                                                           aGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGACACGGACTTGCTCAAAGGAGAGTTATCAATGa                      >  2:1811267/1‑100 (MQ=255)                                                                           aGCCTCATCCCCATTTCTCTGGACACCAACCGGGAAACCACGTAAGCTCCGGCGTAACGCATAACAGATAAGGAATTTATCAAAGGGGAGttcataatta                      >  2:2438255/1‑92 (MQ=255)                                                                            gCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                     >  1:196141/1‑100 (MQ=255)                                                                            gCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                     <  1:2454852/100‑1 (MQ=255)                                                                            gCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                     <  2:274678/100‑1 (MQ=255)                                                                             ccTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                     <  2:117681/99‑1 (MQ=255)                                                                                cATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATAtt                 <  1:405791/100‑1 (MQ=255)                                                                                cATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATAtt                 <  1:159575/100‑1 (MQ=255)                                                                                  tCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg               >  2:2024249/1‑100 (MQ=255)                                                                                  tCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg               >  1:455342/1‑100 (MQ=255)                                                                                  tCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg               >  1:805758/1‑100 (MQ=255)                                                                                    cccATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGt              >  2:1274071/1‑99 (MQ=255)                                                                                    cccATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGt              >  2:2079385/1‑99 (MQ=255)                                                                                    cccATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc             >  2:678378/1‑100 (MQ=255)                                                                                    cccATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc             >  1:573587/1‑100 (MQ=255)                                                                                    cccATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc             >  1:680106/1‑100 (MQ=255)                                                                                      cATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCa           <  2:2110592/100‑1 (MQ=255)                                                                                         ttCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTg        <  2:1877440/100‑1 (MQ=255)                                                                                          tctcTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCGCAGAGGAGAGTTATCAATGAATATTCGTCCAtt         >  2:931701/1‑98 (MQ=255)                                                                                           ctctGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCa      >  1:383448/1‑100 (MQ=255)                                                                                           ctctGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCa      >  1:383368/1‑100 (MQ=255)                                                                                            tctGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCAt     <  1:1078770/100‑1 (MQ=255)                                                                                              tGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGa   >  1:1571914/1‑100 (MQ=255)                                                                                               ggTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt  >  1:2102124/1‑100 (MQ=255)                                                                                                 |                                                                                                  GAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAA‑TTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAT  >  NC_000913/4370519‑4370711

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