Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,697,969 C→G 5.5% 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.5% 454.8 / 12.3 128H504Q (CAC→CAGbcsEcellulose production protein
Reads supporting (aligned to +/- strand):  ref base C (60/61);  new base G (3/4);  total (63/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 = 9.96e-01

Only 100 of 128 total aligned reads displayed.

CAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGA  >  NC_000913/3697824‑3698115
                                                                                                                                                 |                                                                                                                                                  
cAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTg                                                                                                                                                  >  3:559819/1‑148 (MQ=255)
cAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTg                                                                                                                                                  <  3:486286/148‑1 (MQ=255)
  aaCCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGc                                                                                                                                                >  4:735719/1‑148 (MQ=255)
   aCCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcga                                                                                                                                              <  1:969714/149‑1 (MQ=255)
   aCCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcg                                                                                                                                               >  4:732906/1‑148 (MQ=255)
     cGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCCTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTGGcgac                                                                                                                                             <  4:343498/148‑1 (MQ=255)
      gTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTGGcgacca                                                                                                                                           <  2:532350/149‑3 (MQ=255)
       tATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcgacga                                                                                                                                           <  4:745022/148‑1 (MQ=255)
            tCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATAcc                                                                                                                                      <  4:194389/148‑1 (MQ=255)
            tCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCa                                                                                                                                     <  2:188481/149‑1 (MQ=255)
              tGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGa                                                                                                                                   >  2:1002463/1‑149 (MQ=255)
               ggTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGa                                                                                                                                   <  4:477502/148‑1 (MQ=255)
                 tttGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAAc                                                                                                                                 >  3:302371/1‑148 (MQ=255)
                      aGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATg                                                                                                                            <  4:329550/148‑1 (MQ=255)
                       gatgatCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGc                                                                                                                           >  3:161293/1‑148 (MQ=255)
                         tgatCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGa                                                                                                                         >  4:547244/1‑148 (MQ=255)
                           atCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTg                                                                                                                      >  2:482039/1‑149 (MQ=255)
                            tCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTg                                                                                                                      <  4:920238/148‑1 (MQ=255)
                            tCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                     <  1:154182/149‑1 (MQ=255)
                             cAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                    <  1:403377/149‑1 (MQ=255)
                              aaaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACt                                                                                                                       >  2:705073/1‑145 (MQ=255)
                              aaaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACt                                                                                                                       <  1:705073/145‑1 (MQ=255)
                               aaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTa                                                                                                                   >  3:351010/1‑148 (MQ=255)
                               aaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  <  1:482039/149‑1 (MQ=255)
                               aaTCAGTGCAGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCAACTGGCGACGAATACCAGAACCCATGCGACTGTTa                                                                                                                   >  4:688712/1‑148 (MQ=255)
                                aTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  <  3:214852/148‑1 (MQ=255)
                                 tCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgat                                                                                                                <  2:980239/149‑1 (MQ=255)
                                 tCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAga                                                                                                                 <  3:871807/148‑1 (MQ=255)
                                   aGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgatga                                                                                                              <  1:34508/149‑1 (MQ=255)
                                    gTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAgatgat                                                                                                             <  1:208467/149‑1 (MQ=255)
                                          aGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCtgtg                                                                                                       <  1:506933/149‑1 (MQ=255)
                                            cTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCa                                                                                                                              >  3:696075/1‑124 (MQ=255)
                                            cTGGTGCAGATGCGCTTGCTCGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCa                                                                                                                              <  4:696075/124‑1 (MQ=255)
                                               gTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCTACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcg                                                                                                  >  1:497150/1‑149 (MQ=255)
                                               gTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcg                                                                                                  >  1:284712/1‑149 (MQ=255)
                                                 gCAGATGCGCTTGCTTGCCCCAGAACTATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcgc                                                                                                 <  4:161293/148‑1 (MQ=255)
                                                   aGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtca                                                                                              >  2:950850/1‑149 (MQ=255)
                                                    gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTg                                                                                                                                                  >  4:772872/1‑96 (MQ=255)
                                                    gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcat                                                                                             >  2:505829/1‑149 (MQ=255)
                                                    gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcat                                                                                             >  1:169209/1‑149 (MQ=255)
                                                    gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcat                                                                                             >  1:376111/1‑149 (MQ=255)
                                                    gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtca                                                                                              >  4:553058/1‑148 (MQ=255)
                                                     aTGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcatc                                                                                            >  2:392208/1‑149 (MQ=255)
                                                         gcttgcttGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatga                                                                                        <  2:993627/149‑1 (MQ=255)
                                                          cttgcttgCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatgat                                                                                       >  2:6557/1‑149 (MQ=255)
                                                                tgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAATGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCa                                                                                  >  4:1000000/1‑148 (MQ=255)
                                                                    ccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatga                                                                                        >  2:855500/1‑138 (MQ=255)
                                                                    ccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatga                                                                                        <  1:855500/138‑1 (MQ=255)
                                                                    ccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAg                                                                              >  3:567882/1‑148 (MQ=255)
                                                                    ccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAg                                                                              >  3:614967/1‑148 (MQ=255)
                                                                    ccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                             >  2:92306/1‑149 (MQ=255)
                                                                    ccAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                             <  2:284712/149‑1 (MQ=255)
                                                                     cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCAtgatga                                                                                     <  4:604387/140‑1 (MQ=255)
                                                                     cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCg                                                                            >  2:145242/1‑149 (MQ=255)
                                                                      aGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGGCGCCACTGGCGACGAATACCAGAACCCATGCGATTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCg                                                                            >  4:455470/1‑148 (MQ=255)
                                                                              ggggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAtt                                                                    >  3:886463/1‑148 (MQ=255)
                                                                               gggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGa                                                                  >  1:813710/1‑149 (MQ=255)
                                                                                 gCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCGTGATGACCATCAGCGATATCATTGaaa                                                                <  1:112710/149‑1 (MQ=255)
                                                                                 gCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGaaa                                                                <  1:743900/149‑1 (MQ=255)
                                                                                  cATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGaaa                                                                >  3:362202/1‑148 (MQ=255)
                                                                                    tGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAatt                                                              <  4:302371/148‑1 (MQ=255)
                                                                                      ccGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAattat                                                            <  3:11616/148‑1 (MQ=255)
                                                                                        gCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTg                                                          <  3:701654/148‑1 (MQ=255)
                                                                                         cTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGt                                                         <  3:563791/148‑1 (MQ=255)
                                                                                         cTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTc                                                        <  1:397373/149‑1 (MQ=255)
                                                                                              tttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTg                                                   >  1:830799/1‑149 (MQ=255)
                                                                                                  aCGCAAAGTTCTAAACCGGTCATCAATGTCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAt                                                                     >  4:448840/1‑127 (MQ=255)
                                                                                                  aCGCAAAGTTCTAAACCGGTCATCAATGTCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAt                                                                     <  3:448840/127‑1 (MQ=255)
                                                                                                  aCGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTgcgc                                                >  3:671586/1‑148 (MQ=255)
                                                                                                    gCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCAc                                              <  3:756950/148‑1 (MQ=255)
                                                                                                     cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCtgt                                                                                                        <  4:604598/89‑1 (MQ=255)
                                                                                                     cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCtgt                                                                                                        >  3:604598/1‑89 (MQ=255)
                                                                                                     cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcat                                                                                             >  4:587803/1‑100 (MQ=255)
                                                                                                     cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcat                                                                                             <  3:587803/100‑1 (MQ=255)
                                                                                                     cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACt                                             >  4:27784/1‑148 (MQ=255)
                                                                                                     cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACt                                             >  4:747219/1‑148 (MQ=255)
                                                                                                      aaaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTg                                            <  4:614967/148‑1 (MQ=255)
                                                                                                       aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTg                                                          >  4:403831/1‑133 (MQ=255)
                                                                                                       aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTg                                                          <  3:403831/133‑1 (MQ=255)
                                                                                                          ttCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATAtt                                       >  2:653772/1‑149 (MQ=255)
                                                                                                             tAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATAtttt                                     >  3:522092/1‑148 (MQ=255)
                                                                                                              aaaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATAttttt                                    >  4:852523/1‑148 (MQ=255)
                                                                                                               aaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAt                                                                     <  1:237008/114‑1 (MQ=255)
                                                                                                                aCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTcc                                  >  3:726737/1‑148 (MQ=255)
                                                                                                                  cGGTCATGAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGc                               >  2:582000/1‑149 (MQ=255)
                                                                                                                  cGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCg                                >  3:830765/1‑148 (MQ=255)
                                                                                                                  cGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGc                               >  2:635373/1‑149 (MQ=255)
                                                                                                                   ggTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGc                               <  3:619296/148‑1 (MQ=255)
                                                                                                                     tcatcaATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTgg                            <  1:55155/149‑1 (MQ=255)
                                                                                                                        tcaATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGatat                                                                        <  4:538603/102‑1 (MQ=255)
                                                                                                                        tcaATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGatat                                                                        >  3:538603/1‑102 (MQ=255)
                                                                                                                        tcaATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCt                         <  2:300023/149‑1 (MQ=255)
                                                                                                                          aaTGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTa                        <  3:770228/148‑1 (MQ=255)
                                                                                                                                gAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGgcg                 >  1:944514/1‑149 (MQ=255)
                                                                                                                                     cGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  >  2:314843/1‑47 (MQ=255)
                                                                                                                                     cGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                                  <  1:314843/47‑1 (MQ=255)
                                                                                                                                        tcgTCGCCAGTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCActc         >  2:888397/3‑149 (MQ=255)
                                                                                                                                           tCGCCACTGGCGACGAATACCAGAACCCATGGGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCtt       <  4:704941/148‑1 (MQ=255)
                                                                                                                                            cGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCttt      <  3:680645/148‑1 (MQ=255)
                                                                                                                                                aCTGGCTACGAATACCAGAGCCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGa  <  4:157751/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGA  >  NC_000913/3697824‑3698115

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