Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,697,970 T→C 5.4% W505R (TGG→CGG)  bcsE → cellulose production protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,697,9700TC 5.4% 252.4 / 3.4 74W505R (TGG→CGG) bcsEcellulose production protein
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base C (2/2);  total (37/37)
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.06e-01

GGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTG  >  NC_000913/3697839‑3698112
                                                                                                                                   |                                                                                                                                              
ggTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTgg                                                                                                                                              <  2:576142/134‑1 (MQ=255)
ggTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTgg                                                                                                                                              >  1:576142/1‑134 (MQ=255)
        gatgatCAAATCAGTGCCGAGCTGGTGCAGATGCGTTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGc                                                                                                                        >  4:1080100/1‑148 (MQ=255)
            atCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGa                                                                                                                                >  4:805146/1‑136 (MQ=255)
            atCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGa                                                                                                                                <  3:805146/136‑1 (MQ=255)
             tCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCa                                                                                                                           <  1:681817/140‑1 (MQ=255)
             tCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCa                                                                                                                           >  2:681817/1‑140 (MQ=255)
               aaaTCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCTCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGtt                                                                                                                 <  3:1080100/148‑1 (MQ=255)
                        ccGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGc                                                                                                        >  2:650268/1‑148 (MQ=255)
                           aGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                               <  2:27561/138‑1 (MQ=255)
                           aGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                               >  1:27561/1‑138 (MQ=255)
                             cTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACt                                                                                                                                                <  1:893539/103‑1 (MQ=255)
                             cTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACt                                                                                                                                                >  2:893539/1‑103 (MQ=255)
                                 tGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTGGcg                                                                                                                                            >  1:969105/1‑103 (MQ=255)
                                 tGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCAGTGGcg                                                                                                                                            <  2:969105/103‑1 (MQ=255)
                                 tGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTa                                                                                                                >  1:481613/1‑131 (MQ=255)
                                 tGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTa                                                                                                                <  2:481613/131‑1 (MQ=255)
                                     gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCATAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcga                                                                                                                                           >  2:80850/1‑100 (MQ=255)
                                     gATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCATAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcga                                                                                                                                           <  1:80850/100‑1 (MQ=255)
                                         cgcTTGCTTGCCCCAGAACAGTGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                               >  4:733770/1‑124 (MQ=255)
                                         cgcTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAg                                                                                                               <  3:733770/124‑1 (MQ=255)
                                               cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTgg                                                                                                                                              >  2:1025248/1‑87 (MQ=255)
                                               cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTgg                                                                                                                                              <  1:1025248/87‑1 (MQ=255)
                                               cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATAc                                                                                                                                    <  3:126794/97‑1 (MQ=255)
                                               cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATAc                                                                                                                                    >  4:126794/1‑97 (MQ=255)
                                               cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATg                                                                                                                         <  1:424486/108‑1 (MQ=255)
                                               cttgcCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATg                                                                                                                         >  2:424486/1‑108 (MQ=255)
                                                      cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                  <  1:60608/108‑1 (MQ=255)
                                                      cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                  >  2:60608/1‑108 (MQ=255)
                                                      cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcatc                                                                                         >  2:934412/1‑133 (MQ=255)
                                                      cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtcatc                                                                                         <  1:934412/133‑1 (MQ=255)
                                                      cAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                          >  4:346654/1‑148 (MQ=255)
                                                          aCAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATAcc                                                                                                                                   <  3:873160/87‑1 (MQ=255)
                                                          aCAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATAcc                                                                                                                                   >  4:873160/1‑87 (MQ=255)
                                                               ggggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                  >  1:1044095/1‑99 (MQ=255)
                                                               ggggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                  <  2:1044095/99‑1 (MQ=255)
                                                                 ggCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGa                                                               <  1:650268/148‑1 (MQ=255)
                                                                   cATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAt                                                                  <  2:496195/143‑1 (MQ=255)
                                                                   cATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAt                                                                  >  1:496195/1‑143 (MQ=255)
                                                                              cTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGtt                                                  <  3:346654/148‑1 (MQ=255)
                                                                               tttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATa                                                                                                                                     >  1:1084817/1‑64 (MQ=255)
                                                                               tttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATa                                                                                                                                     <  2:1084817/64‑1 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCt                                                                                             >  3:667382/1‑103 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCt                                                                                             <  4:667382/103‑1 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAtt                                                                 >  1:40836/1‑131 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCAtt                                                                 <  2:40836/131‑1 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACCGGCGACGAATACCAGAccc                                                                                                                             <  3:439032/71‑1 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACCGGCGACGAATACCAGAAcc                                                                                                                             <  1:489235/71‑1 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACCGGCGACGAATACCAGAAcc                                                                                                                             >  4:439032/1‑71 (MQ=255)
                                                                                ttAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACCGGCGACGAATACCAGAAcc                                                                                                                             >  2:489235/1‑71 (MQ=255)
                                                                                      cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTg                                                                                                                   <  4:64940/75‑1 (MQ=255)
                                                                                      cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTg                                                                                                                   >  3:64940/1‑75 (MQ=255)
                                                                                      cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                          >  4:28382/1‑116 (MQ=255)
                                                                                      cAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGc                                                                          <  3:28382/116‑1 (MQ=255)
                                                                                        aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCa                                                                   <  3:252799/121‑1 (MQ=255)
                                                                                        aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCa                                                                   >  4:252799/1‑121 (MQ=255)
                                                                                        aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGaa                                                              >  2:656814/1‑126 (MQ=255)
                                                                                        aaGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGaa                                                              <  1:656814/126‑1 (MQ=255)
                                                                                            tCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcgc                                                                                              <  2:66915/90‑1 (MQ=255)
                                                                                            tCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcgc                                                                                              >  1:66915/1‑90 (MQ=255)
                                                                                            tCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAg                                                                           >  3:511369/1‑109 (MQ=255)
                                                                                            tCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAg                                                                           <  4:511369/109‑1 (MQ=255)
                                                                                               aaaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                  >  1:992602/1‑67 (MQ=255)
                                                                                               aaaCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                                                                                                  <  2:992602/67‑1 (MQ=255)
                                                                                                 aCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcgc                                                                                              <  2:911744/85‑1 (MQ=255)
                                                                                                 aCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcgc                                                                                              >  1:911744/1‑85 (MQ=255)
                                                                                                 aCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGa                                                               >  3:1035992/1‑116 (MQ=255)
                                                                                                 aCCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGa                                                               <  4:1035992/116‑1 (MQ=255)
                                                                                                                    cACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAatta                                                          >  2:382674/1‑102 (MQ=255)
                                                                                                                    cACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAatta                                                          <  1:382674/102‑1 (MQ=255)
                                                                                                                                cACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCa                                            >  2:242680/1‑104 (MQ=255)
                                                                                                                                cACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCa                                            <  1:242680/104‑1 (MQ=255)
                                                                                                                                 aCTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTgg                         >  2:1005345/1‑122 (MQ=255)
                                                                                                                                 aCTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTgg                         <  1:1005345/122‑1 (MQ=255)
                                                                                                                                 aCTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTg  <  4:867279/145‑1 (MQ=255)
                                                                                                                                 aCTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTg  >  3:867279/1‑145 (MQ=255)
                                                                                                                                   |                                                                                                                                              
GGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTG  >  NC_000913/3697839‑3698112

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: