Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,716,030 T→G 9.9% intergenic (+20/+36) eamB → / ← grcA cysteine and O‑acetylserine exporter/autonomous glycyl radical cofactor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,716,0300TG 9.9% 222.5 / 15.6 71intergenic (+20/+36)eamB/grcAcysteine and O‑acetylserine exporter/autonomous glycyl radical cofactor
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base G (4/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATT  >  NC_000913/2715885‑2716175
                                                                                                                                                 |                                                                                                                                                 
tttGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGc                                                                                                                                                >  4:424321/1‑149 (MQ=255)
    gCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCctc                                                                                                                                            >  4:32027/1‑149 (MQ=255)
     cAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCctct                                                                                                                                           >  3:90315/1‑149 (MQ=255)
       aTGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCtt                                                                                                                                          <  2:828756/148‑1 (MQ=255)
       aTGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTc                                                                                                                                         <  4:365364/149‑1 (MQ=255)
       aTGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGACCTCtt                                                                                                                                          <  2:77299/148‑1 (MQ=255)
             gCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTa                                                                                                                                   <  4:495830/149‑1 (MQ=255)
                gggCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTc                                                                                                                                         <  2:237240/140‑1 (MQ=255)
                gggCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTc                                                                                                                                         >  1:237240/1‑140 (MQ=255)
                 ggCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGACCTCTTCCGCtttttt                                                                                                                                <  1:171660/148‑5 (MQ=255)
                          ggggCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACt                                                                                                                      >  4:772858/1‑149 (MQ=255)
                             gCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTaa                                                                                                                   <  4:90315/149‑1 (MQ=255)
                             gCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTa                                                                                                                    <  1:24694/148‑1 (MQ=255)
                                tCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTc                                                                                                                                         <  4:541792/124‑1 (MQ=255)
                                tCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTc                                                                                                                                         >  3:541792/1‑124 (MQ=255)
                                tCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGc                                                                                                                 >  1:723673/1‑148 (MQ=255)
                                 cTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGcc                                                                                                                >  2:168407/1‑148 (MQ=255)
                                  tGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCtg                                                                                                              >  4:23421/1‑148 (MQ=255)
                                  tGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCt                                                                                                               >  1:771075/1‑148 (MQ=255)
                                   gTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCtt                                                                                                              <  1:271837/148‑1 (MQ=255)
                                    tttCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTAc                                                                                                            <  3:314039/149‑1 (MQ=255)
                                       cAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                               >  2:769243/1‑111 (MQ=255)
                                       cAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                               <  1:769243/111‑1 (MQ=255)
                                       cAGCGATTGTTTCGCCAGTATGGTCACCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTGCGCATTTTCTATTAACGAAAAAAAGCGGAAGAGATCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAgg                                                                                                         >  4:760615/1‑149 (MQ=255)
                                          cGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCt                                                                                                       >  2:195043/1‑148 (MQ=255)
                                                tttCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTAc                                                                                                                       <  1:882436/126‑1 (MQ=255)
                                                tttCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTAc                                                                                                                       >  2:882436/1‑126 (MQ=255)
                                                   cGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAg                                                                                             <  3:424321/149‑1 (MQ=255)
                                                    gCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAg                                                                                             <  2:415526/148‑1 (MQ=255)
                                                      cAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGt                                                                                           >  2:120817/1‑148 (MQ=255)
                                                       aGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTa                                                                                          <  1:195043/148‑1 (MQ=255)
                                                                 ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAAc                                                                               >  4:379958/1‑149 (MQ=255)
                                                                         atatCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCgcgc                                                                        <  1:120817/148‑1 (MQ=255)
                                                                            tCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCgctgct                                                                    >  4:194805/1‑149 (MQ=255)
                                                                            tCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCgctgc                                                                     <  2:205565/148‑1 (MQ=255)
                                                                             cGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTGTAAGGCtta                                                                                                     >  1:849862/1‑114 (MQ=255)
                                                                               tGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTa                                                                                                             <  2:886760/105‑1 (MQ=255)
                                                                               tGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTa                                                                                                             >  1:886760/1‑105 (MQ=255)
                                                                                gCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGtt                                                                 <  2:397238/148‑1 (MQ=255)
                                                                                    gCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGa                                                            <  3:379958/149‑1 (MQ=255)
                                                                                       cTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGt                                                          <  2:723673/148‑1 (MQ=255)
                                                                                         gTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCag                                                       >  3:522489/1‑149 (MQ=255)
                                                                                           tGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagag                                                     <  3:328166/149‑1 (MQ=255)
                                                                                            gCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagag                                                     <  1:328882/148‑1 (MQ=255)
                                                                                                  cTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGTAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCGTTACAGGCTTTCAGTAAAGGTACGAGCTATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaaa                                              <  3:760615/149‑1 (MQ=255)
                                                                                                  cTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaa                                               <  1:168407/148‑1 (MQ=255)
                                                                                                   tATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAg                                             <  3:772858/149‑1 (MQ=255)
                                                                                                      tGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGa                                          >  3:702320/1‑149 (MQ=255)
                                                                                                      tGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGa                                          <  4:738577/149‑1 (MQ=255)
                                                                                                        cgcgGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAAc                                        <  4:522489/149‑1 (MQ=255)
                                                                                                         gcgGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAAc                                        >  1:721078/1‑148 (MQ=255)
                                                                                                                         aTTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATacac                       >  3:214215/1‑149 (MQ=255)
                                                                                                                            aaCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGa                    <  3:468890/149‑1 (MQ=255)
                                                                                                                                aaaaaaaGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGa                            >  3:258400/1‑137 (MQ=255)
                                                                                                                                aaaaaaaGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGa                            <  4:258400/137‑1 (MQ=255)
                                                                                                                                  aaaaaGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTc               >  2:626904/1‑148 (MQ=255)
                                                                                                                                     aaGCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagaga                                                    >  3:748766/1‑108 (MQ=255)
                                                                                                                                     aaGCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagaga                                                    <  4:748766/108‑1 (MQ=255)
                                                                                                                                       gCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTg          >  2:162512/1‑148 (MQ=255)
                                                                                                                                         ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCgg       >  4:455651/1‑149 (MQ=255)
                                                                                                                                          gAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaaa                                              >  2:821157/1‑109 (MQ=255)
                                                                                                                                          gAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaaa                                              <  1:821157/109‑1 (MQ=255)
                                                                                                                                          gAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCgg       >  2:178722/1‑148 (MQ=255)
                                                                                                                                          gAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGa      >  3:407675/1‑149 (MQ=255)
                                                                                                                                           aaGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagaga                                                    <  4:218047/102‑1 (MQ=255)
                                                                                                                                           aaGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagaga                                                    >  3:218047/1‑102 (MQ=255)
                                                                                                                                            agagGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGat     >  2:895056/1‑148 (MQ=255)
                                                                                                                                            agagGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGat     >  1:188910/1‑148 (MQ=255)
                                                                                                                                             gagGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGatat   <  3:32027/149‑1 (MQ=255)
                                                                                                                                             gagGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGata    >  1:104158/1‑148 (MQ=255)
                                                                                                                                              agGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATAtt  <  4:214215/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATT  >  NC_000913/2715885‑2716175

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