Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,370,616	A→G	100%	intergenic (‑204/‑72)	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,616	0	A	G	94.7%	274.9 / 2.5	88	intergenic (‑204/‑72)	yjeH/groS	putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base G (40/43);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATC  >  NC_000913/4370519‑4370712                                                                                                  |                                                                                                 gAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGcc                                                                                                <  1:1756295/100‑1 (MQ=255)   ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAg                                                                                              <  2:605024/100‑1 (MQ=255)    gTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGc                                                                                             <  1:1024416/100‑1 (MQ=255)    gTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGc                                                                                             >  1:2188930/1‑100 (MQ=255)     tGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGcc                                                                                            >  1:612847/1‑100 (MQ=255)        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                         <  1:655638/100‑1 (MQ=255)        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                         >  1:1941917/1‑100 (MQ=255)        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                         <  1:1322994/100‑1 (MQ=255)        aTCAGGGTTTTCACCCGATTGTGTGCTGATCAGAATTTCTCTTCTTTTTCCCCCGTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                         <  2:849122/100‑1 (MQ=255)          cAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGaa                                                                                       <  1:940663/100‑1 (MQ=255)           aGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGaaa                                                                                      >  1:1147859/1‑100 (MQ=255)            gtgTTTTCACCCGATTTTTTGCGGATCAGAATTTTTTTTCTTTTTCTCCCTTGATGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCTGCCGGGaa                                                                                       >  2:1544820/3‑98 (MQ=255)            gggTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAAc                                                                                     <  1:1594070/100‑1 (MQ=255)              gTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCa                                                                                   <  1:1766550/100‑1 (MQ=255)              gTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCa                                                                                   <  1:1766551/100‑1 (MQ=255)               ttttCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCAc                                                                                  >  1:728423/1‑100 (MQ=255)                 ttCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCTGGAAACCACGt                                                                                >  2:2288483/1‑100 (MQ=255)                 ttCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGt                                                                                >  1:1487394/1‑100 (MQ=255)                  tCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTa                                                                               >  1:459999/1‑100 (MQ=255)                     cccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGc                                                                            >  2:1111522/1‑100 (MQ=255)                          ttttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCgn                                                                       <  1:1170479/100‑2 (MQ=255)                          ttttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCgg                                                                       <  1:1419631/100‑1 (MQ=255)                           tttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGc                                                                      <  2:401848/100‑1 (MQ=255)                            ttGTGCTGATCAGAATTTTTGTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTTGCCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                     <  1:82413/100‑1 (MQ=255)                              gtgCTGATATGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGc                                                                      <  2:1968473/97‑1 (MQ=255)                               tgCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGc                                                                      <  2:116864/96‑1 (MQ=255)                                 cTGATCAGAGTTTTTTTTCTTTTTCCCACTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCAcc                                                                >  2:994885/1‑100 (MQ=255)                                 cTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCAcc                                                                >  1:2153270/1‑100 (MQ=255)                                   gATCAGAATTTTTTTTTTTTTTTTTCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTTTTCGCTAGCCGGGAAACCACGTAAGCTCCGGCGTCa                                                                  >  2:7541/1‑96 (MQ=255)                                      cAGAATTTTTTTTTTTTTTTTTTTTTTAAGGGGCGAAGCCTCATCCCCATTTCTCTTGTCACCAGCCGGGTACCCACGTAAGCTCCGGCGTCACCCATa                                                            >  1:1345806/1‑99 (MQ=255)                                      cAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATaa                                                           <  1:1794223/100‑1 (MQ=255)                                        gAATTTTTTTTTTTTTTTCCCCTTGAAAGGGCGAAGCCTCATCCCTATTTCTCTTGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACTCATaa                                                           >  2:1191259/1‑98 (MQ=255)                                        gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTTGCCTGTCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                         >  2:445987/1‑100 (MQ=255)                                        gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                         <  1:23830/100‑1 (MQ=255)                                        gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                         <  1:615463/100‑1 (MQ=255)                                         aaTTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACACATAACAg                                                        <  2:501943/100‑1 (MQ=255)                                              tttttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAATCTCCGGCGTCACCCATAACAGATACg                                                   <  2:371269/100‑1 (MQ=255)                                                tttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACgg                                                  >  1:393465/1‑99 (MQ=255)                                                tttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGa                                                 <  1:1456918/100‑1 (MQ=255)                                                 ttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGAc                                                <  1:884585/100‑1 (MQ=255)                                                  tctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACt                                               >  2:182438/1‑100 (MQ=255)                                                   ctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACCGATACGGACtt                                              >  1:791701/1‑100 (MQ=255)                                                    tttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACttt                                             >  1:173193/1‑100 (MQ=255)                                                     ttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                            <  1:1111527/100‑1 (MQ=255)                                                     ttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                            >  1:581701/1‑100 (MQ=255)                                                       ttcCCCTTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtct                                           <  2:1951873/99‑1 (MQ=255)                                                         cccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                        >  1:2221023/1‑100 (MQ=255)                                                          ccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaaa                                       <  1:1926553/100‑1 (MQ=255)                                                           cccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                      <  1:667244/100‑1 (MQ=255)                                                            ccgtGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAgg                                     <  2:2184265/97‑1 (MQ=255)                                                              ttGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAgg                                     <  2:980173/98‑1 (MQ=255)                                                                gAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                 <  1:1692921/100‑1 (MQ=255)                                                                 aaGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGt                                <  2:1682715/100‑1 (MQ=255)                                                                   ggggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTa                              <  2:547353/100‑1 (MQ=255)                                                                   ggggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTa                              >  1:44685/1‑100 (MQ=255)                                                                    gggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTAt                             >  2:314155/1‑100 (MQ=255)                                                                    gggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTAt                             <  1:766525/100‑1 (MQ=255)                                                                     ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                            >  2:139944/1‑100 (MQ=255)                                                                     ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCAGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                            <  2:123356/100‑1 (MQ=255)                                                                     ggCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                            >  1:1691650/1‑100 (MQ=255)                                                                      gCGAAGCCTCATCCCCATTTCTCTGGTCGCCCGCCGGGGAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                            <  2:1676561/99‑1 (MQ=255)                                                                      gCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCa                           >  1:352570/1‑100 (MQ=255)                                                                       cGACGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCaa                          >  1:622281/1‑100 (MQ=255)                                                                        gAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAAt                         >  1:59182/1‑100 (MQ=255)                                                                        gAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAAt                         <  1:1176600/100‑1 (MQ=255)                                                                        gAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAAt                         >  1:1144489/1‑100 (MQ=255)                                                                         aaGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATg                        <  2:934440/100‑1 (MQ=255)                                                                           gCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                      <  1:715605/100‑1 (MQ=255)                                                                           gCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                      >  1:676702/1‑100 (MQ=255)                                                                            cgtCATCCCCATTCCTCTGGTCGCCAGCCGGAAAACCACGTATGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAat                     >  1:1609362/3‑100 (MQ=255)                                                                             cTCATCCCCATTTCTCTGGTCGCCAGCCGGGACCCCACTTAAGCTCCCGCCTCACCCATCACAGATACGACCTTTCTCACCGGAGAGATATCAATGaa                      >  2:1447400/1‑98 (MQ=255)                                                                             cTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAata                    >  2:658527/1‑100 (MQ=255)                                                                                   cccATTTCTCTGGTCGCCGGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              >  2:2136413/1‑100 (MQ=255)                                                                                   cccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              >  2:1721871/1‑100 (MQ=255)                                                                                   cccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              >  1:1023974/1‑100 (MQ=255)                                                                                   cccATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              <  2:676700/100‑1 (MQ=255)                                                                                    ccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTcc             <  1:1270354/100‑1 (MQ=255)                                                                                    ccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTAAAAGGAGAGTTATCAATGAATATTCGTcc             >  2:704909/1‑100 (MQ=255)                                                                                       tttCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGt               <  2:1965598/95‑1 (MQ=255)                                                                                       tttCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAAAATTCGTCCAtt          >  1:2056464/1‑100 (MQ=255)                                                                                        ttCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTg         <  1:1212988/100‑1 (MQ=255)                                                                                          ctctGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCa       <  1:256033/100‑1 (MQ=255)                                                                                           tctGGTCGCCAGCCGGGAAAGCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTCTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCAt      <  2:828631/100‑1 (MQ=255)                                                                                            ctGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATg     >  1:1126454/1‑100 (MQ=255)                                                                                            ctGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATg     >  1:55074/1‑100 (MQ=255)                                                                                              ggTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt   >  2:1104513/1‑100 (MQ=255)                                                                                              ggTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAAAAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAAAATTCGACCATTGCATg     >  1:2344023/1‑98 (MQ=255)                                                                                               gTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATc  >  1:1367257/1‑100 (MQ=255)                                                                                                  |                                                                                                 GAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATC  >  NC_000913/4370519‑4370712

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