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

TGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTG  >  NC_000913/3691874‑3692353
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
tggtgCAGACCAACATCCTGCGCACCCGCAAAGTCCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGGGACAGGTCAT‑TCGGCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCTGTTTTTTCGCACACGTCCTGATATTGGCCGACAAACTCCAGGCGTACACGTTTGAAGTCGTTATTAAAc                                                                                                                                                                                                                                         <  2:276447/246‑1 (MQ=255)
 gttgCAGCCCAACATCCGGCGCACCCGCAACGACCATCGGCAAGGTGTTGGTGCGGTCATCGCTCGGGTCAAAGCATGGCACCGGGAAGTGTGCCAGGCCAT‑TCTTCACATTCAGGGTCTGATAGGCCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGCGCTGGCCGGTTTTTCGCACACGTCCTGATACTGGCCGACAACCTCCAGCCGTACACCGCTGAAGCCGTTAATAAAccgt                                                                                                                                                                                                                                      >  2:137527/1‑248 (MQ=255)
   tgCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTg                                                                                                                                                                                                                                     >  1:180090/1‑250 (MQ=255)
   tgCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGGGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTg                                                                                                                                                                                                                                     <  2:180090/250‑1 (MQ=255)
           aaCATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATgg                                                                                                                                                                                                                            <  1:158144/251‑1 (MQ=255)
             cATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGc                                                                                                                                                                                                                          >  1:310492/1‑251 (MQ=255)
                 cGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAttt                                                                                                                                                                                                                      >  2:340219/1‑251 (MQ=255)
                    cgcACCCGCAAAGACTAGCGGCAAGGTGTTGGGGCGGTTATCGCACGGGTCAAAGAATGGCGCCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTTGTGCTGGCCGGTTTTTCGCACACGTCCGTATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgcg                                                                                                                                                                                                                   <  2:371791/251‑1 (MQ=255)
                          cGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGTCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCAAACGACCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGTGCATTTGCGCCAGCg                                                                                                                                                                                                             >  2:370364/1‑251 (MQ=255)
                               aGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                        >  1:166621/1‑251 (MQ=255)
                               aGACCAACGGCAAGGTGTTGGTGCGGCTATCGCTAGGGGCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGCTACAGTCGACTGCCCCGCCCACCATCCAGCCAAAGCGTGGTGCTGCCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACCTGTTTGAAGTCGTTCATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                        >  2:382002/1‑251 (MQ=255)
                                gACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGACAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACAGTCAGGGTCTGATAGGTCAGATCCAGTACACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTTATAATGGCCGACAAACTCCAGCCGTACACGTTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGtttttt                                                                                                                                                                                                       <  1:165188/251‑1 (MQ=255)
                                      cGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGCAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGACAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGCCAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAccc                                                                                                                                                                                                  >  1:408424/1‑250 (MQ=255)
                                      cGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAccc                                                                                                                                                                                                  <  1:371319/250‑1 (MQ=255)
                                      cGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTCCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACATTGGTTTAATGGGCATGTGCGCCAGCGTTTTTTTACCCa                                                                                                                                                                                                 <  1:379680/251‑1 (MQ=255)
                                       ggCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCa                                                                                                                                                                                                 >  2:301036/1‑250 (MQ=255)
                                               gttgttGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGGACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCTGCGTTTTTTTACCCATCTGTTCt                                                                                                                                                                                         <  2:44647/250‑1 (MQ=255)
                                                           aTCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgc                                                                                                                                                                             >  1:377716/1‑250 (MQ=255)
                                                            tCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgc                                                                                                                                                                             >  1:335309/1‑249 (MQ=255)
                                                                         aGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCAACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGATGAAGCCGTTAATAAACAGTGGGTTAATGGGCTTTTGCGCCAGCGTTTTTTTACCCAACTGATCTTTGGTCACTGGCAGCACGCCCATCAGt                                                                                                                                                              >  2:228139/1‑251 (MQ=255)
                                                                              ggCACCGGTAAGGGTGTCAGTTCCT‑GCTTCATATTCAGGGTCTGAAATGTAAGATCCAGTCCACGTCCCCGCCCAACATCCATCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTGCTGATAATGGCCGACAAACTCCAGCCTTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTGCCCACCCGTTCTTTGGTCACTGGCAGCTCGCCCATCAGTccatc                                                                                                                                                         <  2:141627/251‑1 (MQ=255)
                                                                                    gggAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTgag                                                                                                                                                   <  1:288463/251‑1 (MQ=255)
                                                                                        aGTGTGACAGGTCATGTCT‑CACATTCGGGGTCGGACAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGGGGTGCTGGCTGGTTTTTGGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAGTTGCGCCAGCGTTTTTTTACCCAACTGGTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATaa                                                                                                                                               <  2:141906/251‑1 (MQ=255)
                                                                                        aGTGTGACAGGTCAT‑TCTGCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCACCATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATaa                                                                                                                                               <  2:318500/251‑1 (MQ=255)
                                                                                             gACAGGTCAT‑TCTCCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTGTTTTACCCAACGGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCtt                                                                                                                                          <  2:146835/251‑1 (MQ=255)
                                                                                             gACAGGTCAG‑TCTTCACATGCAGGGTCGTATAGGTCACAGCCAGTCCGCTGCTCCGCGCAACATCCAGACCAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGGGGGTTAATGGGCATTTGCTCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAAccgt                                                                                                                                          <  2:2322/251‑3 (MQ=255)
                                                                                               cAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTa                                                                                                                                        <  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)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
TGTTGCAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCAT‑TCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTG  >  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: ‑