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,3450CA96.6% 399.4 / 4.4 131R751L (CGC→CTC) bcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base A (58/68);  total (59/72)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 132 total aligned reads displayed.

GCCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATC  >  CP009273/3686199‑3686485
                                                                                                                                                  |                                                                                                                                            
 ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                            >  1:17437/1‑148 (MQ=255)
 ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                            <  1:2418250/148‑1 (MQ=255)
 ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                            <  1:2901073/148‑1 (MQ=255)
    gagCGCCCATGTATGCGCACGGGCAAACGTACTCTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGAAGTATTCCTGCTGACCGAGTTTa                                                                                                                                         <  1:777176/148‑1 (MQ=255)
    gagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                         <  2:251666/148‑1 (MQ=255)
     agCGCCCATGTATCCGCACGGGCAATCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTaa                                                                                                                                        >  2:825985/1‑148 (MQ=255)
     agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTaa                                                                                                                                        >  2:825987/1‑148 (MQ=255)
       cgcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACGTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGc                                                                                                                                      <  2:2916305/148‑1 (MQ=255)
        gcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCa                                                                                                                                     <  1:1396227/148‑1 (MQ=255)
        gcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCa                                                                                                                                     <  1:2174723/148‑1 (MQ=255)
        gcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCa                                                                                                                                     <  1:2174648/148‑1 (MQ=255)
         cccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGACGACGTATTCCTGCTGACCGAGTTTAAGCaa                                                                                                                                    >  1:467205/1‑148 (MQ=255)
         cccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCaa                                                                                                                                    >  1:2258598/1‑148 (MQ=255)
          ccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAAc                                                                                                                                   <  2:2961558/148‑1 (MQ=255)
          ccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAAc                                                                                                                                   <  2:1565156/148‑1 (MQ=255)
          ccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAAc                                                                                                                                   <  2:1565197/148‑1 (MQ=255)
            aTGTATCCGCACGGGCAAACGTACACTTCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAAATTCATTACCCATCCCGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAg                                                                                                                                 <  2:3203491/148‑1 (MQ=255)
            aTGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAg                                                                                                                                 <  1:2965721/148‑1 (MQ=255)
            aTGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAg                                                                                                                                 <  1:1105951/148‑1 (MQ=255)
            aTGGATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAg                                                                                                                                 <  2:2267988/148‑1 (MQ=255)
              gTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGAt                                                                                                                               <  1:1850030/148‑1 (MQ=255)
               tATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGAtt                                                                                                                              >  2:2469479/1‑148 (MQ=255)
                 tCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCa                                                                                                                            <  2:718269/148‑1 (MQ=255)
                   cGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACt                                                                                                                          <  1:1164802/148‑1 (MQ=255)
                   cGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACt                                                                                                                          <  1:2489739/148‑1 (MQ=255)
                     cACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                        <  2:2872570/148‑1 (MQ=255)
                     cACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                        <  1:129628/148‑1 (MQ=255)
                     cACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                        <  2:2872453/148‑1 (MQ=255)
                          gCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGcc                                                                                                                   <  1:2123680/148‑1 (MQ=255)
                           cAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGccc                                                                                                                  >  2:214108/1‑148 (MQ=255)
                              aCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTc                                                                                                               <  2:1819127/148‑1 (MQ=255)
                                  acacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGa                                                                                                           <  2:107139/148‑1 (MQ=255)
                                  acacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGa                                                                                                           <  1:2283250/148‑1 (MQ=255)
                                  acacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGa                                                                                                           <  1:343269/148‑1 (MQ=255)
                                   cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTc                                                                                                               <  1:1184282/143‑1 (MQ=255)
                                   cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTc                                                                                                               >  2:1184282/1‑143 (MQ=255)
                                    acTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAAt                                                                                                         <  1:1385388/148‑1 (MQ=255)
                                    acTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAAt                                                                                                         <  2:2282683/148‑1 (MQ=255)
                                    acTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAAt                                                                                                         <  2:651183/148‑1 (MQ=255)
                                     cTGCACAAAATCGATATGTTGCTGGGTTGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATc                                                                                                        <  1:1113643/148‑1 (MQ=255)
                                     cTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATc                                                                                                        >  2:585130/1‑148 (MQ=255)
                                      tGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCt                                                                                                       >  1:619546/1‑148 (MQ=255)
                                        cacaAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgc                                                                                                     >  2:3425592/1‑148 (MQ=255)
                                          caAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcgc                                                                                                   >  1:1297218/1‑148 (MQ=255)
                                          caAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcgc                                                                                                   <  1:1422169/148‑1 (MQ=255)
                                          caAAATCGAGATGGTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTGTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcgc                                                                                                   <  2:2403164/148‑1 (MQ=255)
                                             aaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTg                                                                                                >  1:197093/1‑148 (MQ=255)
                                               tCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGAc                                                                                              >  2:1710269/1‑148 (MQ=255)
                                               tCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCCGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGAc                                                                                              <  2:2992719/148‑1 (MQ=255)
                                                 gATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCg                                                                                            <  2:1615692/148‑1 (MQ=255)
                                                  atatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGt                                                                                           >  1:1530506/1‑148 (MQ=255)
                                                   tatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGtt                                                                                          <  2:468701/148‑1 (MQ=255)
                                                   tatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGtt                                                                                          >  2:2084183/1‑148 (MQ=255)
                                                    atGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                         >  2:2046365/1‑148 (MQ=255)
                                                      gTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgat                                                                                       >  1:3145689/1‑148 (MQ=255)
                                                           tGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttga                                                                                  <  2:321393/148‑1 (MQ=255)
                                                           tGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttga                                                                                  <  1:1176240/148‑1 (MQ=255)
                                                            gggTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgat                                                                                 <  2:377524/148‑1 (MQ=255)
                                                             ggtggtGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgatc                                                                                >  2:1996306/1‑148 (MQ=255)
                                                               tggtgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTACGCAACAGATTCACTTTCTGCCCTTCAAGAATCTGCGCCTGACCGTTGATCTTGATccc                                                                              <  1:1494349/148‑1 (MQ=255)
                                                                ggtgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATcccc                                                                             >  1:3503074/1‑148 (MQ=255)
                                                                     gCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAacc                                                                        <  1:956090/148‑1 (MQ=255)
                                                                      cGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAacca                                                                       >  1:1554517/1‑148 (MQ=255)
                                                                         cATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCg                                                                    >  1:1156104/1‑148 (MQ=255)
                                                                            tAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCg                                                                 >  1:3180425/1‑148 (MQ=255)
                                                                            tAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCg                                                                 >  2:1381771/1‑148 (MQ=255)
                                                                            tAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCg                                                                 >  1:2576836/1‑148 (MQ=255)
                                                                                  cAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGaaa                                                           >  1:3006243/1‑148 (MQ=255)
                                                                                    gCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATc                                                         <  1:769716/148‑1 (MQ=255)
                                                                                     cccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                        >  1:156798/1‑148 (MQ=255)
                                                                                     cccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                        >  2:1103204/1‑148 (MQ=255)
                                                                                     cccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                        <  1:908659/148‑1 (MQ=255)
                                                                                     cccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                        >  2:3515630/1‑148 (MQ=255)
                                                                                      ccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCt                                                       <  1:744715/148‑1 (MQ=255)
                                                                                        aaCTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGa                                                     >  2:1006682/1‑148 (MQ=255)
                                                                                        aaCTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGa                                                     <  2:3462669/148‑1 (MQ=255)
                                                                                           ttCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACg                                                  <  1:2046365/148‑1 (MQ=255)
                                                                                               ttACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTac                                              >  1:342488/1‑148 (MQ=255)
                                                                                                  cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACg                                           >  2:2382009/1‑148 (MQ=255)
                                                                                                  cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACg                                           <  2:123233/148‑1 (MQ=255)
                                                                                                  cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACg                                           <  2:3425574/148‑1 (MQ=255)
                                                                                                     aTCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgaga                                        >  1:1778358/1‑148 (MQ=255)
                                                                                                       cACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaag                                      <  2:17437/148‑1 (MQ=255)
                                                                                                        aCGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaaga                                     <  2:1067027/148‑1 (MQ=255)
                                                                                                         cgcgcgCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAagag                                    >  2:1269771/1‑148 (MQ=255)
                                                                                                                 aCCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCAt                            >  2:654812/1‑148 (MQ=255)
                                                                                                                      ggTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcg                       >  1:2739440/1‑148 (MQ=255)
                                                                                                                        tCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcgcg                     <  2:1363969/148‑1 (MQ=255)
                                                                                                                        tCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcgcg                     <  2:870221/148‑1 (MQ=255)
                                                                                                                        tCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcgcg                     <  2:3101477/148‑1 (MQ=255)
                                                                                                                          gggAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCggg                   <  2:1238191/148‑1 (MQ=255)
                                                                                                                             aaGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCaa                <  2:339653/148‑1 (MQ=255)
                                                                                                                              aGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAt               <  1:3395581/148‑1 (MQ=255)
                                                                                                                                aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCACTTg             >  2:3178491/1‑148 (MQ=255)
                                                                                                                                aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTg             >  1:2377762/1‑148 (MQ=255)
                                                                                                                                aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTg             >  1:192552/1‑148 (MQ=255)
                                                                                                                                    aTTCCTGCTGACCGAGTTTAAGCAACAGAATCACTTTCTGCCCTTCCACAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCAGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAAAGCcgc         >  2:528495/1‑148 (MQ=255)
                                                                                                                                     ttCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgcg        >  1:1078803/1‑148 (MQ=255)
                                                                                                                                       ccTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGACGAGGTGACCATCTTCGCGGGCAATTGCCGCggg      >  1:905569/1‑148 (MQ=255)
                                                                                                                                       ccTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTACAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCggg      >  1:3091982/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GCCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATC  >  CP009273/3686199‑3686485

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: ‑