Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,692,114 C→T S165N (AGC→AAC)  bcsB ← cellulose synthase periplasmic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,692,1140CT97.2% 233.6 / ‑1.1 72S165N (AGC→AAC) bcsBcellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (38/31);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-01

AGGCTTGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCG  >  NC_000913/3691869‑3692360
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
agGCTTGTTGCAACCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAAGGGCACCGGGAAGGGTGACAGGTAATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATACAGCCAAGACGTGGTGCTGGCCGTATTTTCGCACACGTCCTGTTAATGTCCGACAAACTCCACCCGTACACGCTTGTAGTCGTtaat                                                                                                                                                                                                                                                    >  2:10205/1‑250 (MQ=255)
 ggCTTGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTtaata                                                                                                                                                                                                                                                   >  1:4953/1‑250 (MQ=255)
     tgttgCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACAGTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCTGCCAAAGCGTTGTGCGGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTTAAGTCGTTAATAAACa                                                                                                                                                                                                                                               <  2:60517/250‑1 (MQ=255)
        tgCAGCCCAACAACCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGGCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACAGCCAGCCAAAGCGGGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACGCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTg                                                                                                                                                                                                                                            <  2:204172/250‑1 (MQ=255)
            gCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTa                                                                                                                                                                                                                                       <  1:245120/251‑1 (MQ=255)
               cAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGAGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATg                                                                                                                                                                                                                                    >  1:542329/1‑251 (MQ=255)
                 aCATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATggg                                                                                                                                                                                                                                  >  1:399793/1‑251 (MQ=255)
                  cATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGAGCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATggg                                                                                                                                                                                                                                  <  1:174624/250‑1 (MQ=255)
                  cATCCGGCGCACCCGAAAAGACCATCGGCAAGGTGTTGGTGCGGTTAGCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGCTAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATACTGGCCGACAAACTCAAGCCGTACACGGTTGAAGTCGTTAATAACAAGTGGGTTAATgggg                                                                                                                                                                                                                                 >  2:174624/1‑250 (MQ=255)
                     ccGGCGCTCCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCTTTCTTCACATTCAGGGTGTGATAGGTCAGATCCAGTCCTCTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGCTTATCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACTGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAtt                                                                                                                                                                                                                              <  1:306355/251‑1 (MQ=255)
                      cGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAttt                                                                                                                                                                                                                             <  1:268898/251‑1 (MQ=255)
                             cccGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAg                                                                                                                                                                                                                      <  1:505720/251‑1 (MQ=255)
                             cccGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGAGCCAGGCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTGTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAg                                                                                                                                                                                                                      <  2:542329/251‑1 (MQ=255)
                                gCAAAGACCATGGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAACGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGCGCTGGCCGGTTATTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGGACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTAGCGCCAGCGt                                                                                                                                                                                                                   >  2:411932/1‑251 (MQ=255)
                                  aaaGACCATCGGCAAGGTGTTGGTTCGGTTATCGCCCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttt                                                                                                                                                                                                                 <  2:144945/251‑1 (MQ=255)
                                    aGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                               <  1:402012/251‑1 (MQ=255)
                                    aGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                               <  1:118328/251‑1 (MQ=255)
                                         aTCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCCTTTGCGCCAGCGTTTTTTTAcc                                                                                                                                                                                                          >  2:28805/1‑251 (MQ=255)
                                            ggCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCa                                                                                                                                                                                                        >  1:532204/1‑250 (MQ=255)
                                              cAAGGTGTTGGTGCGGCCATCGCCCGGGTCAAAGAGTGGCACCGGGACGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATACACAGTGGGTTCATGGGCATTTGCGCCAGCGTTTTTTTACCCAACt                                                                                                                                                                                                     >  2:525978/1‑251 (MQ=255)
                                                    gtTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTTGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGGTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCt                                                                                                                                                                                                <  1:441502/250‑1 (MQ=255)
                                                         tGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTc                                                                                                                                                                                          >  1:268626/1‑251 (MQ=255)
                                                          gCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAATTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAAGGGCCGACAAACGCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATGTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTc                                                                                                                                                                                          <  1:236632/250‑1 (MQ=255)
                                                             gTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTg                                                                                                                                                                                      >  2:393899/1‑251 (MQ=255)
                                                             gTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTg                                                                                                                                                                                      <  1:33843/251‑1 (MQ=255)
                                                                 tCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGGTGAAGTCGGTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGGCACTGgc                                                                                                                                                                                    >  2:410539/1‑249 (MQ=255)
                                                                   gCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAAAAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgcagc                                                                                                                                                                                 >  1:190346/1‑250 (MQ=255)
                                                                       gggTCAAAGAATGGCAACGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGcc                                                                                                                                                                            >  2:519698/1‑251 (MQ=255)
                                                                        ggTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGccc                                                                                                                                                                           >  1:203042/1‑251 (MQ=255)
                                                                        ggTCAAAGAATGGCACCGGGAAGGGTGACAGGTCATTCTTCACATTCAGGGTCTGCTAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGCTACTGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTtaataa                                                                                                                                                                                                                                                  >  1:496912/1‑180 (MQ=255)
                                                                              aGAATGGCACCGGGAAGTGTGCCAGGTCATTCTTCACTTTCAGGGTCTGATAGGTCAGATCCAGTCCCCTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGt                                                                                                                                                                     <  2:4953/251‑1 (MQ=255)
                                                                                 aTGGCACCGGGAAGTAAGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGCTCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGAGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTtcc                                                                                                                                                                  >  2:305983/1‑250 (MQ=255)
                                                                                          ggAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGAGCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTgaga                                                                                                                                                         <  1:80903/251‑1 (MQ=255)
                                                                                               tgtgACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATACACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAAc                                                                                                                                                    >  2:193965/1‑251 (MQ=255)
                                                                                               tgtgACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAAc                                                                                                                                                    >  2:112062/1‑251 (MQ=255)
                                                                                                 tgACAGGACATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACCGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCt                                                                                                                                                  >  1:143527/1‑251 (MQ=255)
                                                                                                     aGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTa                                                                                                                                               >  1:168574/1‑250 (MQ=255)
                                                                                                           ttctacACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGGTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCg                                                                                                                                         >  2:464126/1‑250 (MQ=255)
                                                                                                                cacaTTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGTCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGGTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCAGCAGTTCATCATTGAGATAAACCTTTAACTGCGACTgg                                                                                                                                    <  1:25019/250‑1 (MQ=255)
                                                                                                                     tCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgc                                                                                                                                                                                                                           >  2:370959/1‑158 (MQ=255)
                                                                                                                     tCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgc                                                                                                                                                                                                                           <  1:370959/158‑1 (MQ=255)
                                                                                                                      cAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTAAACTGCGAATGGACAgga                                                                                                                              >  1:400591/1‑249 (MQ=255)
                                                                                                                         ggTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTa                                                                                                                          <  1:112062/251‑1 (MQ=255)
                                                                                                                          gTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTa                                                                                                                          >  2:395022/1‑250 (MQ=255)
                                                                                                                             tGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAAcg                                                                                                                       >  1:454891/1‑250 (MQ=255)
                                                                                                                                     cAGATCCAGTCCACTGCGCCGCCCAACATCCAGCCAAAGCGGGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAAcgatggcgat                                                                                                               <  1:525978/250‑1 (MQ=255)
                                                                                                                                           cAGTCCACAGCTCCGCCCAACATCCAGCCAAAGCGTGGAGCTGGCCGGTTTTTCGCACACGTCATGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTACTAAACAGAGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCACATCCGTTCATCACTGAGATAAACCTTTAACTGCTAATGGACAGGCATTAACGATGGCGATGGGgtgt                                                                                                        >  2:316209/1‑251 (MQ=255)
                                                                                                                                               ccACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTc                                                                                                    <  1:254167/251‑1 (MQ=255)
                                                                                                                                                         gCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGc                                                                                          >  2:290016/1‑251 (MQ=255)
                                                                                                                                                                    aGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGt                                                                               >  1:137239/1‑251 (MQ=255)
                                                                                                                                                                      ccAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGCCAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTc                                                                              >  2:396902/1‑250 (MQ=255)
                                                                                                                                                                        aaaGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACGGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGCGCATCGCCTTCGTcacc                                                                           >  2:277001/1‑251 (MQ=255)
                                                                                                                                                                                gtgCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCAGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTACTGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATc                                                                   >  2:226134/1‑251 (MQ=255)
                                                                                                                                                                                        cGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTgcgc                                                            >  1:104071/1‑250 (MQ=255)
                                                                                                                                                                                         ggTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCAt                                                          >  1:91176/1‑251 (MQ=255)
                                                                                                                                                                                          gTTTGTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGTCGATTGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATa                                                         <  2:238937/251‑1 (MQ=255)
                                                                                                                                                                                             tttCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATAcc                                                       >  2:143180/1‑250 (MQ=255)
                                                                                                                                                                                                     acGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGGGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGCGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATCCCAAACTCAAt                                              >  2:182900/1‑251 (MQ=255)
                                                                                                                                                                                                            gATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCg                                       <  1:193965/251‑1 (MQ=255)
                                                                                                                                                                                                                aTGGCCGACAATCTCCATCCGTACACGGTTGAAGTCGTTAATAATCAGTGTGTTAATGGGCATTTGCGCCAGCGTTTTTCTACCCATCTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCTGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCg                                   <  1:336321/251‑1 (MQ=255)
                                                                                                                                                                                                                aTGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTTAGATaa                                                                                                                                                      <  2:236009/136‑1 (MQ=255)
                                                                                                                                                                                                                aTGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATaa                                                                                                                                                      >  1:236009/1‑136 (MQ=255)
                                                                                                                                                                                                                aTGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCg                                   <  1:264165/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggCCGACAAACTCCAGCCGTACACGGTTGGAGTCGTTAATAAACGGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCATCACGCCCATCGGTTCATCATTGAGATGAACCTTTAACTGCGACTGGACAGGCAGTAAcga                                                                                                                      >  1:81545/1‑166 (MQ=255)
                                                                                                                                                                                                                           aCTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCAGCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCa                        <  2:3458/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTa                    <  1:95483/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTCTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGATATAAACCTTTAACGGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGTTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAATCTCAATGCTGCCGTTCGGGTTAATGCCACGTAGc                  <  1:156526/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ccacGGTTGAATTCGTTAATAAACAGTGGTGTAATTGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACGGGACAGGCAGTAACGATGGCGATGGGGTGTAGTCGAGGTTGAGCATCGCCTTCGGCACCACTTCTTCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTATTGCCACGTAGCTCCAt             <  2:439438/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cacGGTTGAAGTCGTTAATAAACAGGGGGGTAATGGGCATTTGCGCCCGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCCAACTCTATGCTGCCTGTCGGGTTAATGCCACGTTGCACCAt             >  2:126018/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcg     <  1:79028/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCCGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATAGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATCCCAAACTCAATGCTGCCGTTCGGGGTAATGCCACGTAGCACCATGATGCCcggcg  >  2:80692/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcggcg  <  2:399793/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
AGGCTTGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCG  >  NC_000913/3691869‑3692360

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: