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,9330TC99.4% 491.7 / NA 158W273R (TGG→CGG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (65/92);  minor base G (1/0);  total (66/92)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

CCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGC  >  NC_000913/629843‑630022
                                                                                          |                                                                                         
  gACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCgggt                                                                                       >  2:3515651/1‑91 (MQ=255)
   aCCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTa                                                                                      >  2:4846488/1‑93 (MQ=255)
    ccTGGGCACCGTACTTTGACTTTACCGGCGTGCGGTTGACCTGGATGCTGGTGGGTTTCGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTAc                                                                                     >  2:1760635/1‑93 (MQ=255)
    ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAc                                                                                     <  1:5968596/93‑1 (MQ=255)
    ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAc                                                                                     <  2:6600487/93‑1 (MQ=255)
       gggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgc                                                                                  <  2:3637152/93‑1 (MQ=255)
        ggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgct                                                                                 <  2:4574836/93‑1 (MQ=255)
           aCCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGc                                                                              >  2:138936/1‑93 (MQ=255)
           aCCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGc                                                                              <  1:3478048/93‑1 (MQ=255)
            ccGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTggac                                                                             >  2:2537055/1‑91 (MQ=255)
            ccGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGcc                                                                             >  2:5568402/1‑93 (MQ=255)
             cGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGGCGGTGCGGGTACTGCTGGccc                                                                            >  1:623892/1‑93 (MQ=255)
             cGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGgcc                                                                            >  1:453339/1‑93 (MQ=255)
              gTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGGTGGCGGGGCGGTTACTGCTGGcccc                                                                           >  2:64025/1‑93 (MQ=255)
              gTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGcccc                                                                           <  1:4862883/93‑1 (MQ=255)
              gTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGcccc                                                                           <  1:5504967/93‑1 (MQ=255)
              gTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGcccc                                                                           >  1:6594516/1‑93 (MQ=255)
              gTACTTTGACTTTACCGGCGTGCAGGTGACCTGGATGCTGGTGGGTTACGGGTTTGTTGCGGCGGGGATGCCGGTGCGGTGACTGCTGGcccc                                                                           >  2:1473222/1‑93 (MQ=255)
               tACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCcg                                                                          <  1:4908702/93‑1 (MQ=255)
               tACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAATTCTGGGCCcg                                                                          >  1:1893432/1‑93 (MQ=255)
                aCTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCcgc                                                                         <  1:695258/93‑1 (MQ=255)
                  tttGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGGCCCGCGt                                                                       >  1:6303640/1‑93 (MQ=255)
                  tttGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGt                                                                       >  1:1036961/1‑93 (MQ=255)
                    tGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGa                                                                     >  2:1574776/1‑93 (MQ=255)
                     taCTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGAc                                                                    <  2:728742/92‑1 (MQ=255)
                     gACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGAc                                                                    <  2:5090524/93‑1 (MQ=255)
                      aCTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACt                                                                   >  2:2265888/1‑93 (MQ=255)
                       cTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTa                                                                  <  1:3011254/93‑1 (MQ=255)
                        tttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTAc                                                                 >  2:3138542/1‑93 (MQ=255)
                        tttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTAc                                                                 <  2:831581/93‑1 (MQ=255)
                        tttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGCGGTTACTGCTGGCCCCGCGTGACTAc                                                                 >  2:1548748/1‑93 (MQ=255)
                          tACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACct                                                               <  2:204248/93‑1 (MQ=255)
                           aCCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctc                                                              <  1:2790511/93‑1 (MQ=255)
                            ccGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctct                                                             <  2:4774917/93‑1 (MQ=255)
                            ccGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGGCCCGCGTGACTAActct                                                             >  1:3778383/1‑93 (MQ=255)
                            ccGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGATTTGTGGCGGGGGCGGCGCCGGGGGGGTTACTGCTGGGCCCGCGAGACTACctct                                                             >  1:2851574/1‑93 (MQ=255)
                             cGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctc                                                            >  1:3826689/1‑93 (MQ=255)
                              ggCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctct                                                           <  2:1553706/93‑1 (MQ=255)
                              ggCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctct                                                           <  2:2914696/93‑1 (MQ=255)
                              ggCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCCGTGCTGCCGGTGCGGATACTGCTGGCCCCGCGTGGCTAActctct                                                           >  1:4219335/1‑93 (MQ=255)
                                 gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGCCCCGCGTGACTACCTCTCTAcc                                                        >  2:4735806/1‑93 (MQ=255)
                                  tGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCt                                                       <  1:3777984/93‑1 (MQ=255)
                                    cAGTTGACGTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                                     <  2:2409478/93‑1 (MQ=255)
                                    cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                                     >  2:1604869/1‑93 (MQ=255)
                                      gTTGACCTGGATGCTGGGGGGTTACGGTTTTGTGGCGGCGGTGCTTCCGGTGCGGTTACTCCTGGCCCCGCGTGCCTACCTCTCCACcctcct                                                   >  2:4201451/1‑93 (MQ=255)
                                      gTTGACCTAGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                                   <  1:2755595/93‑1 (MQ=255)
                                       ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTg                                                  >  2:1224111/1‑93 (MQ=255)
                                        tGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGa                                                 >  2:1629398/1‑93 (MQ=255)
                                         gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGGCCCGCGTGACTACCTCTCTACCTTCCTGaa                                                >  1:4287852/1‑93 (MQ=255)
                                         gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                                >  1:5342682/1‑93 (MQ=255)
                                         gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                                <  2:5416158/93‑1 (MQ=255)
                                         gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                                >  2:6712344/1‑93 (MQ=255)
                                          aCCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                               >  1:823454/1‑93 (MQ=255)
                                           ccTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaaa                                              <  1:2534506/93‑1 (MQ=255)
                                           ccTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaaa                                              >  1:4300004/1‑93 (MQ=255)
                                            cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAt                                             >  2:4026901/1‑93 (MQ=255)
                                            cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAt                                             <  1:4543019/93‑1 (MQ=255)
                                            cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAt                                             >  2:3061456/1‑93 (MQ=255)
                                            cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAt                                             >  1:5064218/1‑93 (MQ=255)
                                             tGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            >  2:6045722/1‑93 (MQ=255)
                                              ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGGCTACCTCTCTACCTTCCTGCAAATCg                                           >  1:1416780/1‑93 (MQ=255)
                                              ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCg                                           <  1:2998416/93‑1 (MQ=255)
                                               gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCgg                                          <  1:2092850/93‑1 (MQ=255)
                                                 tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGa                                        >  2:6013727/1‑93 (MQ=255)
                                                 tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGa                                        >  1:2296694/1‑93 (MQ=255)
                                                 tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCGCTCTACCTTCCTGAAAATCGGGa                                        <  1:2912916/93‑1 (MQ=255)
                                                   cTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACt                                      >  1:3939592/1‑93 (MQ=255)
                                                     ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAt                                    >  2:2770899/1‑93 (MQ=255)
                                                     ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAt                                    <  2:4964373/93‑1 (MQ=255)
                                                      gtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtc                                   >  1:4970696/1‑93 (MQ=255)
                                                        gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                 >  1:6064343/1‑93 (MQ=255)
                                                        gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                 <  2:3428546/93‑1 (MQ=255)
                                                           ttACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                              <  1:5512148/93‑1 (MQ=255)
                                                           ttACGGTTTTGTGGCGGCGGTGCTGCCGGGGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                              >  2:6309681/1‑93 (MQ=255)
                                                             aCGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTc                            >  1:5583065/1‑93 (MQ=255)
                                                             aCGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTc                            <  2:3052284/93‑1 (MQ=255)
                                                              cGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                           <  2:5929101/93‑1 (MQ=255)
                                                              cGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                           <  2:1874021/93‑1 (MQ=255)
                                                                gTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTgg                         >  2:450738/1‑93 (MQ=255)
                                                                    tgtgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGt                     <  1:1674884/93‑1 (MQ=255)
                                                                    tgtgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGt                     <  1:4884242/93‑1 (MQ=255)
                                                                     gtgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTa                    <  2:3985460/93‑1 (MQ=255)
                                                                     gtgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTa                    <  1:5866763/93‑1 (MQ=255)
                                                                       ggcggcGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAgg                  <  2:2178327/93‑1 (MQ=255)
                                                                       ggcggcGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAgg                  <  1:5778043/93‑1 (MQ=255)
                                                                        gcggcgGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGc                 >  1:5108169/1‑93 (MQ=255)
                                                                          ggcggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAt               >  1:2360711/1‑93 (MQ=255)
                                                                          ggcggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAt               >  1:6469433/1‑93 (MQ=255)
                                                                           gcggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtt              <  2:415215/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            <  2:6211964/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            <  2:1369312/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            <  2:2480842/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            <  1:3980958/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            <  1:3136058/93‑1 (MQ=255)
                                                                               tgctgcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGa          <  2:6483603/93‑1 (MQ=255)
                                                                                  tgcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTa       >  1:3418343/1‑93 (MQ=255)
                                                                                   gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      <  2:4442995/93‑1 (MQ=255)
                                                                                   gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      >  1:1216389/1‑93 (MQ=255)
                                                                                       gtgcGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgcgc  >  2:5025665/1‑93 (MQ=255)
                                                                                       gtgcGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgcgc  <  1:4180126/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGC  >  NC_000913/629843‑630022

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