Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,686,345 C→A 87.7% 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,3450CA87.7% 409.5 / 54.6 169R751L (CGC→CTC) bcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base A (70/77);  total (78/91)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

CCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTC  >  CP009273/3686200‑3686488
                                                                                                                                                 |                                                                                                                                               
ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGt                                                                                                                                               <  1:1199001/148‑1 (MQ=255)
ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                               <  1:305914/148‑1 (MQ=255)
ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                               <  1:305911/148‑1 (MQ=255)
   gagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                            <  1:2404805/148‑1 (MQ=255)
    agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTaa                                                                                                                                           >  1:2223763/1‑148 (MQ=255)
     gcgcCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAg                                                                                                                                          >  2:2291036/1‑148 (MQ=255)
          cATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACa                                                                                                                                     <  1:2792927/148‑1 (MQ=255)
           aTGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAg                                                                                                                                    >  2:2796279/1‑148 (MQ=255)
           aTGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAg                                                                                                                                    >  1:727968/1‑148 (MQ=255)
             gTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGAt                                                                                                                                  >  2:830633/1‑148 (MQ=255)
              tATCCTCACGTGAAATCGTAGACTGCACAAAATCGATATGTTGCTTGGTGGTGAGCGTCTTTAATTTCAGCCCAACTTCATTACCCATCACGCTCGCCACCTTGGTCGGGAAGAAGTATTCCTGCTGTCCGTGTTTAAGCAACAGAtt                                                                                                                                 <  2:2900761/148‑1 (MQ=255)
               aTCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTc                                                                                                                                <  2:2690062/148‑1 (MQ=255)
               aTCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTc                                                                                                                                <  2:1034585/148‑1 (MQ=255)
               aTCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTc                                                                                                                                <  1:3309814/148‑1 (MQ=255)
                tCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCa                                                                                                                               <  1:577329/148‑1 (MQ=255)
                 ccGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCAc                                                                                                                              <  1:1389455/148‑1 (MQ=255)
                 ccGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCAc                                                                                                                              <  2:130402/148‑1 (MQ=255)
                  cGCACGGTCAAACTTACACTGCACAAAATCGATATGTTGATTGGTGGTTAGCGGCATTAATTTCAGCCGAACTTCATTAGCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACt                                                                                                                             <  1:390090/148‑1 (MQ=255)
                  cGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACt                                                                                                                             <  1:2758476/148‑1 (MQ=255)
                     aCGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTc                                                                                                                          <  2:2225815/148‑1 (MQ=255)
                     aCGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTc                                                                                                                          <  2:29558/148‑1 (MQ=255)
                     aCGGGCAAACGTACACTGCACAAAATCGATATGGTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATGACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTc                                                                                                                          <  2:1908325/148‑1 (MQ=255)
                         gCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGcc                                                                                                                      <  1:3016865/148‑1 (MQ=255)
                          cAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGccc                                                                                                                     >  1:2832606/1‑148 (MQ=255)
                            aaCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCtt                                                                                                                   >  2:2098360/1‑148 (MQ=255)
                             aCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTc                                                                                                                  <  1:1699251/148‑1 (MQ=255)
                                  cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGaa                                                                                                             <  2:2192481/148‑1 (MQ=255)
                                    cTGCACAAAATCGATATGTTGCTGTGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATc                                                                                                           <  2:369993/148‑1 (MQ=255)
                                    cTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATc                                                                                                           <  1:1217228/148‑1 (MQ=255)
                                     tGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCt                                                                                                          >  1:2534351/1‑148 (MQ=255)
                                       cacaAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgc                                                                                                        >  1:2928266/1‑148 (MQ=255)
                                         caAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcgc                                                                                                      >  1:2170638/1‑148 (MQ=255)
                                          aaaaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGcc                                                                                                     <  1:2231776/148‑1 (MQ=255)
                                          aaaaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGcc                                                                                                     <  1:2231247/148‑1 (MQ=255)
                                           aaaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCt                                                                                                    >  1:1183853/1‑148 (MQ=255)
                                            aaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTg                                                                                                   >  2:155011/1‑148 (MQ=255)
                                            aaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTg                                                                                                   >  1:3275894/1‑148 (MQ=255)
                                              tCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGAc                                                                                                 >  1:2503506/1‑148 (MQ=255)
                                              tCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGAc                                                                                                 >  1:1861514/1‑148 (MQ=255)
                                                gATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCg                                                                                               <  2:1629756/148‑1 (MQ=255)
                                                 atatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAATACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGt                                                                                              <  2:795749/148‑1 (MQ=255)
                                                  tatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGtt                                                                                             >  2:665639/1‑148 (MQ=255)
                                                   atGTTGCTGGGTGGTTAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                            >  1:719489/1‑148 (MQ=255)
                                                   atGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                            <  2:2467760/148‑1 (MQ=255)
                                                   atGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                            >  2:1414303/1‑148 (MQ=255)
                                                    tGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttga                                                                                           >  1:2650580/1‑148 (MQ=255)
                                                      ttGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatc                                                                                         >  1:943352/1‑148 (MQ=255)
                                                           gggTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgat                                                                                    >  1:3259981/1‑148 (MQ=255)
                                                            ggtggtGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgatc                                                                                   <  1:1675473/148‑1 (MQ=255)
                                                             gtggtgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATcc                                                                                  <  2:144581/148‑1 (MQ=255)
                                                                 tgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCaa                                                                              >  2:2538402/1‑148 (MQ=255)
                                                                 tgAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCaa                                                                              >  2:1661338/1‑148 (MQ=255)
                                                                  gAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCaaa                                                                             >  2:2934456/1‑148 (MQ=255)
                                                                   aGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAac                                                                            <  2:605836/148‑1 (MQ=255)
                                                                       gCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAaccacc                                                                        <  1:916832/148‑1 (MQ=255)
                                                                        cATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCg                                                                       >  1:2134743/1‑148 (MQ=255)
                                                                         aTTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGt                                                                      >  1:615433/1‑148 (MQ=255)
                                                                          ttAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTc                                                                     >  2:1829256/1‑148 (MQ=255)
                                                                          ttAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTc                                                                     >  1:243943/1‑148 (MQ=255)
                                                                           tAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCg                                                                    >  2:1696496/1‑148 (MQ=255)
                                                                                   gCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATc                                                            <  1:1246985/148‑1 (MQ=255)
                                                                                   gCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATc                                                            <  2:85399/148‑1 (MQ=255)
                                                                                    cccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                           <  1:912136/148‑1 (MQ=255)
                                                                                     ccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCt                                                          >  1:748586/1‑148 (MQ=255)
                                                                                       aaCTTCATTACCCATCACGCGCGCAACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGa                                                        <  1:3118997/148‑1 (MQ=255)
                                                                                          ttCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACg                                                     >  1:3060398/1‑148 (MQ=255)
                                                                                          ttCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACg                                                     >  2:447449/1‑148 (MQ=255)
                                                                                           tCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACgg                                                    >  1:1546942/1‑148 (MQ=255)
                                                                                           tCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACgg                                                    <  2:2216931/148‑1 (MQ=255)
                                                                                           tCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACgg                                                    <  2:2645084/148‑1 (MQ=255)
                                                                                             aTTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTa                                                  <  1:2994733/148‑1 (MQ=255)
                                                                                               tACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTaca                                                >  1:1712996/1‑148 (MQ=255)
                                                                                                  ccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACga                                             >  1:1066715/1‑148 (MQ=255)
                                                                                                    aTCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgaga                                           >  2:478885/1‑148 (MQ=255)
                                                                                                      cACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaag                                         <  1:51133/148‑1 (MQ=255)
                                                                                                      cACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaag                                         <  1:2934456/148‑1 (MQ=255)
                                                                                                          cgcgcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGt                                     <  2:2266824/148‑1 (MQ=255)
                                                                                                          cgcgcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGt                                     <  2:526075/148‑1 (MQ=255)
                                                                                                               cACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCa                                <  1:1764037/148‑1 (MQ=255)
                                                                                                                aCCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCAt                               <  2:1640958/148‑1 (MQ=255)
                                                                                                                aCCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCAt                               <  2:959838/148‑1 (MQ=255)
                                                                                                                 ccTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATc                              <  1:1698104/148‑1 (MQ=255)
                                                                                                                 ccTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATc                              <  2:503248/148‑1 (MQ=255)
                                                                                                                     ggTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcg                          >  2:667413/1‑148 (MQ=255)
                                                                                                                     ggTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcg                          >  2:974503/1‑148 (MQ=255)
                                                                                                                     ggTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTcg                          >  2:3250744/1‑148 (MQ=255)
                                                                                                                          ggAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGc                     >  1:1116572/1‑148 (MQ=255)
                                                                                                                           gAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCa                    <  2:2910088/148‑1 (MQ=255)
                                                                                                                             aGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAt                  <  2:3058783/148‑1 (MQ=255)
                                                                                                                             aGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAt                  >  2:2930273/1‑148 (MQ=255)
                                                                                                                               aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTg                >  2:1011092/1‑148 (MQ=255)
                                                                                                                               aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTg                <  1:480815/148‑1 (MQ=255)
                                                                                                                               aCGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTg                <  1:2998892/148‑1 (MQ=255)
                                                                                                                                cGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGc               <  2:426053/148‑1 (MQ=255)
                                                                                                                                cGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGc               >  1:1305755/1‑148 (MQ=255)
                                                                                                                                cGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGc               >  1:1304508/1‑148 (MQ=255)
                                                                                                                                  tATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcg             <  2:1512616/148‑1 (MQ=255)
                                                                                                                                  tATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcg             >  1:2164142/1‑148 (MQ=255)
                                                                                                                                       ctgctgACCGAGATTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGc        <  2:3004285/148‑1 (MQ=255)
                                                                                                                                         gctgACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCAt      <  2:2928392/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTC  >  CP009273/3686200‑3686488

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: