Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,256 T→A 7.2% I272N (ATT→AAT)  cmoB → tRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,954,2560TA 7.2% 370.5 / 14.4 107I272N (ATT→AAT) cmoBtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (50/49);  new base A (4/4);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 107 total aligned reads displayed.

GAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCG  >  NC_000913/1954114‑1954401
                                                                                                                                              |                                                                                                                                                 
gAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                             >  2:201122/1‑149 (MQ=255)
 aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                             >  4:325243/1‑148 (MQ=255)
 aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAt                                                                                                                                            <  1:201122/149‑1 (MQ=255)
 aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAt                                                                                                                                            <  1:78688/149‑1 (MQ=255)
  aaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtg                                                                                                                                           >  1:58596/1‑149 (MQ=255)
   aCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgt                                                                                                                                          <  2:782925/149‑1 (MQ=255)
   aCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgt                                                                                                                                          <  2:28581/149‑1 (MQ=255)
     gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        <  2:494538/149‑1 (MQ=255)
     gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        <  1:444587/149‑1 (MQ=255)
      cTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        >  4:244655/1‑148 (MQ=255)
      cTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAg                                                                                                                                       >  2:725554/1‑149 (MQ=255)
       tGGTTATTGATGGTGATGAGAACACTGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAg                                                                                                                                       >  4:903708/1‑148 (MQ=255)
       tGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGc                                                                                                                                      <  2:1007775/149‑1 (MQ=255)
        ggTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGc                                                                                                                                      >  4:191228/1‑148 (MQ=255)
          ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGtt                                                                                                                                   >  1:472611/1‑149 (MQ=255)
          ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGt                                                                                                                                    <  3:133170/148‑1 (MQ=255)
            aTTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTac                                                                                                                                 <  1:725554/149‑1 (MQ=255)
               gATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccac                                                                                                                              <  2:133802/149‑1 (MQ=255)
                  ggTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTAccacca                                                                                                                            >  3:794726/1‑148 (MQ=255)
                    tgatgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGCAGATATCCGCATTGCGGATGTGAGCGTTACCACcaca                                                                                                                          >  3:174581/1‑148 (MQ=255)
                      atgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaa                                                                                                                       <  2:570522/149‑1 (MQ=255)
                       tgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaa                                                                                                                       >  4:1011851/1‑148 (MQ=255)
                          gaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACcaca                                                                                                                          <  2:967085/142‑1 (MQ=255)
                          gaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACcaca                                                                                                                          >  1:967085/1‑142 (MQ=255)
                           aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                                  >  2:1053119/1‑149 (MQ=255)
                           aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                                  >  2:139660/1‑149 (MQ=255)
                            acacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcag                                                                                                                 <  1:761942/149‑1 (MQ=255)
                            acacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                                  >  4:1026240/1‑148 (MQ=255)
                                 gTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGAc                                                                                                             >  3:799064/1‑148 (MQ=255)
                                 gTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGAc                                                                                                             >  4:173332/1‑148 (MQ=255)
                                  tGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACg                                                                                                            <  3:907493/148‑1 (MQ=255)
                                        tACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCCGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                     >  2:860544/1‑149 (MQ=255)
                                            aGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCgg                                                                                                                                              <  4:585323/104‑1 (MQ=255)
                                            aGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCgg                                                                                                                                              >  3:585323/1‑104 (MQ=255)
                                             ggCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggat                                                                                                <  1:415563/149‑1 (MQ=255)
                                               cGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggatgg                                                                                              <  1:902203/149‑1 (MQ=255)
                                                gATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGt                                                                                             >  2:180660/1‑149 (MQ=255)
                                                    gTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTAcc                                                                                         <  2:58596/149‑1 (MQ=255)
                                                       aCGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGaa                                                                                      <  2:573649/149‑1 (MQ=255)
                                                           gcAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCAc                                                                                                         >  1:1053641/1‑126 (MQ=255)
                                                           gcAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCAc                                                                                                         <  2:1053641/126‑1 (MQ=255)
                                                           gcAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                                  >  2:341599/1‑149 (MQ=255)
                                                            cAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                                  >  4:162441/1‑148 (MQ=255)
                                                             aaaTGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                                <  2:784388/149‑1 (MQ=255)
                                                              aaTGTGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTAGTTTTGTAGATATCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                                >  4:903302/1‑148 (MQ=255)
                                                                   gTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAt                                                                          <  1:341599/149‑1 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGtt                                                                                                                                   <  1:1035237/88‑1 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTcc                                                                      <  1:139660/149‑1 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTc                                                                       >  3:787481/1‑148 (MQ=255)
                                                                        gTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTcc                                                                      >  4:5243/1‑148 (MQ=255)
                                                                         tCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                <  2:1058450/73‑1 (MQ=255)
                                                                         tCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                >  1:1058450/1‑73 (MQ=255)
                                                                          cTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTc                                                                    >  4:188793/1‑148 (MQ=255)
                                                                          cTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTc                                                                    >  4:448787/1‑148 (MQ=255)
                                                                           tATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGa                                                                  >  2:54333/1‑149 (MQ=255)
                                                                            aTTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGa                                                                  >  4:419876/1‑148 (MQ=255)
                                                                              ttCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGAcc                                                                >  3:482919/1‑148 (MQ=255)
                                                                              ttCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGAcc                                                                >  3:45170/1‑148 (MQ=255)
                                                                               tCATCCCTTCCGCGCTGGCGCTGAAAAACTGGATGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATTGTTACCGAATCGCTGGCCGATTTCCTCGACCca                                                              <  2:632503/149‑1 (MQ=255)
                                                                                cATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCca                                                              <  4:787481/148‑1 (MQ=255)
                                                                                   cccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATg                                                          >  1:291852/1‑149 (MQ=255)
                                                                                         cgcgCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCgg                                                    <  1:155568/149‑1 (MQ=255)
                                                                                         cgcgCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCgg                                                    <  2:331447/149‑1 (MQ=255)
                                                                                         cgcgCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCgg                                                    <  1:824741/149‑1 (MQ=255)
                                                                                         cgcgCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCg                                                     <  4:799064/148‑1 (MQ=255)
                                                                                               ggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaa                                               <  3:903708/148‑1 (MQ=255)
                                                                                                    tgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTg                                         >  2:906459/1‑149 (MQ=255)
                                                                                                    tgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTg                                         <  1:394353/149‑1 (MQ=255)
                                                                                                    tgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTACGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                          <  3:903302/148‑1 (MQ=255)
                                                                                                       aaaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                      <  2:787554/149‑1 (MQ=255)
                                                                                                       aaaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGa                                       <  3:86061/148‑1 (MQ=255)
                                                                                                        aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAg                                     <  1:439987/149‑1 (MQ=255)
                                                                                                              gCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTAt                                <  4:482919/148‑1 (MQ=255)
                                                                                                                 gaagaaGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGGAAAACGGTGGAAGGTTATCCg                             >  3:95613/1‑148 (MQ=255)
                                                                                                                  aagaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                           >  2:267435/1‑149 (MQ=255)
                                                                                                                  aagaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                           >  2:1009211/1‑149 (MQ=255)
                                                                                                                    gaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCAc                         <  2:947610/149‑1 (MQ=255)
                                                                                                                     aagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCAc                         >  4:515749/1‑148 (MQ=255)
                                                                                                                      agTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCt                       <  1:906459/149‑1 (MQ=255)
                                                                                                                            gTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgc                 >  2:424736/1‑149 (MQ=255)
                                                                                                                              tttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgcg                <  3:5243/148‑1 (MQ=255)
                                                                                                                               ttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgcgc               <  4:706811/148‑1 (MQ=255)
                                                                                                                                  tAGATATCCGCATTGCGGATGTGAGCGTAACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGt            >  3:859856/1‑148 (MQ=255)
                                                                                                                                   tcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTATCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                          >  1:400485/5‑117 (MQ=255)
                                                                                                                                   tcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTATCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                          <  2:400485/113‑1 (MQ=255)
                                                                                                                                   aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTg           >  3:282450/1‑148 (MQ=255)
                                                                                                                                   aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTg           >  3:373125/1‑148 (MQ=255)
                                                                                                                                    gATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCt         <  1:424736/149‑1 (MQ=255)
                                                                                                                                     acatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGa                                                                  <  1:943680/89‑1 (MQ=255)
                                                                                                                                     acatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc          <  1:828897/145‑1 (MQ=255)
                                                                                                                                     acatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc          >  2:828897/3‑147 (MQ=255)
                                                                                                                                      tatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGa       <  1:826571/149‑1 (MQ=255)
                                                                                                                                        tCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGAt      <  3:419876/148‑1 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTg    >  2:755741/1‑149 (MQ=255)
                                                                                                                                          cGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgc   >  2:512445/1‑149 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGa                                                         >  3:74753/1‑94 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGa                                                         <  4:74753/94‑1 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcg  >  1:873049/1‑149 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcg  >  2:512111/1‑149 (MQ=255)
                                                                                                                                            cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcg  >  3:372394/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCG  >  NC_000913/1954114‑1954401

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