Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,680,910 G→T 66.1% H542N (CAT→AAT)  bcsC ← cellulose synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,680,9100GT66.1% 383.0 / inf 378H542N (CAT→AAT) bcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base G (59/69);  new base T (117/131);  total (178/200)
Fisher's exact test for biased strand distribution p-value = 9.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 378 total aligned reads displayed.

TCTGCCTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGC  >  CP009273/3680762‑3681056
                                                                                                                                                    |                                                                                                                                                  
 cTGCCTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTAc                                                                                                                                                  <  2:4216947/150‑1 (MQ=255)
 cTGCCTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTAc                                                                                                                                                  <  2:1761950/150‑1 (MQ=255)
 cTGCCTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTAc                                                                                                                                                  <  2:333200/150‑1 (MQ=255)
 cTGCCTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGAc                                                                                                                                                  >  1:1351697/1‑150 (MQ=255)
    ccTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTTCTGTTACACTGCGCACGCGGCCGGCGCTTGATATGCGCCAGCGCCGCGCTGTCCTGGTCATTACCag                                                                                                                                               >  1:43556/1‑150 (MQ=255)
     ctctTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCaga                                                                                                                                              <  1:4313457/150‑1 (MQ=255)
     ctctTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCaga                                                                                                                                              <  1:4780434/150‑1 (MQ=255)
       ctTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCagaga                                                                                                                                            >  2:290477/1‑150 (MQ=255)
        tttGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCagagag                                                                                                                                           >  2:4690427/1‑150 (MQ=255)
        tttGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCagagag                                                                                                                                           >  2:582967/1‑150 (MQ=255)
          tGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGt                                                                                                                                         <  2:1973458/150‑1 (MQ=255)
           gCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTa                                                                                                                                        <  2:339850/150‑1 (MQ=255)
           gCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTa                                                                                                                                        <  1:4458283/150‑1 (MQ=255)
           gCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTa                                                                                                                                        <  1:3078145/150‑1 (MQ=255)
              gCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAg                                                                                                                                     >  2:229572/1‑150 (MQ=255)
               cTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGc                                                                                                                                    <  2:4737548/150‑1 (MQ=255)
                  tctcGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTACTGAATAATTCTGTTACACAGCGCACGCCGCAGGCTAATGATATGCGCCAGCGCCGCTCTGTCCAGGTCATAAACACAGAGGTCCAGcccc                                                                                                                                 >  1:2868625/1‑149 (MQ=255)
                   ctcGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGt                                                                                                                                <  2:1910207/150‑1 (MQ=255)
                    tcGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTa                                                                                                                               <  2:1376106/150‑1 (MQ=255)
                    tcGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTa                                                                                                                               <  2:3834938/150‑1 (MQ=255)
                    tcGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTa                                                                                                                               <  1:4538924/150‑1 (MQ=255)
                      gCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAg                                                                                                                             <  2:1786148/150‑1 (MQ=255)
                       cAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGc                                                                                                                            <  1:391588/150‑1 (MQ=255)
                       cAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAAACAGCGCCTGAATATTTCTGTTCCCCTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGc                                                                                                                            >  1:1248000/1‑150 (MQ=255)
                           cGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTaa                                                                                                                        >  1:1533616/1‑150 (MQ=255)
                           cGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTaa                                                                                                                        >  2:2770743/1‑150 (MQ=255)
                              ttAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAAcc                                                                                                                     <  2:1618048/150‑1 (MQ=255)
                               tAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCt                                                                                                                    >  2:393277/1‑150 (MQ=255)
                               tAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCt                                                                                                                    <  2:1824750/150‑1 (MQ=255)
                                aGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTg                                                                                                                   <  1:2620627/150‑1 (MQ=255)
                                   ggTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTc                                                                                                                <  2:1569994/150‑1 (MQ=255)
                                       tCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCggg                                                                                                            >  1:637429/1‑150 (MQ=255)
                                         cAGCACCTGAACGCTTAGCAATCGATTAACCACAACCTTAATAATGCAGTTCCAATGCGCACGCGGCAAGCAATAGAAATGCGCAAGCGCCGCTCTGTCCAGGACCTAACCACAGAGGGAAAGCCCGTAAGCGTAAACCTGCTCAGGGTc                                                                                                          >  2:1177751/1‑150 (MQ=255)
                                             aCCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTc                                                                                                      <  2:26751/150‑1 (MQ=255)
                                                tGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACATCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGAc                                                                                                   >  1:802222/1‑150 (MQ=255)
                                                tGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGGAAGCGTAAACCTGCTCCGGGGCGTTCGAc                                                                                                   >  2:4393599/1‑150 (MQ=255)
                                                 gATCGCTTTGCAGTCGATTAAACAGCGCCTGAATATTTCTGTTCCACAGCGCACGCGGCAGGCGATTGATATGCGCCACCGCCGCGCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAACGTAAACCAGCTCCGGGTCGtggggat                                                                                                  >  1:458969/1‑144 (MQ=255)
                                                   tCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTc                                                                                                <  1:1990767/150‑1 (MQ=255)
                                                      cTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTAGACTTatgc                                                                                             >  2:3778762/1‑150 (MQ=255)
                                                       tttGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTctgct                                                                                            <  1:4237915/150‑1 (MQ=255)
                                                       tttGCAGTCGATTAACCAGCTCCTGAATATTGCTGGTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTctgct                                                                                            <  2:4772521/150‑1 (MQ=255)
                                                        ttGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTctgctg                                                                                           <  2:4857395/150‑1 (MQ=255)
                                                        ttGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTctgctg                                                                                           <  2:2068609/150‑1 (MQ=255)
                                                        ttGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTctgctg                                                                                           <  2:122613/150‑1 (MQ=255)
                                                           cAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTgcgc                                                                                        >  2:4515738/1‑150 (MQ=255)
                                                            aGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGcc                                                                                       >  2:2119367/1‑150 (MQ=255)
                                                             gTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCa                                                                                      <  2:2596309/150‑1 (MQ=255)
                                                                gATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGAt                                                                                   >  1:3710296/1‑150 (MQ=255)
                                                                  ttAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTg                                                                                 <  2:3611268/150‑1 (MQ=255)
                                                                   tAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTgc                                                                                <  2:4284300/150‑1 (MQ=255)
                                                                    aaCCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTgcg                                                                               >  1:4026456/1‑150 (MQ=255)
                                                                    aaCCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTgcg                                                                               >  2:2984071/1‑150 (MQ=255)
                                                                    aaCCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTgcg                                                                               >  2:2195695/1‑150 (MQ=255)
                                                                     aCCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTgcgc                                                                              >  2:903993/1‑150 (MQ=255)
                                                                        aGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGAAtt                                                                           <  1:1125623/150‑1 (MQ=255)
                                                                          cTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTaa                                                                         >  1:1012685/1‑150 (MQ=255)
                                                                               gAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGCCATATGCTGCGCCCGATTGCGCCTTAAAGTAt                                                                    >  1:1600086/1‑150 (MQ=255)
                                                                                aaTATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATc                                                                   >  1:3638689/1‑150 (MQ=255)
                                                                                 atatTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCg                                                                  >  1:3650365/1‑150 (MQ=255)
                                                                                 atatTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCg                                                                  >  1:2058871/1‑150 (MQ=255)
                                                                                   atTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAAAGTATCGGc                                                                >  2:4080407/1‑150 (MQ=255)
                                                                                    ttGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGcc                                                               >  2:4352426/1‑150 (MQ=255)
                                                                                    ttGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTGCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGcc                                                               >  1:2663655/1‑150 (MQ=255)
                                                                                     tGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCCTTAACGTATCGGCCt                                                              >  2:3959565/1‑150 (MQ=255)
                                                                                     tGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCt                                                              <  2:1351697/150‑1 (MQ=255)
                                                                                     tGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCt                                                              <  1:2315092/150‑1 (MQ=255)
                                                                                        tGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGc                                                           >  1:2808658/1‑150 (MQ=255)
                                                                                        tGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGc                                                           >  1:963691/1‑150 (MQ=255)
                                                                                        tGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCAGCTGCGCCAGATTGCGCATTAACCTATCGGCCAGGc                                                           >  1:1410355/1‑150 (MQ=255)
                                                                                         gTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCt                                                          <  1:2229128/150‑1 (MQ=255)
                                                                                         gTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCt                                                          <  1:4502635/150‑1 (MQ=255)
                                                                                           tCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTGCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGc                                                        >  2:1928959/1‑150 (MQ=255)
                                                                                            ccACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCg                                                       >  1:308037/1‑150 (MQ=255)
                                                                                            ccACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCg                                                       <  1:346543/150‑1 (MQ=255)
                                                                                              aCTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGtt                                                     <  1:4056715/150‑1 (MQ=255)
                                                                                               cTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTATGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTg                                                    >  1:2197938/1‑150 (MQ=255)
                                                                                                   gcACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCg                                                <  1:548469/150‑1 (MQ=255)
                                                                                                   gcACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCg                                                <  2:127439/150‑1 (MQ=255)
                                                                                                    cACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCgg                                               >  1:3190163/1‑150 (MQ=255)
                                                                                                       gcgGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCCGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCt                                            >  1:2706173/1‑150 (MQ=255)
                                                                                                        cgGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTATGCTGCGCCAGATTGCGCATTAACGTATAGGCCTGGCTGCGTTGTCCGGccgg                                           >  1:2044488/1‑148 (MQ=255)
                                                                                                         ggCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGc                                          >  1:1780564/1‑150 (MQ=255)
                                                                                                          gCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGcc                                         <  2:2739406/150‑1 (MQ=255)
                                                                                                           cAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCa                                        <  2:2045689/150‑1 (MQ=255)
                                                                                                               cTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAgaga                                    <  2:350140/150‑1 (MQ=255)
                                                                                                                 aTTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATc                                  <  1:1871465/150‑1 (MQ=255)
                                                                                                                     atatGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGc                              <  2:307520/150‑1 (MQ=255)
                                                                                                                       atGCGCCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGa                            <  2:2412636/150‑1 (MQ=255)
                                                                                                                         gcgcCAGCGCCGCTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCACATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGaaa                          <  1:1091830/150‑1 (MQ=255)
                                                                                                                          cgcCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAg                         <  2:417194/150‑1 (MQ=255)
                                                                                                                           gcCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGt                        <  2:4456385/150‑1 (MQ=255)
                                                                                                                                cgccgcTCTGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGgta                   <  2:70796/150‑1 (MQ=255)
                                                                                                                                      tctGTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAAt             <  2:167017/150‑1 (MQ=255)
                                                                                                                                         gTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCa          <  2:740564/150‑1 (MQ=255)
                                                                                                                                         gTCCTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCa          <  2:15018/150‑1 (MQ=255)
                                                                                                                                         gTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCa          <  1:382224/150‑1 (MQ=255)
                                                                                                                                          tCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCAt         <  2:4418962/150‑1 (MQ=255)
                                                                                                                                           ccTGGTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATa        <  1:2429826/150‑1 (MQ=255)
                                                                                                                                              ggTCATTACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCCTACGc     >  1:715696/1‑150 (MQ=255)
                                                                                                                                                 cATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGc  >  2:1079911/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCTGCCTCTTTGCCGCTTTCTCGCAGGCGGTTAGCGGTTTCCAGCACCTGATCGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGC  >  CP009273/3680762‑3681056

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