Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,697,854 A→C Q466P (CAA→CCA)  bcsE → cellulose production protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,697,8540AC98.7% 551.8 / ‑4.3 154Q466P (CAA→CCA) bcsEcellulose production protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (78/74);  total (80/74)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

TGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGC  >  NC_000913/3697757‑3697950
                                                                                                 |                                                                                                
tGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCcaa                                                                                                <  1:944777/100‑1 (MQ=255)
 gTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAAt                                                                                               >  2:2986043/1‑100 (MQ=255)
 gTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAAt                                                                                               >  1:2632956/1‑100 (MQ=255)
 gTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAAt                                                                                               >  2:4260/1‑100 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCTCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  2:3826063/100‑1 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  1:2895713/100‑1 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  1:3316187/100‑1 (MQ=255)
   cGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCa                                                                                             <  1:2145284/100‑1 (MQ=255)
    ggATTAACGATCTGGATACCGCGTTGAATCATATTCTCCCATTGCCTACTGGCGACATCTTCTCAAACCGTCTGGTCTGGTTTGAAGATGATCCAATCAg                                                                                            >  2:1555154/1‑100 (MQ=255)
      aTTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTg                                                                                          >  1:46487/1‑100 (MQ=255)
       ttAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGc                                                                                         <  1:2196163/100‑1 (MQ=255)
       ttAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGc                                                                                         >  1:2547193/1‑100 (MQ=255)
       ttAAAGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGc                                                                                         >  2:3385164/1‑100 (MQ=255)
       ntAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGc                                                                                         >  1:435507/2‑100 (MQ=255)
        tAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGcc                                                                                        <  2:866330/100‑1 (MQ=255)
          aCGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGa                                                                                      >  2:662417/1‑100 (MQ=255)
           cGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAg                                                                                     >  2:2915850/1‑100 (MQ=255)
           cGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAg                                                                                     <  2:3836283/100‑1 (MQ=255)
           cGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAg                                                                                     >  2:1040720/1‑100 (MQ=255)
           cGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAg                                                                                     >  2:3222065/1‑100 (MQ=255)
              tCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTg                                                                                  <  2:931081/100‑1 (MQ=255)
               cTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTgg                                                                                 <  2:2121868/100‑1 (MQ=255)
                tGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGt                                                                                <  1:2783358/100‑1 (MQ=255)
                 ggATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTg                                                                               >  2:3423520/1‑100 (MQ=255)
                      ncgcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGAt                                                                          >  1:3238818/2‑100 (MQ=255)
                       cgcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCCTATGGTCTGGCTTGAAGATGATCCAATCAGTGCCGAGCCGGTCCAGacg                                                                         >  2:3245316/1‑98 (MQ=255)
                        gctttGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATgc                                                                        <  1:3068888/97‑1 (MQ=255)
                             gAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAAATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgcttg                                                                   <  2:2794956/100‑1 (MQ=255)
                               aTCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgcttgct                                                                 >  2:464845/1‑100 (MQ=255)
                               aTCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgcttgct                                                                 >  2:1280767/1‑100 (MQ=255)
                                  atatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCttgcttgc                                                              <  2:3441457/100‑1 (MQ=255)
                                  atatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGAGCttgcttgc                                                              >  2:2920729/1‑100 (MQ=255)
                                   tatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcc                                                             <  2:725688/100‑1 (MQ=255)
                                   tatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcc                                                             >  2:708932/1‑100 (MQ=255)
                                   tatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcc                                                             <  1:373696/100‑1 (MQ=255)
                                   tatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcc                                                             <  1:3654112/100‑1 (MQ=255)
                                    atTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGccc                                                            >  1:2287744/1‑100 (MQ=255)
                                    atTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGccc                                                            <  1:1216798/100‑1 (MQ=255)
                                     ttttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcccc                                                           <  1:423723/100‑1 (MQ=255)
                                     ttttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcccc                                                           <  2:2578397/100‑1 (MQ=255)
                                     ttttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcccc                                                           <  2:2439875/100‑1 (MQ=255)
                                     ttttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGccac                                                           >  1:1215836/1‑98 (MQ=255)
                                       ttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAg                                                         <  2:3359845/100‑1 (MQ=255)
                                       ttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAg                                                         <  1:3172883/100‑1 (MQ=255)
                                         cccATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGaa                                                       >  2:1729531/1‑100 (MQ=255)
                                          ccATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAAc                                                      >  1:447768/1‑100 (MQ=255)
                                           cATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACa                                                     >  1:2628773/1‑100 (MQ=255)
                                            aTTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACaa                                                    >  1:680574/1‑100 (MQ=255)
                                             ttGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAAt                                                   <  1:1533774/100‑1 (MQ=255)
                                              tGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATg                                                  >  2:2350275/1‑100 (MQ=255)
                                               gCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATgg                                                 >  2:3190266/1‑100 (MQ=255)
                                                 cTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATgggg                                               <  1:1007429/100‑1 (MQ=255)
                                                    cTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCAt                                            <  1:793542/100‑1 (MQ=255)
                                                      ggCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGc                                          >  1:3584975/1‑100 (MQ=255)
                                                      ggCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGc                                          <  1:2735063/100‑1 (MQ=255)
                                                       gCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGcc                                         >  1:1119380/1‑100 (MQ=255)
                                                        cGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCg                                        <  1:2233506/100‑1 (MQ=255)
                                                        cGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCg                                        >  2:1965463/1‑100 (MQ=255)
                                                        cGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCg                                        <  1:471439/100‑1 (MQ=255)
                                                           cATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTg                                     <  2:1659736/100‑1 (MQ=255)
                                                             ttttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGcc                                   <  1:1944581/100‑1 (MQ=255)
                                                             ttttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGcc                                   <  1:3454614/100‑1 (MQ=255)
                                                             ttttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGcc                                   <  1:3462184/100‑1 (MQ=255)
                                                             ttttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGcc                                   <  1:734008/100‑1 (MQ=255)
                                                             ttttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGcc                                   <  2:1442839/100‑1 (MQ=255)
                                                              tttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTTCCCCAGAACAATTGGGCATGCCGCTGCCt                                  >  2:755089/1‑100 (MQ=255)
                                                                 ctcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCGGCTGCCTTTa                               >  2:2543249/1‑100 (MQ=255)
                                                                 ctcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTa                               <  2:1828752/100‑1 (MQ=255)
                                                                  tcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTaa                              >  2:3081490/1‑100 (MQ=255)
                                                                  tcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTaa                              >  1:2375428/1‑100 (MQ=255)
                                                                  tcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTaa                              >  2:111039/1‑100 (MQ=255)
                                                                   cAAACCGTATGGTCTGGTTTGAAGATGATCCAATCCGTGCCGAGCTGGTTCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAAc                             <  2:3942906/100‑1 (MQ=255)
                                                                   cAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAAc                             >  2:1618817/1‑100 (MQ=255)
                                                                   cAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAAc                             <  2:1893449/100‑1 (MQ=255)
                                                                    taaCCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACg                            <  1:1999580/99‑1 (MQ=255)
                                                                    agaCCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACg                            <  2:1799521/98‑1 (MQ=255)
                                                                    aaaCCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACg                            >  2:738106/1‑100 (MQ=255)
                                                                    aaaCCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACg                            <  1:2193735/100‑1 (MQ=255)
                                                                     aaCCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACga                           >  2:1968191/1‑99 (MQ=255)
                                                                      aCCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCa                          >  1:1756915/1‑100 (MQ=255)
                                                                         gTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAg                       <  1:1762655/100‑1 (MQ=255)
                                                                             ggTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAACGTTCt                   >  2:497622/1‑100 (MQ=255)
                                                                             ggTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCt                   <  2:2793451/100‑1 (MQ=255)
                                                                             ggTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCGGCCTGTAACGCAAAGTCCt                   >  2:1556973/1‑100 (MQ=255)
                                                                                  ggTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAAcc              >  1:1565072/1‑100 (MQ=255)
                                                                                  ggTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAAcc              >  2:1884992/1‑100 (MQ=255)
                                                                                      tGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGTCGCTGCCTTTAACGCAAAGTTCTAAAccagct          >  1:3211910/1‑96 (MQ=255)
                                                                                      tGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtc          <  2:1669323/100‑1 (MQ=255)
                                                                                       gAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtca         <  1:3286443/100‑1 (MQ=255)
                                                                                       gAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtca         <  2:1434519/100‑1 (MQ=255)
                                                                                        aaGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcat        <  2:88506/100‑1 (MQ=255)
                                                                                         aGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcatc       >  1:3694673/1‑100 (MQ=255)
                                                                                         aGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcatc       <  2:177014/100‑1 (MQ=255)
                                                                                         aGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcatc       <  2:2588653/100‑1 (MQ=255)
                                                                                            tgatCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAAt    >  2:2044677/1‑100 (MQ=255)
                                                                                            tgatCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAAt    >  1:991553/1‑100 (MQ=255)
                                                                                            tgatCCAATCAGTGCCGAGCGGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCGAAGTTCTAGGCCGGTCATCAAt    <  2:2815152/100‑1 (MQ=255)
                                                                                             gatCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATg   >  1:409748/1‑100 (MQ=255)
                                                                                              atCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGc  >  1:2002760/1‑100 (MQ=255)
                                                                                              atCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGc  >  1:3770899/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGC  >  NC_000913/3697757‑3697950

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑