Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,686,345 C→A 100% R751L (CGC→CTC)  bcsA ← cellulose synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,686,3450CA97.1% 879.3 / 13.6 278R751L (CGC→CTC) bcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base A (137/133);  total (142/136)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 278 total aligned reads displayed.

TGCCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTCA  >  CP009273/3686198‑3686489
                                                                                                                                                   |                                                                                                                                                
  ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGttt                                                                                                                                              <  2:42484/150‑1 (MQ=255)
   cAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                             >  1:4253146/1‑150 (MQ=255)
   cAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                             >  2:5110588/1‑150 (MQ=255)
    agagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTaa                                                                                                                                            <  2:1232928/150‑1 (MQ=255)
      agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGc                                                                                                                                          >  1:3527746/1‑150 (MQ=255)
      agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGc                                                                                                                                          >  1:4703288/1‑150 (MQ=255)
      agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGc                                                                                                                                          <  1:5146026/150‑1 (MQ=255)
         gcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAAc                                                                                                                                       <  1:3842990/150‑1 (MQ=255)
         gcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAAc                                                                                                                                       <  1:2387611/150‑1 (MQ=255)
          cccATGTATCCGCACGGGCAAACTTACACTGCACAAAATCGATTTGTTGCTGGGTGGTGAGCGGCATTAAGTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCATGCTGACCGAGTTTAAGCAACa                                                                                                                                      <  1:4666797/150‑1 (MQ=255)
           ccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGACGACGTATTCCTGCTGACCGAGTTTAAGCAACAg                                                                                                                                     >  2:1513000/1‑150 (MQ=255)
              tGTATCCGCACGGGCAAACGTACACTGCACAAACTCGATATGTTGCTGGGTGGTGAGCGGCATTACTTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGAtt                                                                                                                                  <  2:4281600/150‑1 (MQ=255)
                tATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCa                                                                                                                                >  1:4936232/1‑150 (MQ=255)
                tATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCa                                                                                                                                >  2:2684235/1‑150 (MQ=255)
                  tCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACt                                                                                                                              >  1:3816207/1‑150 (MQ=255)
                    cGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                            <  1:1025860/150‑1 (MQ=255)
                    cGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                            <  2:4054875/150‑1 (MQ=255)
                    cGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                            >  1:1977682/1‑150 (MQ=255)
                     gCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGGGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTc                                                                                                                           >  1:1812892/1‑150 (MQ=255)
                      cACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCt                                                                                                                          >  1:3895148/1‑150 (MQ=255)
                       aCGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTg                                                                                                                         <  2:710116/150‑1 (MQ=255)
                         gggCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGcc                                                                                                                       >  1:3532644/1‑150 (MQ=255)
                             aaaCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTCCCCAGCACGCGCGCCACCTGGGTCGGGACGACGTCTTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCtgc                                                                                                                   >  2:2811714/1‑148 (MQ=255)
                              aaCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTcc                                                                                                                  <  1:3139444/150‑1 (MQ=255)
                                cGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAg                                                                                                                >  2:3019015/1‑150 (MQ=255)
                                cGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAg                                                                                                                >  2:3761688/1‑150 (MQ=255)
                                    cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATc                                                                                                            >  2:588402/1‑150 (MQ=255)
                                    cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATc                                                                                                            >  2:907259/1‑150 (MQ=255)
                                     tcTGCACAAAATCGATATGGTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCt                                                                                                           <  2:3475336/149‑1 (MQ=255)
                                       tGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgc                                                                                                         <  1:3870060/150‑1 (MQ=255)
                                        gCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcg                                                                                                        >  1:4329449/1‑150 (MQ=255)
                                        gCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcg                                                                                                        >  2:1274286/1‑150 (MQ=255)
                                           caAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCt                                                                                                     >  2:4330217/1‑150 (MQ=255)
                                               aTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGAcc                                                                                                 <  2:1025278/150‑1 (MQ=255)
                                                  gATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGtt                                                                                              <  1:1333502/150‑1 (MQ=255)
                                                   atatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                             <  1:4535886/150‑1 (MQ=255)
                                                   atatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                             <  1:61050/150‑1 (MQ=255)
                                                     atGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgat                                                                                           <  2:3354456/150‑1 (MQ=255)
                                                        ttGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatctt                                                                                        >  2:1583038/1‑150 (MQ=255)
                                                              tgtggtgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATccc                                                                                  >  2:189835/2‑150 (MQ=255)
                                                              ggtggtGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATccc                                                                                  <  2:3434202/150‑1 (MQ=255)
                                                              ggtggtGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATccc                                                                                  <  1:3424579/150‑1 (MQ=255)
                                                               gtggtgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATcccc                                                                                 >  1:359999/1‑150 (MQ=255)
                                                                tggtgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCa                                                                                <  1:4087289/150‑1 (MQ=255)
                                                                   tgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAac                                                                             >  1:4228450/1‑150 (MQ=255)
                                                                    gAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAacc                                                                            <  1:2178327/150‑1 (MQ=255)
                                                                    gAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAacc                                                                            <  1:3798815/150‑1 (MQ=255)
                                                                     aGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAacca                                                                           <  1:2930598/150‑1 (MQ=255)
                                                                       cGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAaccacc                                                                         <  1:680053/150‑1 (MQ=255)
                                                                        ggCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCg                                                                        <  1:3171098/150‑1 (MQ=255)
                                                                           aTTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCg                                                                     >  2:783206/1‑150 (MQ=255)
                                                                            ttAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCgg                                                                    >  2:919032/1‑150 (MQ=255)
                                                                             tAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGga                                                                   >  1:1191369/1‑150 (MQ=255)
                                                                                  tCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAAt                                                              >  2:437595/1‑150 (MQ=255)
                                                                                  tCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAAt                                                              >  2:4741879/1‑150 (MQ=255)
                                                                                   cAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATc                                                             >  1:2265324/1‑150 (MQ=255)
                                                                                    aGCCCAACTTCATGACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                            <  2:4666578/150‑1 (MQ=255)
                                                                                       ccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGa                                                         <  1:2051294/150‑1 (MQ=255)
                                                                                        cAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGaa                                                        <  2:3963945/150‑1 (MQ=255)
                                                                                         aaCTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAAc                                                       <  1:5068338/150‑1 (MQ=255)
                                                                                         aaCTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAAc                                                       >  1:4296327/1‑150 (MQ=255)
                                                                                               aTTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTaca                                                 >  2:4173161/1‑150 (MQ=255)
                                                                                                ttACCCCTCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTacac                                                <  2:4892453/150‑1 (MQ=255)
                                                                                                ttACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTacac                                                >  1:3270608/1‑150 (MQ=255)
                                                                                                ttACACATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAATCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCACACCACCGTCGGAGAAATCCTGAACGGTacac                                                <  2:3545289/150‑1 (MQ=255)
                                                                                                   cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgag                                             >  1:858008/1‑150 (MQ=255)
                                                                                                   cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgag                                             <  1:2131804/150‑1 (MQ=255)
                                                                                                    ccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgaga                                            >  2:3656143/1‑150 (MQ=255)
                                                                                                     cATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaa                                           >  1:4385063/1‑150 (MQ=255)
                                                                                                       tCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaaga                                         <  1:3946270/150‑1 (MQ=255)
                                                                                                       tCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaaga                                         <  2:2065360/150‑1 (MQ=255)
                                                                                                            cgcgcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGa                                    <  1:927332/150‑1 (MQ=255)
                                                                                                              cgcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGAcc                                  >  1:2414956/1‑150 (MQ=255)
                                                                                                               gcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCa                                 >  2:3293579/1‑150 (MQ=255)
                                                                                                                 cACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATc                               >  2:2247015/1‑150 (MQ=255)
                                                                                                                 cACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATc                               <  1:4174807/150‑1 (MQ=255)
                                                                                                                    cTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTc                            >  1:1001928/1‑150 (MQ=255)
                                                                                                                    cTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTc                            >  2:1625908/1‑150 (MQ=255)
                                                                                                                    cTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTc                            <  1:3903902/150‑1 (MQ=255)
                                                                                                                    cTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTc                            <  2:2098141/150‑1 (MQ=255)
                                                                                                                     tGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcg                           >  2:1399635/1‑150 (MQ=255)
                                                                                                                       ggTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTCCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcgcg                         <  1:2862542/150‑1 (MQ=255)
                                                                                                                         tCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCggg                       <  2:3078444/150‑1 (MQ=255)
                                                                                                                             gAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAt                   >  2:1088267/1‑150 (MQ=255)
                                                                                                                              aaGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAtt                  <  1:86312/150‑1 (MQ=255)
                                                                                                                               atacttATTCCTGCTGACCTAGTTTAAGCAAATTTTTCACTTTCTGCCCTTCCAGAATATGCGCCTTACGTTTTTTATTTATCCCAAAACCACCGTCGTATAAATCGTGAACGTTACACTATAATTGTTGACCATCTTCGGGGGAAATTg                 <  2:1481868/145‑1 (MQ=255)
                                                                                                                                gACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGc                <  2:2886509/150‑1 (MQ=255)
                                                                                                                                 aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGcc               <  2:663568/150‑1 (MQ=255)
                                                                                                                                  cGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcg              <  1:4433127/150‑1 (MQ=255)
                                                                                                                                   gTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTCTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTTATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTTACCATCTTCGCGGGCAATTGCcgc             >  1:3721528/1‑150 (MQ=255)
                                                                                                                                    tATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgcg            <  2:2394023/150‑1 (MQ=255)
                                                                                                                                    tATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgcg            <  1:1486043/150‑1 (MQ=255)
                                                                                                                                     aTTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCgg           >  2:2542328/1‑150 (MQ=255)
                                                                                                                                          tgctgACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATc      <  1:1619541/150‑1 (MQ=255)
                                                                                                                                          tgctgACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATc      <  1:155629/150‑1 (MQ=255)
                                                                                                                                           gctgACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCg     <  2:3325694/150‑1 (MQ=255)
                                                                                                                                           gctgACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCg     >  2:3115128/1‑150 (MQ=255)
                                                                                                                                             tgACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTc   >  2:1214402/1‑150 (MQ=255)
                                                                                                                                              gACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTNCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTCa  >  1:2682366/1‑150 (MQ=255)
                                                                                                                                              gACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTCa  >  2:2615116/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGCCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTCA  >  CP009273/3686198‑3686489

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