Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 629,933 T→C W273R (TGG→CGG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913629,9330TC98.9% 270.4 / ‑4.7 91W273R (TGG→CGG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (43/47);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

Only 100 of 112 total aligned reads displayed.

CGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGC  >  NC_000913/629844‑630022
                                                                                         |                                                                                          
cGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGt                                                                                         <  1:314447/93‑1 (MQ=255)
cGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGt                                                                                         <  2:1834115/93‑1 (MQ=255)
  gccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTa                                                                                       <  2:1994350/92‑1 (MQ=255)
   ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAc                                                                                      >  1:2977090/1‑93 (MQ=255)
   ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAc                                                                                      >  2:3673782/1‑93 (MQ=255)
    cTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAct                                                                                     <  1:4510996/93‑1 (MQ=255)
     tGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTTctg                                                                                    >  2:1819267/1‑93 (MQ=255)
       ggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgct                                                                                  <  1:4036475/93‑1 (MQ=255)
        gCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgctg                                                                                 <  2:730098/93‑1 (MQ=255)
         cACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTgg                                                                                >  2:181383/1‑93 (MQ=255)
         cACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTgg                                                                                >  2:181384/1‑93 (MQ=255)
          aCCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGc                                                                               >  2:573972/1‑93 (MQ=255)
           ccGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTActgcttggc                                                                              >  2:171698/1‑89 (MQ=255)
             gTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGcccc                                                                            <  2:4234474/93‑1 (MQ=255)
             gTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGcccc                                                                            <  2:3091317/93‑1 (MQ=255)
                cTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCcgcg                                                                         <  2:2325199/93‑1 (MQ=255)
                 tttGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGt                                                                       <  2:3500704/93‑1 (MQ=255)
                  ttGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGTCCCCCGTg                                                                      >  2:495247/1‑93 (MQ=255)
                   tGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGGTACTGCTGGCCCCGCGTGa                                                                     >  1:1403443/1‑93 (MQ=255)
                    gACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGCCCCGCGTGAc                                                                    >  2:1279661/1‑93 (MQ=255)
                      cTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGGTGCGGGTGCGGTTACTGCTGGCCCCGGGGGACtc                                                                  >  1:1462094/1‑92 (MQ=255)
                      cTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGCTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTgcata                                                                  >  1:4578423/1‑89 (MQ=255)
                       tttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGGCCCGCGTGACTAc                                                                 >  2:1569364/1‑93 (MQ=255)
                       tttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGCGGTTACTGCTGGCCCCGCGGGACTAc                                                                 >  1:621697/1‑93 (MQ=255)
                        ttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGCGGGTACTCGTGGCCCCGCGGGACTAcc                                                                >  2:962242/1‑93 (MQ=255)
                         tACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACct                                                               >  1:3317707/1‑93 (MQ=255)
                          aCCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctc                                                              <  2:2768856/93‑1 (MQ=255)
                            cGGCGTGCGGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctc                                                            <  1:830037/93‑1 (MQ=255)
                            cGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctc                                                            <  1:1835465/93‑1 (MQ=255)
                            cGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctc                                                            <  1:155058/93‑1 (MQ=255)
                                gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAcc                                                        <  2:1174197/93‑1 (MQ=255)
                                gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAcc                                                        >  2:3011758/1‑93 (MQ=255)
                                gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGCGGTTACTGGTGGCCCCGCGGGACTACCTGTCTAcc                                                        >  2:1115679/1‑93 (MQ=255)
                                 tGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGCCCCGCGTGACTAACTCTCTACCt                                                       >  2:2783102/1‑93 (MQ=255)
                                   cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                                     >  1:283354/1‑93 (MQ=255)
                                   cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                                     <  1:3446198/93‑1 (MQ=255)
                                     gTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGGCCCGCGTGACTACCTCTCTACcctcct                                                   >  2:4199173/1‑93 (MQ=255)
                                     gTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                                   >  1:3284866/1‑93 (MQ=255)
                                      ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGGCCCGCGTGACTACCTCTCTAACTTCCTg                                                  >  2:2492577/1‑93 (MQ=255)
                                      ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTg                                                  <  2:4269887/93‑1 (MQ=255)
                                      ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTg                                                  <  1:4012266/93‑1 (MQ=255)
                                       tGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGa                                                 <  2:1923366/93‑1 (MQ=255)
                                        gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCCTGGCTACCTCTCTCCCTTCCTGaa                                                >  1:4022762/1‑93 (MQ=255)
                                         aCCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                               >  1:462133/1‑93 (MQ=255)
                                         aCCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGACCCGCGTGACTACCTCTCTACCTTCCTGaaa                                               >  1:557999/1‑93 (MQ=255)
                                          ccTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaaa                                              <  1:3684457/93‑1 (MQ=255)
                                           cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAt                                             <  2:3263456/93‑1 (MQ=255)
                                            tGGATGGTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            <  1:3650379/93‑1 (MQ=255)
                                            tGGATGCTGGTGTGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            >  2:3037763/1‑93 (MQ=255)
                                            tGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTATCTTCCTGAAAATc                                            >  2:125174/1‑93 (MQ=255)
                                            tGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaaaac                                            >  1:1914332/1‑91 (MQ=255)
                                            tGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            >  1:182458/1‑93 (MQ=255)
                                             ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCg                                           >  1:2247288/1‑93 (MQ=255)
                                              gATGCTGGTTGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTCCCTCTCTACCTTCCTGAAAATCgg                                          <  2:3818012/93‑1 (MQ=255)
                                              gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCgg                                          >  2:3647543/1‑93 (MQ=255)
                                              gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCgg                                          <  2:4712728/93‑1 (MQ=255)
                                              gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGCGGGTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCgg                                          >  2:4066352/1‑93 (MQ=255)
                                               aTGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCggg                                         <  2:3118160/93‑1 (MQ=255)
                                                tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGa                                        <  2:3680577/93‑1 (MQ=255)
                                                tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGa                                        <  2:3680576/93‑1 (MQ=255)
                                                tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGa                                        >  1:2733190/1‑93 (MQ=255)
                                                tGCTGGTGGGTTACGGTTCTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTGCCTTCCTGAAAATCGGGa                                        <  2:1955532/93‑1 (MQ=255)
                                                 gCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGAc                                       <  2:204278/93‑1 (MQ=255)
                                                 gCTGGTGCGGTACGGTTTTGTGTCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAGCGGGAc                                       <  2:499583/93‑1 (MQ=255)
                                                   tggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTa                                     >  1:3116941/1‑93 (MQ=255)
                                                   tggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTa                                     >  1:3364487/1‑93 (MQ=255)
                                                    ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGCGGACTCCTCTCTCCCTTCCTGCAAATCAGGACtct                                    >  1:2602509/1‑91 (MQ=255)
                                                    ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAt                                    <  1:549129/93‑1 (MQ=255)
                                                     gtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtc                                   <  2:2891494/93‑1 (MQ=255)
                                                      tggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcg                                  >  1:2939111/1‑93 (MQ=255)
                                                         gTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgtcg                               >  1:2787022/1‑93 (MQ=255)
                                                           tACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGt                             >  2:2666590/1‑93 (MQ=255)
                                                           tACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGt                             >  2:556309/1‑93 (MQ=255)
                                                            aCGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTc                            <  2:4373689/93‑1 (MQ=255)
                                                             cGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                           <  1:3158395/93‑1 (MQ=255)
                                                             cGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                           <  2:2188179/93‑1 (MQ=255)
                                                                ttttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGc                        <  1:4062698/93‑1 (MQ=255)
                                                                ttttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGc                        <  2:3868085/93‑1 (MQ=255)
                                                                ttttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGc                        <  2:384089/93‑1 (MQ=255)
                                                                ttttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAACTTCCTGAAAATCGGGACTATCGTCGGTCTGGc                        >  1:4594133/1‑93 (MQ=255)
                                                                  ttatgGCGACGATGCTGCCGGTGCGGTTACTACAGGCCAAGCGTGACTACCTCTCTACCTTCCTGACAATCAGGACTATCGTCGGTCTGGCgg                      <  2:4222727/90‑1 (MQ=255)
                                                                  ttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCgg                      >  1:1373769/1‑93 (MQ=255)
                                                                  ttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCgg                      <  2:3485173/93‑1 (MQ=255)
                                                                   tgtgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGt                     <  1:1032451/93‑1 (MQ=255)
                                                                   tggggcggcGGTGCTCCCGGTGCGGTTACTGCGGGCCCCGCGTGACTACCTCTCTCCCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGt                     <  2:1091306/90‑1 (MQ=38)
                                                                    gtgGCGGCGGGGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTa                    >  1:4286058/1‑93 (MQ=255)
                                                                     tgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAg                   <  2:2365637/93‑1 (MQ=255)
                                                                     tgGCGGCGGCGCGGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCCCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAg                   <  2:3997214/93‑1 (MQ=255)
                                                                       gcggcgGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGc                 <  2:1523208/93‑1 (MQ=255)
                                                                            ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCAGGACTATCGTCGGTCTGGCGGTAGGCAtttt            >  2:2732891/1‑93 (MQ=255)
                                                                             gTGCTGCCGGTGCGGTTACTGCTTGGCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTg           >  2:1523778/1‑93 (MQ=255)
                                                                             gTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTg           >  2:1523777/1‑93 (MQ=255)
                                                                              tgctgcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGa          >  1:1509934/1‑93 (MQ=255)
                                                                               gctgcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAt         <  1:3104562/93‑1 (MQ=255)
                                                                                  gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      >  2:990436/1‑93 (MQ=255)
                                                                                  gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      >  2:1735140/1‑93 (MQ=255)
                                                                                  gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      <  2:3168148/93‑1 (MQ=255)
                                                                                     ggTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGAATTTTGATTATgcg   >  1:838801/1‑93 (MQ=255)
                                                                                      gtgcGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgcgc  >  1:3634480/1‑93 (MQ=255)
                                                                                      gtgcGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgcgc  >  1:2814062/1‑93 (MQ=255)
                                                                                         |                                                                                          
CGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGC  >  NC_000913/629844‑630022

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