Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,466 T→G 16.6% Y384* (TAT→TAG clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4660TG16.6% 325.8 / 63.1 145Y384* (TAT→TAGclcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base T (57/63);  new base G (10/14);  total (68/77)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 147 total aligned reads displayed.

CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAT  >  NC_000913/1666368‑1666565
                                                                                                  |                                                                                                   
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                    >  1:3833069/1‑100 (MQ=255)
  tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                  <  1:1492588/100‑1 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              >  1:710693/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              >  1:1338809/1‑100 (MQ=255)
         ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                           <  1:3397676/100‑1 (MQ=255)
         ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                           <  1:2072500/100‑1 (MQ=255)
          gATTGACCGGGATGGCGGCAGTGTGGGAGGCGACCACGCACGCGCCGATTATGTCGGCGTTGATGTTATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                            <  2:927893/98‑1 (MQ=255)
          gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTccc                                                                                          <  1:2123536/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATAAGCTACTCCCCgg                                                                                       >  1:870836/1‑100 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:49475/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGACGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:1791397/100‑1 (MQ=38)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                       >  1:2440462/1‑99 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                      >  2:1250281/1‑100 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                      <  2:1675572/100‑1 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAATATCAGCTACTCCCCGGt                                                                                      >  1:2322261/1‑100 (MQ=255)
              cttccgatctgCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                      <  1:3737986/90‑1 (MQ=37)
                ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                    >  1:2818973/1‑100 (MQ=255)
                 cGGTATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATAGGTGAAATCACTGGTGTGTAGCTGATACTCCCCGGtcat                                                                                   <  2:2418281/100‑4 (MQ=255)
                 cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTa                                                                                   <  2:1720208/100‑1 (MQ=255)
                   ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                 >  2:1478528/1‑100 (MQ=255)
                   ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatc                                                                                 >  1:3927712/1‑99 (MQ=255)
                     aTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattga                                                                               <  1:2469197/100‑1 (MQ=255)
                      tGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                              <  2:424/100‑1 (MQ=255)
                       ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                             >  1:1881742/1‑100 (MQ=255)
                        gcgcCACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                             >  2:1210466/1‑99 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGACATGACCGGGGAGTATCAGCTACTCCCCGGTTCattgatt                                                                             >  2:413198/1‑99 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            >  1:1755804/1‑100 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            <  2:2974354/100‑1 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            <  2:2611915/100‑1 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtagatcggaag                                                                            >  2:3737947/1‑90 (MQ=38)
                         cGACACTGTTCGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                           <  2:2734964/100‑1 (MQ=255)
                           acacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                         >  1:861874/1‑100 (MQ=255)
                            cacTGTTGGCGGCGACCACGCACGCGTCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                        <  1:447197/100‑1 (MQ=255)
                            cacTGTTGGCGGCGACCACGCACGAGCAAATTATGTCGACGTAGATGATATGTGAAATTACCGGGGAGTATCAGTTACTCCCCAGTTTattgaatgac                                                                          >  2:3750064/1‑96 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                      <  1:397546/100‑1 (MQ=255)
                                   ggcggcGACCACGCACGCGCAGATTATGTCGACGTTGATGATATTTGAAATTACCGGGGAGTATCAGCTACTCCCTGGTTTATTGATTGCCTGCGTAAt                                                                  <  2:349708/99‑1 (MQ=255)
                                    gcggctACCACGCACGCGCCTATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTTCCTGCGTAATTg                                                                <  1:1139798/100‑1 (MQ=255)
                                    gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAAAATAAGATACTCCACGGTTTATTGATTGCCTGCGTAATTg                                                                >  1:2214535/1‑100 (MQ=255)
                                        cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                            <  1:3427785/100‑1 (MQ=255)
                                        cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                            <  1:2582791/100‑1 (MQ=255)
                                         gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCg                                                           <  2:3858676/100‑1 (MQ=255)
                                          aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                          <  2:2044913/100‑1 (MQ=255)
                                          aCCACGCACGCGCCGAATATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                          >  2:1644880/1‑100 (MQ=255)
                                           ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGTGAGTATCAGCTACTCCCCGGTTTATTGATTTCCTGCGTAATTGCGTc                                                            <  2:3278368/97‑1 (MQ=255)
                                           ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                         >  1:47299/1‑100 (MQ=255)
                                                cacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAttt                                                    >  1:417548/1‑100 (MQ=255)
                                                   gcgcCGTTTATGTCGACGATGATGAAATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTAGGTGATTTcgc                                                 >  2:1271382/1‑100 (MQ=255)
                                                    cgcCGATTATGTGGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                <  1:3624670/100‑1 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                <  1:181846/100‑1 (MQ=255)
                                                     gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                               <  1:2493708/100‑1 (MQ=255)
                                                     gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                               <  2:2828279/100‑1 (MQ=255)
                                                       cGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                             >  2:3847950/1‑100 (MQ=255)
                                                         aTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGt                                           <  1:2426299/100‑1 (MQ=255)
                                                          ttATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGtt                                          <  2:3036948/100‑1 (MQ=255)
                                                            aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                        >  1:146019/1‑100 (MQ=255)
                                                            aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                        <  2:2426502/100‑1 (MQ=255)
                                                            aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                        <  2:2121052/100‑1 (MQ=255)
                                                             tGTCGACGTTGATGATATTTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                       <  1:1250295/100‑1 (MQ=255)
                                                               tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAcc                                     <  2:1400923/100‑1 (MQ=255)
                                                                cGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                    >  1:2227329/1‑100 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                <  1:1381155/100‑1 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                <  2:3555053/100‑1 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                <  1:283735/100‑1 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                <  2:870825/100‑1 (MQ=255)
                                                                       gatgatATGTGAAATGACCGGGGAGTATCAGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAAAGTGActct                             >  2:195812/1‑100 (MQ=255)
                                                                        atgatATGTGAAATGACCGTGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                            <  1:2061134/100‑1 (MQ=255)
                                                                        atgatATGTGAAATGACCGGGGAGTATAAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                            <  2:3576951/100‑1 (MQ=255)
                                                                         tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                           >  2:1842569/1‑100 (MQ=255)
                                                                         tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGATAATTTCGCGAACGTTACACCGTGACtctat                           >  1:646557/1‑100 (MQ=255)
                                                                              tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                      >  2:58049/1‑100 (MQ=255)
                                                                               gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                     >  2:3150639/1‑100 (MQ=255)
                                                                               gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                     >  1:1251018/1‑100 (MQ=255)
                                                                               gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCAAACGTTACACCGTGACTCTATCTAcc                      >  1:893951/1‑99 (MQ=255)
                                                                               gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATAGCGTCGGTAATATCGCGGACGTTACACCGTGACTCTATCTACCg                     >  1:121773/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  2:968237/1‑100 (MQ=255)
                                                                                  aataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                  >  1:3350458/6‑100 (MQ=38)
                                                                                   aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                 <  2:1991486/100‑1 (MQ=255)
                                                                                   aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                 <  2:2991689/100‑1 (MQ=255)
                                                                                   aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                 <  2:1881727/100‑1 (MQ=255)
                                                                                    aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                <  2:2322239/100‑1 (MQ=255)
                                                                                    aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                <  2:1755789/100‑1 (MQ=255)
                                                                                      gACCGGGGAGTATCAGCTACTCCCTGGTTTATTGATTGTCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcac              <  1:3660833/100‑1 (MQ=255)
                                                                                        ncGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  2:107952/2‑100 (MQ=255)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:2969777/1‑100 (MQ=39)
                                                                                        ccGGGGAGTAGCTGATACTCCACAGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTAACACCGTGACTCTATCTACCGCCAGCagac            >  2:229814/1‑100 (MQ=255)
                                                                                         cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt           <  2:1210781/100‑1 (MQ=255)
                                                                                         cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt           >  1:2557051/1‑100 (MQ=255)
                                                                                         cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCTCCAGCACACt           >  2:2563163/1‑100 (MQ=255)
                                                                                         cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  2:1385326/1‑99 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          >  2:3389448/1‑100 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          >  1:979611/1‑100 (MQ=255)
                                                                                            ggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg        <  1:318822/100‑1 (MQ=255)
                                                                                              aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCTGACGTTACACCGTGAATCTATTTACCGACAGCACACTacg        >  1:789882/1‑98 (MQ=255)
                                                                                               gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag     <  2:506599/100‑1 (MQ=255)
                                                                                                tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcgc       <  1:195817/97‑1 (MQ=255)
                                                                                                 aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCCCACTGCgcagca   >  2:966601/1‑100 (MQ=255)
                                                                                                 aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca   >  1:722690/1‑100 (MQ=255)
                                                                                                 aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca   >  1:1019842/1‑100 (MQ=255)
                                                                                                  tCAGCTACTCCCCGGTTTATTGATTGCTTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt  <  1:1328354/100‑1 (MQ=255)
                                                                                                  tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt  <  1:2541388/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAT  >  NC_000913/1666368‑1666565

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