Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,615,859 T→G 100% T403P (ACC→CCC)  yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,615,8590TG98.8% 537.7 / ‑3.1 157T403P (ACC→CCC) yjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (81/73);  total (82/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 160 total aligned reads displayed.

CCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTT  >  NC_000913/4615763‑4615954
                                                                                                |                                                                                               
ccTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGAt                                                                                              <  2:1618366/100‑1 (MQ=255)
ccTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGAt                                                                                              >  2:237847/1‑100 (MQ=255)
ccTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGAt                                                                                              >  2:2374670/1‑100 (MQ=255)
  tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                            >  2:3265525/1‑100 (MQ=255)
  tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                            <  1:951330/100‑1 (MQ=255)
    cGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGa                                                                                          <  1:1998767/100‑1 (MQ=255)
     aggatACGCTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGa                                                                                          >  1:362006/5‑99 (MQ=255)
      ggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGggatcggatc                                                                                        >  1:2463420/1‑100 (MQ=255)
        ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAATCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                      <  2:4241795/100‑1 (MQ=255)
        ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                      <  1:3697212/100‑1 (MQ=255)
          aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    <  2:389975/100‑1 (MQ=255)
          aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    >  2:2634197/1‑100 (MQ=255)
           cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                   >  2:2694378/1‑100 (MQ=255)
            gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    <  1:3307241/98‑1 (MQ=255)
            gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTc                                                                                  >  1:3175317/1‑100 (MQ=255)
            gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGAGGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                   >  2:1211166/1‑99 (MQ=255)
              ttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCAt                                                                                <  2:51022/100‑1 (MQ=255)
               tGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATc                                                                               >  1:239672/1‑100 (MQ=255)
               tGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATc                                                                               >  2:20456/1‑100 (MQ=255)
                gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                              >  1:760315/1‑100 (MQ=255)
                gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGAGGCGACAGTTTGCAGATGGGGGATCGGATCTGGcaca                                                                                 >  1:243487/1‑97 (MQ=255)
                    gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                          <  1:1614478/100‑1 (MQ=255)
                       tcgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACg                                                                       <  2:1174741/100‑1 (MQ=255)
                        cgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACgg                                                                      >  1:527848/1‑100 (MQ=255)
                        cgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACg                                                                       >  1:2169050/1‑99 (MQ=255)
                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACgac                                                                     >  1:1437731/1‑98 (MQ=255)
                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     <  1:4268884/100‑1 (MQ=255)
                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     <  2:830503/100‑1 (MQ=255)
                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     >  1:579223/1‑100 (MQ=255)
                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     >  2:317618/1‑100 (MQ=255)
                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     <  2:3798383/100‑1 (MQ=255)
                          tcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCa                                                                    >  2:2571094/1‑100 (MQ=255)
                          tcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCa                                                                    <  1:1664775/100‑1 (MQ=255)
                             aGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAc                                                                 >  2:3666797/1‑100 (MQ=255)
                              gCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcg                                                                <  2:253913/100‑1 (MQ=255)
                                 tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                             >  1:738834/1‑100 (MQ=255)
                                  cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                            <  2:3707941/100‑1 (MQ=255)
                                     aaTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGccc                                                         <  1:3059288/100‑1 (MQ=255)
                                        gTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGc                                                      >  1:2206274/1‑100 (MQ=255)
                                         tCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCg                                                     >  2:1730209/1‑100 (MQ=255)
                                          cGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgt                                                    >  2:1824068/1‑100 (MQ=255)
                                          cGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgt                                                    >  2:72237/1‑100 (MQ=255)
                                            cTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtgg                                                  <  2:163427/100‑1 (MQ=255)
                                             tGATGCCGGAATAATAATAAGCATGATGCGGGGCGACGGTTTGCAGATGGGTGATCGGATCTGGCTCATCACCATACGGCAAACGTGCGCCCGGCgtgg                                                  >  1:4171093/1‑99 (MQ=255)
                                              gATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggtg                                                >  2:2814733/1‑100 (MQ=255)
                                               aTGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                               >  2:454411/1‑100 (MQ=255)
                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              <  1:3966298/100‑1 (MQ=255)
                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              >  2:3276170/1‑100 (MQ=255)
                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              >  2:1012084/1‑100 (MQ=255)
                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              >  1:2508883/1‑100 (MQ=255)
                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGAGGGGGGATCGGATCGGGCTCAGCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              <  1:969308/100‑1 (MQ=255)
                                                 gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                             >  1:1750902/1‑100 (MQ=255)
                                                  ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                            <  1:1193998/100‑1 (MQ=255)
                                                   cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                           >  1:1897666/1‑100 (MQ=255)
                                                    ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGata                                          >  2:2230640/1‑100 (MQ=255)
                                                       ataataataAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCg                                       <  1:4233269/100‑1 (MQ=255)
                                                        taataataaGCATGATGCGGGGCGACAGTTTGCAGATGTGCGATAGGATCTGGTTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                      <  1:88010/100‑1 (MQ=255)
                                                        taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                      <  2:613232/100‑1 (MQ=255)
                                                         aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGa                                     >  2:89786/1‑100 (MQ=255)
                                                          ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                    <  2:760317/100‑1 (MQ=255)
                                                          ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                    >  1:3060796/1‑100 (MQ=255)
                                                           taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                   >  1:1049006/1‑100 (MQ=255)
                                                           taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                   <  1:2619363/100‑1 (MQ=255)
                                                            aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                  <  2:716908/100‑1 (MQ=255)
                                                            aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                  >  1:3085062/1‑100 (MQ=255)
                                                            aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                  <  1:3307253/100‑1 (MQ=255)
                                                               aaGCAGGATGCGGGGCGACAGTTTGCAGATGGGGGCTCCGATCTGGCTCATCGCCCTACGGCAAACGCGCGCCCGGCGTTGTGCCGATATCGGAACTg                                 <  1:1519621/98‑1 (MQ=255)
                                                                aGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATc                              <  1:1352300/100‑1 (MQ=255)
                                                                 gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                             >  1:1841014/1‑100 (MQ=255)
                                                                 gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                             >  2:2711145/1‑100 (MQ=255)
                                                                 gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                             >  2:2318656/1‑100 (MQ=255)
                                                                 gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                             >  2:4276228/1‑100 (MQ=255)
                                                                 gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                             >  1:2619816/1‑100 (MQ=255)
                                                                  cATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATccc                            <  2:523644/100‑1 (MQ=255)
                                                                    tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                          >  2:3101308/1‑100 (MQ=255)
                                                                      atgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGAc                        <  1:1426088/100‑1 (MQ=255)
                                                                      atgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGAc                        >  2:2928014/1‑100 (MQ=255)
                                                                        gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                      <  1:1054508/100‑1 (MQ=255)
                                                                         cGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgt                     >  1:4000020/1‑100 (MQ=255)
                                                                          ggggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtg                    <  1:89786/100‑1 (MQ=255)
                                                                           gggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGc                   <  2:4195809/100‑1 (MQ=255)
                                                                           gggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGAGATCGGAACTGATCCCCGAccatac                   >  1:3899280/1‑95 (MQ=255)
                                                                            ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                  <  1:411554/100‑1 (MQ=255)
                                                                            ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                  <  2:3989435/100‑1 (MQ=255)
                                                                              cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                <  2:1610337/100‑1 (MQ=255)
                                                                                aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa              <  2:2173408/100‑1 (MQ=255)
                                                                                aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa              >  2:236532/1‑100 (MQ=255)
                                                                                aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa              <  1:310846/100‑1 (MQ=255)
                                                                                aCAGGTTGCAGATCGGGGATCGGCGCTGGTGCATCGCCAGACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa              <  1:254463/100‑1 (MQ=255)
                                                                                 cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa             <  1:4095919/100‑1 (MQ=255)
                                                                                 cAGTTTGCAGATGGGGGATCGGATCTGGCGCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa              <  1:1170753/99‑1 (MQ=255)
                                                                                  aGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAAc            >  2:3340577/1‑100 (MQ=255)
                                                                                  aGTTTGCAGATGGGGGAGCGGATCTGGCTCAGCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAAc            <  1:1533516/100‑1 (MQ=255)
                                                                                   gTTTGCAGATGGGGGATCGGATCTGGCTTATCGCTATACGGCAACCGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa           <  1:520108/100‑1 (MQ=255)
                                                                                   gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa           >  2:3475798/1‑100 (MQ=255)
                                                                                   gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa           <  2:2124953/100‑1 (MQ=255)
                                                                                      tGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg        <  2:3995439/100‑1 (MQ=255)
                                                                                        cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCTGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGtcg      >  2:379502/1‑100 (MQ=255)
                                                                                        cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg      <  2:4214246/100‑1 (MQ=255)
                                                                                            tGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGtt  >  1:2665070/1‑100 (MQ=255)
                                                                                                |                                                                                               
CCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTT  >  NC_000913/4615763‑4615954

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