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,8540AC97.8% 777.1 / 4.3 222Q466P (CAA→CCA) bcsEcellulose production protein
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (99/118);  total (102/120)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 223 total aligned reads displayed.

TGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGC  >  NC_000913/3697757‑3697950
                                                                                                 |                                                                                                
tGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCcaa                                                                                                <  2:1668567/100‑1 (MQ=255)
 gTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAAt                                                                                               <  1:1199198/100‑1 (MQ=255)
 gTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAAt                                                                                               >  2:5424780/1‑100 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  1:5498174/100‑1 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  1:1482684/100‑1 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  1:5133975/100‑1 (MQ=255)
  tCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATc                                                                                              <  2:3937635/100‑1 (MQ=255)
   cGGGTTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCa                                                                                             <  2:458624/100‑1 (MQ=255)
   cGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCa                                                                                             >  1:2581295/1‑100 (MQ=255)
    ggATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAg                                                                                            >  1:416710/1‑100 (MQ=255)
       ttAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGGCTGGTTTGAAGATGATCCAATCAGTGc                                                                                         >  2:1714900/1‑100 (MQ=255)
        tAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGcc                                                                                        <  1:1937185/100‑1 (MQ=255)
         aaCGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCg                                                                                       >  2:1870640/1‑100 (MQ=255)
         aaCGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCg                                                                                       >  2:3755501/1‑100 (MQ=255)
         aaCGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCg                                                                                       <  2:5098134/100‑1 (MQ=255)
         aaCGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCg                                                                                       >  2:2574985/1‑100 (MQ=255)
           cGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAg                                                                                     >  1:1247987/1‑100 (MQ=255)
            gATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGc                                                                                    <  2:3905634/100‑1 (MQ=255)
            gATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGc                                                                                    >  1:3675445/1‑100 (MQ=255)
             aTCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCt                                                                                   <  2:4286445/100‑1 (MQ=255)
             aTCTGGATACCGCGTTGAATCATACTTTCCCATCGCCTACTGGCGACATTTTCTCAAACCGCATGGCCTGGTTTGAAGATGATCCAATCAGTCCCGAGcg                                                                                   >  2:471872/1‑99 (MQ=255)
               cTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTgg                                                                                 <  1:5014980/100‑1 (MQ=255)
               cTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTgg                                                                                 <  1:4325342/100‑1 (MQ=255)
               cTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTgg                                                                                 <  1:3634397/100‑1 (MQ=255)
               cTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTgg                                                                                 <  1:3612532/100‑1 (MQ=255)
                tGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGt                                                                                <  1:2046079/100‑1 (MQ=255)
                tGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGt                                                                                <  2:4632358/100‑1 (MQ=255)
                gggATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGt                                                                                <  2:4327682/99‑1 (MQ=255)
                 ggATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTg                                                                               >  2:2928065/1‑100 (MQ=255)
                  gATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGc                                                                              <  2:3352120/100‑1 (MQ=255)
                   aTACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCa                                                                             <  1:1708644/100‑1 (MQ=255)
                    tACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAg                                                                            <  1:3888590/100‑1 (MQ=255)
                     aCCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGa                                                                           <  2:1183496/100‑1 (MQ=255)
                     aCCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGa                                                                           <  1:1073475/100‑1 (MQ=255)
                      ccGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGAt                                                                          >  2:4403321/1‑100 (MQ=255)
                      ccGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGAt                                                                          >  2:163747/1‑100 (MQ=255)
                       cgcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATg                                                                         >  2:2205818/1‑100 (MQ=255)
                       cgcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATg                                                                         <  1:2193988/100‑1 (MQ=255)
                       cgcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATg                                                                         <  2:4532240/100‑1 (MQ=255)
                        gcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATgc                                                                        >  1:1836153/1‑100 (MQ=255)
                        gcgTTGAATCATATTTTCCCATTGCCTACTGGCGACATCTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATgc                                                                        >  2:3297479/1‑100 (MQ=255)
                         cgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCGGATgcg                                                                       <  2:3655527/100‑1 (MQ=255)
                         cgTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATgcg                                                                       <  2:639684/100‑1 (MQ=255)
                          gatGAATAATATTTTCCCATCGACTACTGGCAACATTTTCTCAAACCGTAAGGTCTAATTTGTAGATCATTCAATCAGTGCCGTGCTGATGCAGATgcgc                                                                      >  2:303780/3‑100 (MQ=255)
                            tGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCGGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgctt                                                                    >  2:2688870/1‑100 (MQ=255)
                             gAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgcttg                                                                   >  1:3296804/1‑100 (MQ=255)
                               ntCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgcttgct                                                                 >  1:2341268/2‑100 (MQ=255)
                                tCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCgcttgctt                                                                >  1:4603649/1‑100 (MQ=255)
                                  atatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGcttgcttgn                                                              <  1:1808906/100‑2 (MQ=255)
                                  atatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCttgcttgc                                                              <  1:3274943/100‑1 (MQ=255)
                                   tatTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcc                                                             <  1:2129565/100‑1 (MQ=255)
                                     ttttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcccc                                                           <  2:3195553/100‑1 (MQ=255)
                                     ttttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGcccc                                                           <  1:1419025/100‑1 (MQ=255)
                                       ttCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAg                                                         >  2:3072856/1‑100 (MQ=255)
                                         cccATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGaa                                                       >  2:2461047/1‑100 (MQ=255)
                                          ccATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAAc                                                      >  2:2720335/1‑100 (MQ=255)
                                           cATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACa                                                     <  1:571074/100‑1 (MQ=255)
                                           cATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACa                                                     <  2:4758262/100‑1 (MQ=255)
                                               gCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATgg                                                 >  1:2665436/1‑100 (MQ=255)
                                               gCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATgg                                                 <  1:2240080/100‑1 (MQ=255)
                                               gCCTACTGGCGACATATTCCCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATgg                                                 <  1:5086824/100‑1 (MQ=255)
                                                 ccaCTGGCGAAATTTTCTCCAACCGTATGGTCTGGTTTGAAGCTGATCCAATCAGTGCCGAGCTGGTGCAGACGCGCCTGCTAGCCCCAGAACCATgggg                                               <  2:1363088/98‑1 (MQ=255)
                                                  tACTGTCGACATTTTCTCACACCGTACGGTCTGGTTTGAAGATGATCCAACCAGTGCCGATCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATggggn                                              <  1:4448149/100‑2 (MQ=255)
                                                  tACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGc                                              >  2:1020617/1‑100 (MQ=255)
                                                  tACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGc                                              <  2:1415647/100‑1 (MQ=255)
                                                  tACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGc                                              >  1:4157472/1‑100 (MQ=255)
                                                   aCTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCa                                             >  2:2282821/1‑100 (MQ=255)
                                                    cTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCAt                                            >  2:788438/1‑100 (MQ=255)
                                                    cTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCAt                                            <  1:2778598/100‑1 (MQ=255)
                                                       gCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGcc                                         >  1:4009393/1‑100 (MQ=255)
                                                        cGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCg                                        >  1:4406215/1‑100 (MQ=255)
                                                            aTTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGTTTGCCCCAGAACAATGGGGCATGCCGCTGc                                    >  1:811186/1‑100 (MQ=255)
                                                            aTTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGc                                    >  2:1249132/1‑100 (MQ=255)
                                                            aTTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGc                                    >  2:3922562/1‑100 (MQ=255)
                                                             ttttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGcc                                   <  2:2532722/100‑1 (MQ=255)
                                                               ttCTCAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCtt                                 >  2:1987278/1‑100 (MQ=255)
                                                                 ctcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTa                               >  2:5222917/1‑100 (MQ=255)
                                                                  tcAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTaa                              >  1:3376218/1‑100 (MQ=255)
                                                                   cAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAAc                             <  2:3586325/100‑1 (MQ=255)
                                                                   cAAACCGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAAc                             <  2:3863324/100‑1 (MQ=255)
                                                                       ccGTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCaa                         <  2:2102105/100‑1 (MQ=255)
                                                                         gTATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAg                       <  2:2707029/100‑1 (MQ=255)
                                                                          tATGGTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGt                      >  2:4064436/1‑100 (MQ=255)
                                                                            tGGTCTGGTTTGAATATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTCGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTc                    <  2:5023596/100‑1 (MQ=255)
                                                                             ggTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCt                   <  1:4748021/100‑1 (MQ=255)
                                                                             ggTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCt                   >  1:430881/1‑100 (MQ=255)
                                                                              gTCTGGTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTa                  >  1:2708840/1‑100 (MQ=255)
                                                                                  ggTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAAcc              >  1:2452915/1‑100 (MQ=255)
                                                                                  ggTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAAcc              <  1:3241555/100‑1 (MQ=255)
                                                                                   gTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCg             <  2:3265946/100‑1 (MQ=255)
                                                                                   gTTTGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCACTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCg             >  1:3416062/1‑100 (MQ=255)
                                                                                    tttGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCgg            <  1:2485679/100‑1 (MQ=255)
                                                                                      tGAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtc          >  1:1629077/1‑100 (MQ=255)
                                                                                       gAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCTCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtca         <  2:487908/100‑1 (MQ=255)
                                                                                       gAAGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtca         <  1:3597081/100‑1 (MQ=255)
                                                                                         aGATGATCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcatc       >  1:4967852/1‑100 (MQ=255)
                                                                                          gatgctCCAATCAGTGCCGAGCTGGTGCAGCTGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCTCTGCCTTTAACGCAAAGTTCTAAACCGGtcatca      <  1:3058214/100‑1 (MQ=255)
                                                                                            tgatCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGtcagatcg    >  2:3543164/1‑95 (MQ=255)
                                                                                             gatCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATg   >  2:3423182/1‑100 (MQ=255)
                                                                                              atCCAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGc  <  2:2115609/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAACCGTATGGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGC  >  NC_000913/3697757‑3697950

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