Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,170 +C intergenic (+172/+279) hemN → / ← glnG O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,1701.C100.0% 970.5 / NA 351intergenic (+172/+279)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (188/163);  total (188/163)

Only 100 of 382 total aligned reads displayed.

TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACA  >  NC_000913/4051077‑4051256
                                                                                              |                                                                                      
tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                                                         >  1:403145/1‑93 (MQ=255)
tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                                                         >  2:5931394/1‑93 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                                        >  2:1426575/1‑92 (MQ=255)
  ggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgct                                                                                       >  1:3644975/1‑91 (MQ=255)
        aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                 >  1:4896748/1‑93 (MQ=25)
        aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                 <  1:587731/93‑1 (MQ=21)
          aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                                               >  1:4773600/1‑93 (MQ=25)
           aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTc                                                                              >  2:6115004/1‑93 (MQ=21)
             cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                            >  2:40678/1‑93 (MQ=21)
             cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                            >  1:2937837/1‑93 (MQ=21)
              tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                                           >  1:2551322/1‑93 (MQ=21)
               gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                                          >  2:1608396/1‑93 (MQ=21)
                ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                                         <  1:7638748/93‑1 (MQ=21)
                 tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                        <  2:1044368/93‑1 (MQ=21)
                 tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                        <  2:470343/93‑1 (MQ=21)
                 tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                        <  2:462644/93‑1 (MQ=18)
                  gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                                       <  2:6555577/93‑1 (MQ=18)
                  gAAATTATTCCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                                       >  2:5238448/1‑93 (MQ=18)
                   aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                                                                      <  1:7170501/93‑1 (MQ=18)
                    aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                     <  2:1196254/93‑1 (MQ=18)
                     attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                                    <  1:5130382/93‑1 (MQ=18)
                     attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                                    <  1:4364784/93‑1 (MQ=18)
                      ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCTGCGCTGGTCGCGTCTTATCAGGCCt                                                                   >  2:6518690/1‑93 (MQ=18)
                      ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                   <  1:924041/93‑1 (MQ=18)
                      ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                   >  2:7239589/1‑93 (MQ=18)
                       tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                  >  1:2132772/1‑93 (MQ=17)
                       tattGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCATGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                  >  1:5315838/1‑93 (MQ=17)
                        attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                                 >  1:7418423/1‑93 (MQ=18)
                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                               <  2:5769260/93‑1 (MQ=18)
                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                               <  1:7072314/93‑1 (MQ=18)
                           gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                                              >  2:3383981/1‑93 (MQ=17)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            >  2:2931412/1‑93 (MQ=17)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            >  2:7902399/1‑93 (MQ=17)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            <  1:3343706/93‑1 (MQ=17)
                             ttttGCAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            <  1:7743821/93‑1 (MQ=17)
                              tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                                           >  1:6537498/1‑93 (MQ=17)
                               ttGTAGTCCGGATAAGTCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                                          <  1:6365350/93‑1 (MQ=17)
                                  tAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                       >  2:486187/1‑93 (MQ=18)
                                    ggCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                                                     <  1:3852731/93‑1 (MQ=18)
                                     gCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                                    >  2:3027617/1‑93 (MQ=18)
                                       cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                                  >  1:290205/1‑93 (MQ=18)
                                       cGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                                  >  1:3456803/1‑93 (MQ=18)
                                        ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                                 >  1:5767232/1‑93 (MQ=18)
                                        ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                                 <  2:5105281/93‑1 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                >  2:2354109/1‑93 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGCCGCTCCTTATCAGGCCTCCGCTTTTCTCCGACTCCg                                                >  1:4554583/1‑93 (MQ=9)
                                         gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCCGCGCTGGTCGCGTCTTATCAGGCCTACCACTTTCTCCGCATCCg                                                >  2:4029761/1‑93 (MQ=18)
                                          aTAAGGCGTTTACGCCGCATCAGG‑CATCTTAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                                               <  1:6206312/93‑1 (MQ=18)
                                           tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                              >  1:320390/1‑93 (MQ=18)
                                           tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                              >  2:637286/1‑93 (MQ=21)
                                            aaGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCa                                             >  1:396199/1‑93 (MQ=21)
                                              ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                                           >  2:2440664/1‑93 (MQ=21)
                                              ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                                           >  2:6986175/1‑93 (MQ=21)
                                               gCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                                          <  2:1377290/93‑1 (MQ=21)
                                                 gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                        >  1:5206303/1‑93 (MQ=25)
                                                     ccgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                                    <  2:5661747/92‑1 (MQ=25)
                                                     aCGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                                    <  1:3560743/93‑1 (MQ=25)
                                                     aCGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                                    >  2:3713275/1‑93 (MQ=25)
                                                      cgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa                                   <  2:3470275/93‑1 (MQ=25)
                                                      cgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa                                   <  2:6621085/93‑1 (MQ=25)
                                                      cgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa                                   <  1:6730841/93‑1 (MQ=25)
                                                       gccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg                                  >  2:827923/1‑93 (MQ=25)
                                                        ccgcATCAGG‑CATCTGAGCATGATGCCTTATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                 <  1:1709941/93‑1 (MQ=25)
                                                        ccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                 <  2:7132502/93‑1 (MQ=25)
                                                         cgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCt                                <  1:1076059/93‑1 (MQ=25)
                                                         cgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCt                                <  1:2394228/93‑1 (MQ=25)
                                                          gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                               >  1:1029816/1‑93 (MQ=37)
                                                           cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                              <  1:2765066/93‑1 (MQ=37)
                                                           cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                              <  2:2215351/93‑1 (MQ=37)
                                                           cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTACGCATCCGGCATAAACAACGAGCTTc                              <  2:7392046/93‑1 (MQ=25)
                                                            aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                             >  2:3519736/1‑93 (MQ=37)
                                                             tCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                            >  2:1948652/1‑93 (MQ=37)
                                                             tCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                            >  1:1992777/1‑93 (MQ=37)
                                                                 ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                         >  2:421273/3‑93 (MQ=255)
                                                                   cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa                       <  1:6956384/93‑1 (MQ=255)
                                                                    aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                      <  1:1468115/93‑1 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                     >  2:1876072/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                     >  1:7389843/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                     >  2:2688563/1‑93 (MQ=255)
                                                                      cTGAGCATGATGCTTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAg                    <  2:7806322/93‑1 (MQ=255)
                                                                      cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAg                    >  2:3846448/1‑93 (MQ=255)
                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                   >  1:6639440/1‑93 (MQ=255)
                                                                        gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                  <  1:4439657/93‑1 (MQ=255)
                                                                        gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                  >  2:3174637/1‑93 (MQ=255)
                                                                          gaatgatcCTTGATGCGGCGCTGGTCGCGCCTTATCAGGACTACGACTTTCTCCGCATCCGGCATAGACAACGAGATTCAGCTAAACAGcaca                >  2:5395675/3‑93 (MQ=255)
                                                                            atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt              <  1:5505025/93‑1 (MQ=255)
                                                                            atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt              <  1:1877153/93‑1 (MQ=255)
                                                                            atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGCCATAAAAAACGAGCTTCAGCTAAACAGCCCaac              >  2:5888815/1‑92 (MQ=37)
                                                                             tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc             >  1:3105830/1‑93 (MQ=255)
                                                                             tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc             >  2:5828223/1‑93 (MQ=255)
                                                                              gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCa            <  1:7879342/93‑1 (MQ=255)
                                                                               atgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAt           >  1:5560640/1‑93 (MQ=255)
                                                                               atgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAt           >  2:4463001/1‑93 (MQ=255)
                                                                                tgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt          <  2:2522696/93‑1 (MQ=255)
                                                                                tgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt          <  1:5645188/93‑1 (MQ=255)
                                                                                tgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt          <  2:6826815/93‑1 (MQ=255)
                                                                                  ccTGATGCGGCGCTTGGCGGGTTTTTTCCGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc        <  1:6003830/93‑1 (MQ=37)
                                                                                  ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc        <  2:7418422/93‑1 (MQ=255)
                                                                                    tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCgg      >  2:7634406/1‑93 (MQ=255)
                                                                                        gcggcgcTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcaca  >  1:1258611/1‑93 (MQ=255)
                                                                                              |                                                                                      
TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACA  >  NC_000913/4051077‑4051256

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