Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,427 C→A 5.5% intergenic (+11/+26) fhlA → / ← ygbA formate hydrogenlyase transcriptional activator/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,4270CA 5.5% 240.7 / 3.7 73intergenic (+11/+26)fhlA/ygbAformate hydrogenlyase transcriptional activator/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base A (2/2);  total (33/40)
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.71e-01

TGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCC  >  NC_000913/2856280‑2856570
                                                                                                                                                   |                                                                                                                                               
tGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCTAATTGc                                                                                                                                                 >  2:393067/1‑148 (MQ=255)
 gATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGc                                                                                                                                                 <  8:425377/147‑1 (MQ=255)
  aTTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCg                                                                                                                                               >  2:44196/1‑148 (MQ=255)
   ttGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCgg                                                                                                                                              <  5:224848/148‑1 (MQ=255)
   ttGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGCGAATTGATAAATGGGCATTGATTTAACTGCAAATTGCCgg                                                                                                                                              <  3:219340/148‑1 (MQ=255)
            tGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTg                                                                                                                                     <  1:257980/148‑1 (MQ=255)
               tGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGAACTGCCt                                                                                                                                  >  4:309309/1‑148 (MQ=255)
                gaaagaaaCCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTg                                                                                                                                 >  3:32952/1‑148 (MQ=255)
                gaaagaaaCCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTg                                                                                                                                 >  7:236881/1‑148 (MQ=255)
                gaaagaaaCCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTg                                                                                                                                 >  7:356531/1‑148 (MQ=255)
                        ccAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGTGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCAGTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCAt                                                                                                                         <  3:94317/148‑1 (MQ=255)
                         cAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGAACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATa                                                                                                                        >  2:395858/1‑148 (MQ=255)
                             ggCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTAt                                                                                                                    <  3:82790/148‑1 (MQ=255)
                                gTGGTTGCCGGGCCTAAAGGCGCGGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCa                                                                                                                 <  7:6187/148‑1 (MQ=255)
                                  ggTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGCCAGATCTGCCTGTCCGGCATACTATTCAt                                                                                                                >  3:296491/1‑147 (MQ=255)
                                  ggTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATg                                                                                                               <  1:201246/148‑1 (MQ=255)
                                   gTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGa                                                                                                              <  7:74239/148‑1 (MQ=255)
                                   gTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGa                                                                                                              <  8:386953/148‑1 (MQ=255)
                                               ataGGGGCTGCGCAACGTCTGGGGCTGAAACGCTCGACCCTGCTGTCACGGATGAAGCGGCTGGGTATTGATAAATTGTCATTGATTTTTCTGCAAATTGCCGGa                                                                                                                                             <  4:62721/103‑1 (MQ=255)
                                               aaaGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGa                                                                                                  >  7:98077/1‑148 (MQ=255)
                                                  ggCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGa                                                                                               >  3:335417/1‑148 (MQ=255)
                                                   gcgcTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGat                                                                                              <  1:173449/148‑1 (MQ=255)
                                                         cgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGAt                                                                                                                                        <  2:304668/100‑1 (MQ=255)
                                                         cgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGAt                                                                                                                                        >  1:304668/1‑100 (MQ=255)
                                                         cgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATAttttt                                                                                        >  5:145795/1‑148 (MQ=255)
                                                            aaCGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCg                                                                                     >  1:120652/1‑148 (MQ=255)
                                                             aCGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                    >  4:291418/1‑148 (MQ=255)
                                                              cGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                   <  3:156050/148‑1 (MQ=255)
                                                                    gggCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTc                                                                              <  3:291418/147‑1 (MQ=255)
                                                                          aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAcc                                                                       >  7:239326/1‑148 (MQ=255)
                                                                          aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAcc                                                                       >  1:375221/1‑148 (MQ=255)
                                                                          aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAcc                                                                       >  4:252403/1‑148 (MQ=255)
                                                                                 cGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCt                                                                <  4:32952/148‑1 (MQ=255)
                                                                                      ctgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCa                                                           >  1:353529/1‑148 (MQ=255)
                                                                                         ctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATc                                                         <  3:309309/147‑1 (MQ=255)
                                                                                         ctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTAGGACGATAATTTTCAGGCAGTTATGGCACCTGACGCTTGTCa                                                           >  7:427291/1‑145 (MQ=255)
                                                                                             cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGa                                                                                                                                         >  8:189163/1‑63 (MQ=255)
                                                                                             cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGa                                                                                                                                         <  7:189163/63‑1 (MQ=255)
                                                                                              aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgag                                                   <  4:232513/148‑1 (MQ=255)
                                                                                               cGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgat                                                      >  8:133342/1‑144 (MQ=255)
                                                                                                ggATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat                                                 <  7:112593/148‑1 (MQ=255)
                                                                                                      tgCGGCTTGTAATTGATAATTCGGCATTGTTTTAACTTCAAATTGCCGTACAGATCAGCCTGTCAGGCATACTATTCTTGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                            <  8:427291/146‑1 (MQ=255)
                                                                                                             gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                            >  6:149099/1‑140 (MQ=255)
                                                                                                             gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                            <  5:149099/140‑1 (MQ=255)
                                                                                                              ggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa                                     <  8:330553/146‑1 (MQ=255)
                                                                                                               gAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGa                                  <  7:322316/148‑1 (MQ=255)
                                                                                                                 aTTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTc                                                                                                     <  3:281850/79‑1 (MQ=255)
                                                                                                                 aTTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTc                                                                                                     >  4:281850/1‑79 (MQ=255)
                                                                                                                  ttGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCg                               <  3:97092/148‑1 (MQ=255)
                                                                                                                    gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATAc                                                                                                                       <  5:174468/58‑1 (MQ=255)
                                                                                                                    gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATAc                                                                                                                       >  6:174468/1‑58 (MQ=255)
                                                                                                                    gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCgg                              <  4:335417/147‑1 (MQ=255)
                                                                                                                     aTAAATCGGCATTGATTTAACTGCTAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGa                             <  1:393067/147‑1 (MQ=255)
                                                                                                                      tAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATg                           <  2:129808/148‑1 (MQ=255)
                                                                                                                       aaaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGa                          >  1:27312/1‑148 (MQ=255)
                                                                                                                        aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCAGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGAc                         >  5:420129/1‑148 (MQ=255)
                                                                                                                        aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGGTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGAc                         <  4:41672/148‑1 (MQ=255)
                                                                                                                                ttGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCgg                              <  6:22423/135‑1 (MQ=255)
                                                                                                                                ttGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCgg                              >  5:22423/1‑135 (MQ=255)
                                                                                                                                 tGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTATTCCGGCAGt                                                                                >  4:105556/1‑84 (MQ=255)
                                                                                                                                 tGATTTAACTGCAAATTGCCGGACAGATCTGCCGGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGt                                                                                <  3:105556/84‑1 (MQ=255)
                                                                                                                                    tttAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAt               >  6:382494/1‑146 (MQ=255)
                                                                                                                                    tttAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAt               <  5:382494/146‑1 (MQ=255)
                                                                                                                                      tAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtct            <  8:94455/147‑1 (MQ=255)
                                                                                                                                      tAACTGCAAATTGACGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa                                     >  4:31449/1‑122 (MQ=255)
                                                                                                                                      tAACTGCAAATTGACGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa                                     <  3:25848/122‑1 (MQ=255)
                                                                                                                                      tAACTGCAAATTGACGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa                                     >  4:25848/1‑122 (MQ=255)
                                                                                                                                      tAACTGCAAATTGACGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa                                     <  3:31449/122‑1 (MQ=255)
                                                                                                                                        aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcgg         >  5:61815/1‑148 (MQ=255)
                                                                                                                                        aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcgg         <  2:27312/148‑1 (MQ=255)
                                                                                                                                        aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcgg         <  1:28276/148‑1 (MQ=255)
                                                                                                                                        aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcg          <  7:133342/147‑1 (MQ=255)
                                                                                                                                          tGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGTTTAAGATAGGATGACGCCACAGCATTCTcggcc       >  1:91366/1‑148 (MQ=255)
                                                                                                                                            cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcc       <  6:145795/146‑1 (MQ=255)
                                                                                                                                                ttGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcccggcc  >  3:285126/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCC  >  NC_000913/2856280‑2856570

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