Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,256 T→A 5.4% 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 5.4% 456.8 / 12.2 130I272N (ATT→AAT) cmoBtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (62/61);  new base A (3/4);  total (65/65)
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.04e-01

Only 100 of 135 total aligned reads displayed.

GTGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGC  >  NC_000913/1954108‑1954402
                                                                                                                                                    |                                                                                                                                                  
gTGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAt                                                                                                                                                    >  1:484290/1‑149 (MQ=255)
 tGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAtt                                                                                                                                                   <  2:71284/149‑1 (MQ=255)
   cTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                 <  3:286205/149‑1 (MQ=255)
    tGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTGTTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCg                                                                                                                                                >  2:922011/1‑149 (MQ=255)
       aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCTGGCGATCGTTACGCGCAAATGCGTAGTGTCTATCTCATCCCTTCCGCGCTGGCCCTGAAAAACTGGCTGAAGAAGTGTGTTTTTGTAGATATCCGCATTGCGTAt                                                                                                                                             <  4:777020/149‑1 (MQ=255)
         aCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgt                                                                                                                                           <  3:550082/149‑1 (MQ=255)
           gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAtatct                                                                                                                                         <  4:36321/149‑5 (MQ=255)
           gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAtatct                                                                                                                                         <  3:265655/149‑5 (MQ=255)
            cTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAg                                                                                                                                        >  4:561554/1‑149 (MQ=255)
              ggTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCg                                                                                                                                      >  4:375168/1‑149 (MQ=255)
                ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGtt                                                                                                                                    <  2:864363/149‑1 (MQ=255)
                ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGtt                                                                                                                                    <  2:717142/149‑1 (MQ=255)
                   ttGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacc                                                                                                                                 >  4:826749/1‑149 (MQ=255)
                   ttGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacc                                                                                                                                 >  4:295594/1‑149 (MQ=255)
                    tGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacca                                                                                                                                <  4:549963/149‑1 (MQ=255)
                       tGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTAccacca                                                                                                                             >  3:309922/1‑149 (MQ=255)
                       tGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTAccacca                                                                                                                             <  4:514647/149‑1 (MQ=255)
                       tGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTAccacca                                                                                                                             <  4:437877/149‑1 (MQ=255)
                         gTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACcaca                                                                                                                           <  2:520299/149‑1 (MQ=255)
                          tgatgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACag                                                                                                                          >  2:829457/1‑149 (MQ=255)
                           gatgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACaga                                                                                                                         >  4:411810/1‑149 (MQ=255)
                            atgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTGTTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaa                                                                                                                        <  1:922011/149‑1 (MQ=255)
                            atgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaa                                                                                                                        <  1:687312/149‑1 (MQ=255)
                            atgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaa                                                                                                                        <  4:639312/149‑1 (MQ=255)
                            atgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaa                                                                                                                        <  4:175573/149‑1 (MQ=255)
                              gagaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaaga                                                                                                                      <  3:28832/149‑1 (MQ=255)
                                gaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagc                                                                                                                    >  3:104534/1‑149 (MQ=255)
                                gaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGAGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagc                                                                                                                    >  2:234469/1‑149 (MQ=255)
                                 aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                                   >  2:485540/1‑149 (MQ=255)
                                 aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAt                                                                                                                                             >  3:703858/1‑123 (MQ=255)
                                 aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAt                                                                                                                                             <  4:703858/123‑1 (MQ=255)
                                  acacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcag                                                                                                                  >  1:498302/1‑149 (MQ=255)
                                   cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACTGAAGagcagc                                                                                                                 <  2:726011/149‑1 (MQ=255)
                                   cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAATagcagc                                                                                                                 <  3:201584/149‑1 (MQ=255)
                                   cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcagc                                                                                                                 >  3:573032/1‑149 (MQ=255)
                                      ggTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGAc                                                                                                              <  2:627614/149‑1 (MQ=255)
                                       gTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACg                                                                                                             >  3:877160/1‑149 (MQ=255)
                                           tGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCAcc                                                                                                         >  1:581259/1‑149 (MQ=255)
                                             gTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                 >  3:402027/1‑107 (MQ=255)
                                             gTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGa                                                                                                       >  4:627224/1‑149 (MQ=255)
                                              tACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                      >  1:263665/1‑149 (MQ=255)
                                              tACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                      <  3:826749/149‑1 (MQ=255)
                                               aCCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAat                                                                                                     <  3:451923/149‑1 (MQ=255)
                                                  aGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgga                                                                                                  >  2:692591/1‑149 (MQ=255)
                                                   ggCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggat                                                                                                 >  1:693573/1‑149 (MQ=255)
                                                   ggCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggat                                                                                                 <  4:124632/149‑1 (MQ=255)
                                                   ggCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggat                                                                                                 >  4:517053/1‑149 (MQ=255)
                                                      gATCGTTACGCGCAAATGCGTAATGCCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGt                                                                                              >  3:478819/1‑149 (MQ=255)
                                                        tCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTa                                                                                            >  2:651219/1‑149 (MQ=255)
                                                         cGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTAc                                                                                           <  2:380200/149‑1 (MQ=255)
                                                          gTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTAcc                                                                                          >  1:140000/1‑149 (MQ=255)
                                                          gTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTAcc                                                                                          <  4:881341/149‑1 (MQ=255)
                                                           ttACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggat                                                                                                 <  3:450326/141‑1 (MQ=255)
                                                             aCGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGaa                                                                                       <  4:600334/149‑1 (MQ=255)
                                                             aCGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGaa                                                                                       >  2:472439/1‑149 (MQ=255)
                                                                 gcAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                                   <  3:81178/149‑1 (MQ=255)
                                                                    aaTGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTgg                                                                                <  4:901189/149‑1 (MQ=255)
                                                                       gCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCg                                                                             >  4:674349/1‑149 (MQ=255)
                                                                         gTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAt                                                                           >  1:13662/1‑149 (MQ=255)
                                                                               tCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTc                                                                     >  1:861608/1‑149 (MQ=255)
                                                                                   tttCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGAcc                                                                 <  4:104534/149‑1 (MQ=255)
                                                                                       aTCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcag                                                                                                                  >  2:126824/1‑96 (MQ=255)
                                                                                       aTCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcag                                                                                                                  <  1:126824/96‑1 (MQ=255)
                                                                                        tCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggatgg                                                                                               <  3:644302/114‑1 (MQ=255)
                                                                                        tCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACAt                                                            >  3:71012/1‑149 (MQ=255)
                                                                                        tCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACAt                                                            >  1:168734/1‑149 (MQ=255)
                                                                                        tCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGGTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggatgg                                                                                               >  4:644302/1‑114 (MQ=255)
                                                                                         cccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATg                                                           >  4:718250/1‑149 (MQ=255)
                                                                                            ttCCGCGCTGGCGCGGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATc                                                        >  3:69794/1‑149 (MQ=255)
                                                                                              ccGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGGGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTAACACCACAGAAGAGCAGGGGCGCACCGAATGGATGGTTACCGAATCGCTCGGCGGTTTCCTCGACCCACATGATCCg                                                      >  2:166195/1‑149 (MQ=255)
                                                                                                    tggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaa                                                <  2:636912/149‑1 (MQ=255)
                                                                                                    tggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaa                                                >  4:231866/1‑149 (MQ=255)
                                                                                                     ggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                               <  1:367996/149‑1 (MQ=255)
                                                                                                     ggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                               >  3:192597/1‑149 (MQ=255)
                                                                                                         ctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                           <  1:460364/149‑1 (MQ=255)
                                                                                                         ctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                           <  3:298066/149‑1 (MQ=255)
                                                                                                             aaaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                       <  4:291226/149‑1 (MQ=255)
                                                                                                             aaaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                       <  1:829457/149‑1 (MQ=255)
                                                                                                             aaaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                       <  3:626967/149‑1 (MQ=255)
                                                                                                              aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAg                                      <  1:234469/149‑1 (MQ=255)
                                                                                                              aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAg                                      <  4:309922/149‑1 (MQ=255)
                                                                                                                 cTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                                   <  2:276755/149‑1 (MQ=255)
                                                                                                                    gCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATc                                <  1:835287/149‑1 (MQ=255)
                                                                                                                        aagaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGCCCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                            <  4:461704/149‑1 (MQ=255)
                                                                                                                        aagaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                            >  4:561009/1‑149 (MQ=255)
                                                                                                                              tgtgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTaa                      >  4:160938/1‑149 (MQ=255)
                                                                                                                              tgtgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTaa                      >  1:496805/1‑149 (MQ=255)
                                                                                                                                tgGTTTTTTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCgg                                                     >  2:507593/1‑116 (MQ=255)
                                                                                                                                tgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAAt                                                                                      <  2:556613/83‑1 (MQ=255)
                                                                                                                                tgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCgg                                                     <  1:507593/116‑1 (MQ=255)
                                                                                                                                  gTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgc                  >  1:666362/1‑149 (MQ=255)
                                                                                                                                     ttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgcgcg               <  2:484290/149‑1 (MQ=255)
                                                                                                                                       gTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGt             >  2:305386/1‑149 (MQ=255)
                                                                                                                                         tcacatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                               >  3:541751/5‑113 (MQ=255)
                                                                                                                                         tcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                               <  4:541751/109‑1 (MQ=255)
                                                                                                                                         tcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc           >  3:185837/5‑149 (MQ=255)
                                                                                                                                            tatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGa        <  3:95316/149‑1 (MQ=255)
                                                                                                                                                  cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcgc  <  4:478819/149‑1 (MQ=255)
                                                                                                                                                    ttGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATcc                                                       <  3:911494/94‑1 (MQ=255)
                                                                                                                                                    ttGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATcc                                                       >  4:911494/1‑94 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGC  >  NC_000913/1954108‑1954402

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