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,1140CT98.5% 232.6 / NA 69S165N (AGC→AAC) bcsBcellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/34);  minor base G (0/1);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTG  >  NC_000913/3691874‑3692353
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
tggtgCAGACCAACATCCTGCGCACCCGCAAAGTCCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGGGACAGGTCATTCGGCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCTGTTTTTTCGCACACGTCCTGATATTGGCCGACAAACTCCAGGCGTACACGTTTGAAGTCGTTATTAAAc                                                                                                                                                                                                                                         <  2:276447/246‑1 (MQ=255)
 gttgCAGCCCAACATCCGGCGCACCCGCAACGACCATCGGCAAGGTGTTGGTGCGGTCATCGCTCGGGTCAAAGCATGGCACCGGGAAGTGTGCCAGGCCATTCTTCACATTCAGGGTCTGATAGGCCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGCGCTGGCCGGTTTTTCGCACACGTCCTGATACTGGCCGACAACCTCCAGCCGTACACCGCTGAAGCCGTTAATAAAccgt                                                                                                                                                                                                                                      >  2:137527/1‑248 (MQ=255)
   tgCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTg                                                                                                                                                                                                                                     >  1:180090/1‑250 (MQ=255)
   tgCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGGGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTg                                                                                                                                                                                                                                     <  2:180090/250‑1 (MQ=255)
           aaCATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATgg                                                                                                                                                                                                                            <  1:158144/251‑1 (MQ=255)
             cATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGc                                                                                                                                                                                                                          >  1:310492/1‑251 (MQ=255)
                 cGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAttt                                                                                                                                                                                                                      >  2:340219/1‑251 (MQ=255)
                    cgcACCCGCAAAGACTAGCGGCAAGGTGTTGGGGCGGTTATCGCACGGGTCAAAGAATGGCGCCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTTGTGCTGGCCGGTTTTTCGCACACGTCCGTATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgcg                                                                                                                                                                                                                   <  2:371791/251‑1 (MQ=255)
                          cGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGTCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCAAACGACCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGTGCATTTGCGCCAGCg                                                                                                                                                                                                             >  2:370364/1‑251 (MQ=255)
                               aGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                        >  1:166621/1‑251 (MQ=255)
                               aGACCAACGGCAAGGTGTTGGTGCGGCTATCGCTAGGGGCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGCTACAGTCGACTGCCCCGCCCACCATCCAGCCAAAGCGTGGTGCTGCCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACCTGTTTGAAGTCGTTCATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                        >  2:382002/1‑251 (MQ=255)
                                gACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGACAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACAGTCAGGGTCTGATAGGTCAGATCCAGTACACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTTATAATGGCCGACAAACTCCAGCCGTACACGTTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGtttttt                                                                                                                                                                                                       <  1:165188/251‑1 (MQ=255)
                                      cGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGCAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGACAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGCCAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAccc                                                                                                                                                                                                  >  1:408424/1‑250 (MQ=255)
                                      cGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAccc                                                                                                                                                                                                  <  1:371319/250‑1 (MQ=255)
                                      cGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTCCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACATTGGTTTAATGGGCATGTGCGCCAGCGTTTTTTTACCCa                                                                                                                                                                                                 <  1:379680/251‑1 (MQ=255)
                                       ggCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCa                                                                                                                                                                                                 >  2:301036/1‑250 (MQ=255)
                                               gttgttGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGGACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCTGCGTTTTTTTACCCATCTGTTCt                                                                                                                                                                                         <  2:44647/250‑1 (MQ=255)
                                                           aTCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgc                                                                                                                                                                             >  1:377716/1‑250 (MQ=255)
                                                            tCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgc                                                                                                                                                                             >  1:335309/1‑249 (MQ=255)
                                                                         aGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCAACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGATGAAGCCGTTAATAAACAGTGGGTTAATGGGCTTTTGCGCCAGCGTTTTTTTACCCAACTGATCTTTGGTCACTGGCAGCACGCCCATCAGt                                                                                                                                                              >  2:228139/1‑251 (MQ=255)
                                                                              ggCACCGGTAAGGGTGTCAGTTCCTGCTTCATATTCAGGGTCTGAAATGTAAGATCCAGTCCACGTCCCCGCCCAACATCCATCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTGCTGATAATGGCCGACAAACTCCAGCCTTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTGCCCACCCGTTCTTTGGTCACTGGCAGCTCGCCCATCAGTccatc                                                                                                                                                         <  2:141627/251‑1 (MQ=255)
                                                                                    gggAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTgag                                                                                                                                                   <  1:288463/251‑1 (MQ=255)
                                                                                        aGTGTGACAGGTCATGTCTCACATTCGGGGTCGGACAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGGGGTGCTGGCTGGTTTTTGGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAGTTGCGCCAGCGTTTTTTTACCCAACTGGTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATaa                                                                                                                                               <  2:141906/251‑1 (MQ=255)
                                                                                        aGTGTGACAGGTCATTCTGCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCACCATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATaa                                                                                                                                               <  2:318500/251‑1 (MQ=255)
                                                                                             gACAGGTCATTCTCCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTGTTTTACCCAACGGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCtt                                                                                                                                          <  2:146835/251‑1 (MQ=255)
                                                                                             gACAGGTCAGTCTTCACATGCAGGGTCGTATAGGTCACAGCCAGTCCGCTGCTCCGCGCAACATCCAGACCAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGGGGGTTAATGGGCATTTGCTCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAAccgt                                                                                                                                          <  2:2322/251‑3 (MQ=255)
                                                                                               cAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTa                                                                                                                                        <  1:86744/251‑1 (MQ=255)
                                                                                                            cacaTTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATACACCTTTAACTGCGACTgg                                                                                                                             >  1:221718/1‑250 (MQ=255)
                                                                                                               aTTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCACAACATCCAGCCAAAGAGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGCAAAACTCCAGCCGTACACGGTTGAAGTCGATAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGAc                                                                                                                           >  1:358886/1‑249 (MQ=255)
                                                                                                                     ggTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTa                                                                                                                   <  1:345084/251‑1 (MQ=255)
                                                                                                                     ggTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTGTAACTGCGACTGGACAGGCAGTa                                                                                                                   <  1:380272/251‑1 (MQ=255)
                                                                                                                                 cAGATCCCGTCCTCTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGTTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAACCCTTTAACTGCGACTGGACAGGCAGTAAcgatggcgat                                                                                                        <  2:79451/250‑1 (MQ=255)
                                                                                                                                 cAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGACGGGTTTTCGCACACGTCCTGATAATGGCCGCCAAACTCCAGCCGTACACGGTTGAGGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTTAGATACACCTTTAACTGCGCCTGGACCGGCCGTAActatggcgatg                                                                                                       >  2:255073/1‑251 (MQ=255)
                                                                                                                                        agaccACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCAGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCCTTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCCCTGGCAGCACGCCCACCAGTTCATCATTGAGATACACCTTTAACTGCGACTGGACAGCCCGTAACGACGGCGATGGGgtgt                                                                                                 >  2:362024/4‑250 (MQ=255)
                                                                                                                                         gTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTAt                                                                                               <  1:360373/251‑1 (MQ=255)
                                                                                                                                          tCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTAtc                                                                                              >  1:436847/1‑250 (MQ=255)
                                                                                                                                           ccACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCa                                                                                                                      <  1:399946/226‑1 (MQ=255)
                                                                                                                                           ccACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCa                                                                                                                      >  2:399946/1‑226 (MQ=255)
                                                                                                                                           ccACTGCTCCGCCAAACATCCAGCCACAGCGTGGTGCAGTCCGCTTTTTGGCACACGTCCTGATCATGGCCGACAAACTCCAGCCGTACACGGGTGCAGTCGTTAATAAACAGAGGGTTAAAGGGCATTTGCGCCAGCGTTTTTTTACCCAACTATTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCTTTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACCATTGCGATGGGTAGTATTc                                                                                             >  2:82983/1‑251 (MQ=255)
                                                                                                                                                   ccGCCCATCATCCAGCCAAAGCGTGTGGCGGCCCGGTTTGTCGCACCCGTACTAAGAATGGCCGACAACCTCCAGCCTGACACGGTAGAATTCGGGAATGACCAGTGGGTTATTGGGCATTTGGGCCAGCGTTTTTTTGCCCAACTTTTCGTTGGTCACTTGCATCACGCCCATCAGTTCGTCATTGAGATAAACCTTTAACTGCGACTGGACAGGCATTAACGATGGCGATGGGGTGTATTCGAGGt                                                                                        <  1:23791/248‑1 (MQ=255)
                                                                                                                                                   ccGCCCAGCATCCAGCCAAAGCGTGGGGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACGGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTg                                                                                      <  1:394698/250‑1 (MQ=255)
                                                                                                                                                         aaCATCCAGCCGAAGGGTGTTGCTGGCCGGTTTTTCGCACACGTCCTGAGAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACTCCCATCAGTTCATCATTGAGAGAAGCCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCTAGGGTGAGCATCg                                                                               <  2:245399/251‑1 (MQ=255)
                                                                                                                                                             tCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCtt                                                                           >  2:188679/1‑251 (MQ=255)
                                                                                                                                                               cAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTGTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCg                                                                         <  1:249246/251‑1 (MQ=255)
                                                                                                                                                                   cAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTATTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCCTCGCCTTCGTca                                                                      >  2:345773/1‑250 (MQ=255)
                                                                                                                                                                    aaaGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGGCAAACTCAAGCCGTACACAGTTGAAGTCGTTAATACACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACACAACTGCACTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGCCAGTAACGATGGCGATGGGGTGAATTCGAGGTTGAGCATCGCCTTCGTcacc                                                                    >  1:392501/1‑251 (MQ=255)
                                                                                                                                                                     aaGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAAc                                                                                                                                      <  1:20131/184‑1 (MQ=255)
                                                                                                                                                                     aaGCGTGGTGCTGGCCGGTTTTTCGCACACGCCCTGCTAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGACGTTAATAAACAGTGGGTTAATGGGCATTTGAGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAACCCTTTAAc                                                                                                                                      >  2:20131/1‑184 (MQ=255)
                                                                                                                                                                            gtgCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTGCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACGTCATc                                                            <  2:212501/251‑1 (MQ=255)
                                                                                                                                                                              gCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCg                                                           <  1:105824/250‑1 (MQ=255)
                                                                                                                                                                                    cGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCa                                                    <  1:217288/251‑1 (MQ=255)
                                                                                                                                                                                     ggTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGCGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCAt                                                   >  2:297337/1‑251 (MQ=255)
                                                                                                                                                                                         tttCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTcatcat                                                                                                                                                       >  2:167290/1‑147 (MQ=255)
                                                                                                                                                                                         tttCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTcatcat                                                                                                                                                       <  1:167290/147‑1 (MQ=255)
                                                                                                                                                                                              cacaCGTCCTGATAATGGCCGACAAACTCCAGCCGTAAACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGATGCGCATACCAAACTc                                          >  1:401385/1‑251 (MQ=255)
                                                                                                                                                                                                cacGTCCAGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGGCGTTAATAAACAGTGGGTTAATGAGCATTTGCGCCAGCGTTTTTTTACCCAACTGGTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGGGTATTCGAGGTTGAGCATCGCCTTCGGCACCACTTCATCGCTGCGCATACCAAACTCaa                                        <  2:364981/251‑1 (MQ=255)
                                                                                                                                                                                                         agaatgacagACAAACTCAAGCCGTACACGGTTGAAGTCGTTACTAAACAGTGGATTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCTTCCACACTTCATCGCTGCGCATACCAAACTCAAACCTGCCGt                               >  2:358689/10‑251 (MQ=255)
                                                                                                                                                                                                           aaTGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTc                             <  1:301036/251‑1 (MQ=255)
                                                                                                                                                                                                                  gACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTaa                      <  2:310492/251‑1 (MQ=255)
                                                                                                                                                                                                                    cAAACTCCAGCCGTACACGGTTGAAGTCGTTACTAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATg                    >  2:385416/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAACCTCAATGCTGCCGTTCGGGTTAATg                    >  1:311190/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATg                    >  1:244368/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaCTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATTGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCAGCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGcc                  <  2:53881/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCa          <  2:166621/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCACTGCTGCCGTTCGGGTTAATGCCACGTAGCAcc        >  1:138573/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGGTGCGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAACCTCAATTCTGCCGTTCGGGTTGATGCCACGTAACAcc        >  2:294060/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTACACGGTTGAAGTCGTTAATAAACAGTGGGTAAATGGACATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGAAAGACAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTCGCAcc        >  1:171012/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCAtgctg  <  1:297337/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGACAGCGTTTTTTTACCGAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATAATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGCTGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGACACGTAGCACCAtgctg  >  1:286245/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
TGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTG  >  NC_000913/3691874‑3692353

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