Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	383,477	G→A	100%	G153G (GGC→GGT)	yaiP ←	putative family 2 glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	383,477	0	G	A	95.8%	299.5 / 2.3	96	G153G (GGC→GGT)	yaiP	putative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (50/42);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

TAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACA  >  NC_000913/383380‑383573                                                                                                  |                                                                                                 taataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCAcc                                                                                                >  1:360449/1‑100 (MQ=255)   ataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCtt                                                                                              >  2:2089156/1‑100 (MQ=255)   ataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCtt                                                                                              >  2:1963021/1‑100 (MQ=255)     aaaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTc                                                                                            >  1:2264408/1‑100 (MQ=255)     aaaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTc                                                                                            >  2:1889541/1‑100 (MQ=255)       aCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCaa                                                                                          <  1:649948/100‑1 (MQ=255)         ggTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACg                                                                                        <  1:925941/100‑1 (MQ=255)         ggTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACg                                                                                        <  2:637423/100‑1 (MQ=255)            gcgcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAg                                                                                     >  1:1707332/1‑100 (MQ=255)             cgcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACACTCATCGGCACCTTTACGGTCGCGCGGATGTCCGGTAACACACCCGCACCTTTCAACGCAGt                                                                                    >  2:1084291/1‑100 (MQ=255)             cgcCACCCACGAGCTGCTGTAGCTTGCGCCTCATAACAATCATCGGCAACTTTACGTTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGt                                                                                    <  1:2136971/100‑1 (MQ=255)              gcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTa                                                                                   >  2:1277675/1‑100 (MQ=255)               acACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTag                                                                                  >  1:930147/2‑100 (MQ=255)                cACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGAGGTAACAGACCCGCACCTTTCAACGCAGTaga                                                                                 >  2:93515/1‑100 (MQ=255)                   ccAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAggg                                                                              <  2:1050345/100‑1 (MQ=255)                    cAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGa                                                                             <  1:2089170/100‑1 (MQ=255)                    cAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGa                                                                             >  2:256405/1‑100 (MQ=255)                     aGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTCGAGGGaa                                                                            >  2:725465/1‑100 (MQ=255)                       gAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATg                                                                          <  2:1024355/100‑1 (MQ=255)                        aGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATgc                                                                         <  2:1332475/100‑1 (MQ=255)                         gctgctGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATgcc                                                                        <  1:2335784/100‑1 (MQ=255)                          ctgctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATgccg                                                                       <  2:1707320/100‑1 (MQ=255)                            gctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATgccgcc                                                                     <  2:2389690/100‑1 (MQ=255)                             ctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCt                                                                    <  2:1434900/100‑1 (MQ=255)                             ctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCt                                                                    <  2:1831224/100‑1 (MQ=255)                               gTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTAc                                                                  <  2:81001/100‑1 (MQ=255)                               gTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTAc                                                                  >  1:286812/1‑100 (MQ=255)                                 aGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACgg                                                                >  1:1297352/1‑100 (MQ=255)                                 aGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATACCGCCTACg                                                                 >  2:362658/1‑99 (MQ=255)                                 aGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAAACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACgg                                                                >  1:594446/1‑100 (MQ=255)                                   cGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCa                                                              <  2:664107/100‑1 (MQ=255)                                   cGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCAGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACgg                                                                <  2:2102997/98‑1 (MQ=255)                                     tGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATc                                                            <  2:895226/100‑1 (MQ=255)                                      gcgcTTCATAACAATCACCGGCCACTTTACGGTCGCGCGGATGTGCGGCAACAGACCCGCACCTTTCATCGCAGTATAGGGAATGCCGCCTACGGCATCg                                                           <  2:2341537/100‑1 (MQ=255)                                       cgcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCgg                                                          >  2:2010990/1‑100 (MQ=255)                                        gcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGc                                                         >  1:719472/1‑100 (MQ=255)                                        gcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGc                                                         >  1:1522471/1‑100 (MQ=255)                                           tCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCAcc                                                      <  1:2098193/100‑1 (MQ=255)                                           tCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCAcc                                                      <  2:285368/100‑1 (MQ=255)                                            cATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcg                                                     >  2:1693787/1‑100 (MQ=255)                                            cATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcg                                                     >  2:191634/1‑100 (MQ=255)                                             aTAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcgc                                                    <  1:344085/100‑1 (MQ=255)                                              tAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcgcg                                                   >  1:656526/1‑100 (MQ=255)                                              tAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcgcg                                                   >  2:739941/1‑100 (MQ=255)                                              tAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcgcg                                                   >  2:2342371/1‑100 (MQ=255)                                              tAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcgcg                                                   <  2:2264394/100‑1 (MQ=255)                                               aaCAATCATCGGCAAGTTTACGGTCGCGCGGATATAAGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCg                                                           >  2:1954474/1‑91 (MQ=255)                                                aCAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCt                                                 >  2:850188/1‑100 (MQ=255)                                                    tcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAt                                             >  1:363122/1‑100 (MQ=255)                                                    tcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAt                                             <  2:957689/100‑1 (MQ=255)                                                    tcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAt                                             >  1:1892237/1‑100 (MQ=255)                                                      atcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAttt                                           <  2:125261/100‑1 (MQ=255)                                                      atcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAttt                                           >  2:167329/1‑100 (MQ=255)                                                       tcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTc                                          <  1:167331/100‑1 (MQ=255)                                                       tcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTc                                          <  2:1882866/100‑1 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        <  1:711872/100‑1 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        <  2:2291111/100‑1 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        >  2:1747649/1‑100 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        <  2:1830677/100‑1 (MQ=255)                                                            aaCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCa                                     <  1:850195/100‑1 (MQ=255)                                                             aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                    <  1:1671309/100‑1 (MQ=255)                                                             aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                    >  2:427631/1‑100 (MQ=255)                                                              cTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGc                                   >  2:2464544/1‑100 (MQ=255)                                                                   cGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCAtata                              <  2:534010/100‑1 (MQ=255)                                                                   cGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCAtata                              >  1:495727/1‑100 (MQ=255)                                                                   cGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCAtata                              >  1:289959/1‑100 (MQ=255)                                                                    ggTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAg                             >  1:1663166/1‑100 (MQ=255)                                                                     gTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCGTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGc                            <  2:1402975/100‑1 (MQ=255)                                                                       cgcgcgGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGccc                          >  1:181786/1‑100 (MQ=255)                                                                        gcgcgGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCa                         >  1:586470/1‑100 (MQ=255)                                                                         cgcgGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCAt                        >  2:2114164/1‑100 (MQ=255)                                                                         cgcgGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCAt                        <  1:508583/100‑1 (MQ=255)                                                                         cgcgGATGTGCGGTAACAGACCCGCACCTTGCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCGTATAGccc                          >  2:1686571/1‑98 (MQ=255)                                                                               tgtgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTc                  <  2:2196937/100‑1 (MQ=255)                                                                               tgtgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTc                  <  2:1358446/100‑1 (MQ=255)                                                                                gtgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCt                 <  2:43554/100‑1 (MQ=255)                                                                                 tgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCtt                >  1:544425/1‑100 (MQ=255)                                                                                 tgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCtt                >  1:2514053/1‑100 (MQ=255)                                                                                 tgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCt                 >  2:1213147/1‑99 (MQ=255)                                                                                  gCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg               <  1:178025/100‑1 (MQ=255)                                                                                  gCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg               <  1:826165/100‑1 (MQ=255)                                                                                   cGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGa              >  2:919622/1‑100 (MQ=255)                                                                                    ggTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGAt             >  2:1715805/1‑100 (MQ=255)                                                                                    ggTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGAt             >  2:1443665/1‑100 (MQ=255)                                                                                     gTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATc            <  1:1747132/100‑1 (MQ=255)                                                                                     gTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATc            >  1:154721/1‑100 (MQ=255)                                                                                      tAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCg           <  1:1776703/100‑1 (MQ=255)                                                                                        aCAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCggg         <  2:194607/100‑1 (MQ=255)                                                                                          aGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCg       <  1:2158734/100‑1 (MQ=255)                                                                                           gACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCgg      >  1:2103174/1‑100 (MQ=255)                                                                                            aCCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGa     >  2:1268448/1‑100 (MQ=255)                                                                                             cccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa    >  2:1276130/1‑100 (MQ=255)                                                                                             cccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa    >  1:1229449/1‑100 (MQ=255)                                                                                             cccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa    >  2:1597655/1‑100 (MQ=255)                                                                                              ccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAAc   <  2:1522460/100‑1 (MQ=255)                                                                                              ccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGCAAc   <  2:754133/100‑1 (MQ=255)                                                                                               cgcaCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACa  >  1:1422405/1‑100 (MQ=255)                                                                                                  |                                                                                                 TAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACA  >  NC_000913/383380‑383573

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