Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,697,969 C→G 5.7% H504Q (CAC→CAG bcsE → cellulose production protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,697,9690CG 5.7% 308.0 / 6.6 86H504Q (CAC→CAGbcsEcellulose production protein
Reads supporting (aligned to +/- strand):  ref base C (36/45);  new base G (2/3);  total (38/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 = 9.94e-01

CTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGA  >  NC_000913/3697822‑3698115
                                                                                                                                                   |                                                                                                                                                  
ctcAAACCGTAAGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACt                                                                                                                                                   >  3:362355/1‑149 (MQ=255)
 tcAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTg                                                                                                                                                  >  2:350027/1‑149 (MQ=255)
 tcAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTg                                                                                                                                                  <  3:1490/149‑1 (MQ=255)
 tcAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTg                                                                                                                                                  <  2:311251/149‑1 (MQ=255)
 tcAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTg                                                                                                                                                  >  2:401817/1‑149 (MQ=255)
 tcAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTg                                                                                                                                                  >  3:251613/1‑149 (MQ=255)
   aaaCCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGc                                                                                                                                                >  3:309220/1‑149 (MQ=255)
     aCCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcga                                                                                                                                              >  4:263938/1‑149 (MQ=255)
      ccGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGATTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcgac                                                                                                                                             <  3:587321/149‑1 (MQ=255)
       cGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTGGcgacc                                                                                                                                            <  2:363924/149‑2 (MQ=255)
        gTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcgacga                                                                                                                                           >  4:171938/1‑149 (MQ=255)
          aTGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTGGcgaccatc                                                                                                                                         <  1:350027/149‑5 (MQ=255)
          aTGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAAt                                                                                                                                         <  2:277633/149‑1 (MQ=255)
             gTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATAcc                                                                                                                                      <  1:133479/149‑1 (MQ=255)
              tCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCa                                                                                                                                     <  1:357725/149‑1 (MQ=255)
                 ggTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGaa                                                                                                                                  <  3:389115/149‑1 (MQ=255)
                   tttGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAAcc                                                                                                                                >  4:519799/1‑149 (MQ=255)
                      gAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATAGTCGCCACTGGCGACGAATACCTGAACCCAt                                                                                                                             <  1:439078/149‑1 (MQ=255)
                       aaGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATg                                                                                                                            <  4:11703/149‑1 (MQ=255)
                               cAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                    <  1:341989/149‑1 (MQ=255)
                               cAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                    <  3:552706/149‑1 (MQ=255)
                               cAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                    <  3:603010/149‑1 (MQ=255)
                                 aaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  >  1:400582/1‑149 (MQ=255)
                                 aaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  <  2:426958/149‑1 (MQ=255)
                                 aaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  >  1:142092/1‑149 (MQ=255)
                                    cAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgatg                                                                                                               >  1:227565/1‑149 (MQ=255)
                                      gTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgatgat                                                                                                             <  2:400582/149‑1 (MQ=255)
                                          cGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCtg                                                                                                         <  1:299634/149‑1 (MQ=255)
                                           gAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                     >  3:269111/1‑136 (MQ=255)
                                           gAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                     <  4:269111/136‑1 (MQ=255)
                                                 gTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcg                                                                                                  <  4:309220/149‑1 (MQ=255)
                                                    cAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtc                                                                                               >  4:256182/1‑149 (MQ=255)
                                                           gcttgcttGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgat                                                                                                                <  1:125603/125‑1 (MQ=255)
                                                           gcttgcttGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatga                                                                                        >  3:234644/1‑149 (MQ=255)
                                                           gcttgcttGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACAGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgat                                                                                                                >  2:125603/1‑125 (MQ=255)
                                                            cttgcttgCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatgat                                                                                       >  3:60267/1‑149 (MQ=255)
                                                              tgcttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCAtgatga                                                                                     <  4:362355/149‑1 (MQ=255)
                                                                cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGAcc                                                                                   <  2:304844/149‑1 (MQ=255)
                                                                cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGAcc                                                                                   <  3:475147/149‑1 (MQ=255)
                                                                cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGAcc                                                                                   <  1:501133/149‑1 (MQ=255)
                                                                    ccccAGAACAATGGGGCATGCCGCTGCATTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCGCGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCa                                                                               <  2:294063/149‑1 (MQ=255)
                                                                     cccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                    >  3:178082/1‑111 (MQ=255)
                                                                     cccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                    <  4:178082/111‑1 (MQ=255)
                                                                       cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCg                                                                            >  2:520237/1‑149 (MQ=255)
                                                                          aaCAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGata                                                                         >  1:129365/1‑149 (MQ=255)
                                                                              aTGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAt                                                                     <  2:129365/149‑1 (MQ=255)
                                                                               tGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAtt                                                                    <  4:119287/149‑1 (MQ=255)
                                                                                 gggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGa                                                                  >  1:412074/1‑149 (MQ=255)
                                                                                  ggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGaa                                                                 <  3:263938/149‑1 (MQ=255)
                                                                                      tGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAatta                                                             <  3:148194/149‑1 (MQ=255)
                                                                                       gCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAattat                                                            <  4:251613/149‑1 (MQ=255)
                                                                                       gCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAattat                                                            >  4:67253/1‑149 (MQ=255)
                                                                                         cGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTg                                                          >  1:273718/1‑149 (MQ=255)
                                                                                         cGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTg                                                          >  1:427703/1‑149 (MQ=255)
                                                                                         cGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTg                                                          >  2:165835/1‑149 (MQ=255)
                                                                                                tttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTg                                                   >  2:596628/1‑149 (MQ=255)
                                                                                                tttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTg                                                   >  4:56774/1‑149 (MQ=255)
                                                                                                     cGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                             >  2:373630/1‑118 (MQ=255)
                                                                                                     cGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                             <  1:373630/118‑1 (MQ=255)
                                                                                                     cGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCAc                                              <  2:273718/149‑1 (MQ=255)
                                                                                                      gCAAAGTTCTAAACCGGTCCTCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACt                                             >  4:525942/1‑149 (MQ=255)
                                                                                                      gCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatg                                                                                         <  1:174335/105‑1 (MQ=255)
                                                                                                      gCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatg                                                                                         >  2:174335/1‑105 (MQ=255)
                                                                                                      gCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACt                                             >  3:60741/1‑149 (MQ=255)
                                                                                                       cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATCCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTg                                            >  1:392812/1‑149 (MQ=255)
                                                                                                         aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGat                                          <  2:588244/149‑1 (MQ=255)
                                                                                                           gTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGatat                                        <  3:56774/149‑1 (MQ=255)
                                                                                                           gTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGatat                                        >  2:271238/1‑149 (MQ=255)
                                                                                                             tCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATAttt                                      <  1:403560/149‑1 (MQ=255)
                                                                                                             tCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATAttt                                      >  3:466987/1‑149 (MQ=255)
                                                                                                                aaaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCAtgatga                                                                                     <  2:347042/99‑1 (MQ=255)
                                                                                                                aaaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCAtgatga                                                                                     >  1:347042/1‑99 (MQ=255)
                                                                                                                aaaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTc                                   <  3:171938/149‑1 (MQ=255)
                                                                                                                 aaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTcc                                  <  1:466363/149‑1 (MQ=255)
                                                                                                                  aCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTccc                                 <  2:205298/149‑1 (MQ=255)
                                                                                                                     ggTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCt                              <  1:271238/149‑1 (MQ=255)
                                                                                                                      gTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTg                             <  2:412074/149‑1 (MQ=255)
                                                                                                                       tcatcaATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTgg                            <  1:596628/149‑1 (MQ=255)
                                                                                                                       tcatcaATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTgg                            <  2:142092/149‑1 (MQ=255)
                                                                                                                             aTGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATc                      <  1:520237/149‑1 (MQ=255)
                                                                                                                                  gAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGgcg                 <  4:3369/149‑1 (MQ=255)
                                                                                                                                      aCGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatgat                                                                                       >  3:415817/1‑75 (MQ=255)
                                                                                                                                      aCGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatgat                                                                                       <  4:415817/75‑1 (MQ=255)
                                                                                                                                       tattcgTCGCCAGTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCa            >  4:524763/6‑149 (MQ=255)
                                                                                                                                             tCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCttt      >  2:275839/1‑149 (MQ=255)
                                                                                                                                                 cACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGa  >  3:563050/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGA  >  NC_000913/3697822‑3698115

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: