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,3450CA99.4% 551.1 / ‑4.2 168R751L (CGC→CTC) bcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (90/77);  total (90/78)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 168 total aligned reads displayed.

CCAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTC  >  CP009273/3686200‑3686488
                                                                                                                                                 |                                                                                                                                                
ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                                <  1:1079505/148‑1 (MQ=255)
ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                                <  1:3366394/148‑1 (MQ=255)
ccAGAGCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGGTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGt                                                                                                                                                <  2:1046166/148‑1 (MQ=255)
   gagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                             <  1:2916319/148‑1 (MQ=255)
   gagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                             <  1:247162/148‑1 (MQ=255)
   gagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                             <  1:3499648/148‑1 (MQ=255)
   gagCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTa                                                                                                                                             <  2:3088969/148‑1 (MQ=255)
    agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTaa                                                                                                                                            >  1:2838676/1‑148 (MQ=255)
    agCGCCCATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTaa                                                                                                                                            >  2:260223/1‑148 (MQ=255)
        cccAGGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATGACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCaa                                                                                                                                        <  2:1939650/148‑1 (MQ=255)
         ccATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAAc                                                                                                                                       <  2:1888327/148‑1 (MQ=255)
          cATGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACa                                                                                                                                      >  2:3424100/1‑148 (MQ=255)
            tGTATCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGa                                                                                                                                    <  2:2107161/148‑1 (MQ=255)
               aTCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTc                                                                                                                                 >  1:1776316/1‑148 (MQ=255)
               aTCCGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTc                                                                                                                                 >  2:2060075/1‑148 (MQ=255)
                 ccGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCACCAGATTCAc                                                                                                                               >  2:2603125/1‑148 (MQ=255)
                 ccGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCAc                                                                                                                               >  1:313062/1‑148 (MQ=255)
                 ccGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCAc                                                                                                                               >  2:638715/1‑148 (MQ=255)
                 ccGCACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCAc                                                                                                                               <  2:1434065/148‑1 (MQ=255)
                    cACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                            >  2:1706102/1‑148 (MQ=255)
                    cACGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACttt                                                                                                                            <  1:2759635/148‑1 (MQ=255)
                     aCGGGCAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTc                                                                                                                           <  1:1050044/148‑1 (MQ=255)
                         gCAAACGTACACTGCACAAAATCGATCTGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGcc                                                                                                                       <  1:3377145/148‑1 (MQ=255)
                          cAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGccc                                                                                                                      <  2:2265459/148‑1 (MQ=255)
                          cAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGccc                                                                                                                      <  2:1392066/148‑1 (MQ=255)
                          cAAACGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGccc                                                                                                                      <  1:2679701/148‑1 (MQ=255)
                             aCGTACACTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTc                                                                                                                   <  1:1513875/148‑1 (MQ=255)
                                 acacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGa                                                                                                               >  1:279550/1‑148 (MQ=255)
                                  cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGaa                                                                                                              <  2:1252808/148‑1 (MQ=255)
                                  cacTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGaa                                                                                                              >  1:1093315/1‑148 (MQ=255)
                                   acTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAAt                                                                                                             <  1:877812/148‑1 (MQ=255)
                                   acTGCACAAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAAt                                                                                                             <  1:1553886/148‑1 (MQ=255)
                                         caAAATCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGCGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTgcgc                                                                                                       <  1:428874/148‑1 (MQ=255)
                                            aaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTg                                                                                                    <  1:260264/148‑1 (MQ=255)
                                            aaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTg                                                                                                    <  1:398980/148‑1 (MQ=255)
                                            aaTCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTg                                                                                                    <  1:1619393/148‑1 (MQ=255)
                                              tCGATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCACCTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGAc                                                                                                  <  2:1558601/148‑1 (MQ=255)
                                                gATATGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCg                                                                                                <  1:3434557/148‑1 (MQ=255)
                                                 atatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGt                                                                                               >  2:812401/1‑148 (MQ=255)
                                                 atatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGt                                                                                               >  2:3137063/1‑148 (MQ=255)
                                                  tatGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACAGtt                                                                                              >  2:1799212/1‑148 (MQ=255)
                                                   atGTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttg                                                                                             <  1:672769/148‑1 (MQ=255)
                                                     gTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgat                                                                                           >  2:3080452/1‑148 (MQ=255)
                                                     gTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgat                                                                                           >  2:2918285/1‑148 (MQ=255)
                                                     gTTGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgat                                                                                           >  2:2918141/1‑148 (MQ=255)
                                                      ttGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatc                                                                                          >  2:1699313/1‑148 (MQ=255)
                                                       tGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatct                                                                                         >  1:1264530/1‑148 (MQ=255)
                                                       tGCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatct                                                                                         >  1:1264523/1‑148 (MQ=255)
                                                        gCTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatctt                                                                                        <  1:2683945/148‑1 (MQ=255)
                                                         cTGGGTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttg                                                                                       <  2:1309354/148‑1 (MQ=255)
                                                           gggTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgat                                                                                     >  2:1072567/1‑148 (MQ=255)
                                                           gggTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgat                                                                                     >  1:1454011/1‑148 (MQ=255)
                                                           gggTGGTGAGCGGCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGttgatcttgat                                                                                     >  1:1616060/1‑148 (MQ=255)
                                                                      ggCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAaccac                                                                         <  1:542477/148‑1 (MQ=255)
                                                                      ggCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAaccac                                                                         <  2:1086219/148‑1 (MQ=255)
                                                                       gCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAaccacc                                                                        >  1:934560/1‑148 (MQ=255)
                                                                       gCATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAaccacc                                                                        >  1:2256365/1‑148 (MQ=255)
                                                                        cATTAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCg                                                                       >  1:862466/1‑148 (MQ=255)
                                                                           tAATTTCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGGTGATCTTGATCCCCAAACCACCGTCg                                                                    >  2:470589/1‑148 (MQ=255)
                                                                              tttCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGgag                                                                 >  1:2192488/1‑148 (MQ=255)
                                                                              tttCAGCCCAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGgag                                                                 >  1:2193923/1‑148 (MQ=255)
                                                                                    cccAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATcc                                                           <  2:1269675/148‑1 (MQ=255)
                                                                                      cAACTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTg                                                         >  2:840223/1‑148 (MQ=255)
                                                                                       aaCTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGa                                                        >  2:1798753/1‑148 (MQ=255)
                                                                                         cTTCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAAc                                                      >  1:3404048/1‑148 (MQ=255)
                                                                                          ttCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACg                                                     >  1:2009042/1‑148 (MQ=255)
                                                                                          ttCATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACg                                                     <  2:2838676/148‑1 (MQ=255)
                                                                                            cATTACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGt                                                   <  2:2660332/148‑1 (MQ=255)
                                                                                              ttACCCATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTac                                                 <  2:3213560/148‑1 (MQ=255)
                                                                                                 cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACg                                              >  1:2993400/1‑148 (MQ=255)
                                                                                                 cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACg                                              >  1:2993393/1‑148 (MQ=255)
                                                                                                 cccATCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACg                                              >  1:2121980/1‑148 (MQ=255)
                                                                                                    aTCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAAACCACCGTCGGAGAAATCCTGAACGGTACACgag                                            <  2:52597/148‑1 (MQ=255)
                                                                                                    aTCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgaga                                           <  2:1002333/148‑1 (MQ=255)
                                                                                                    aTCACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACgaga                                           <  1:3266126/148‑1 (MQ=255)
                                                                                                      cACGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaag                                         >  2:2759797/1‑148 (MQ=255)
                                                                                                       aCGCGCGCCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGagaaga                                        <  1:3201720/148‑1 (MQ=255)
                                                                                                           gcgcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTg                                    <  1:1804463/148‑1 (MQ=255)
                                                                                                           gcgcCACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTg                                    <  2:2770714/148‑1 (MQ=255)
                                                                                                              ccACCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGAcc                                 >  2:2806831/1‑148 (MQ=255)
                                                                                                                aCCTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCAt                               <  2:149384/148‑1 (MQ=255)
                                                                                                                  cTGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCt                             <  1:2362981/148‑1 (MQ=255)
                                                                                                                   tGGGTCGGTAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCtt                            <  2:3047439/148‑1 (MQ=255)
                                                                                                                   tGGGTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCtt                            <  2:3045573/148‑1 (MQ=255)
                                                                                                                     ggTCGGGAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCttag                          >  2:1609014/1‑146 (MQ=255)
                                                                                                                         gggAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCggg                      <  1:3311663/148‑1 (MQ=255)
                                                                                                                          ggAAGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGc                     >  1:1411884/1‑148 (MQ=255)
                                                                                                                             aGACGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAAt                  >  2:2614580/1‑148 (MQ=255)
                                                                                                                                cGTATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGc               >  1:3332818/1‑148 (MQ=255)
                                                                                                                                  tATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcg             <  2:1307827/148‑1 (MQ=255)
                                                                                                                                  tATTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcg             >  1:1744507/1‑148 (MQ=255)
                                                                                                                                   aTTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgc            >  2:3428243/1‑148 (MQ=255)
                                                                                                                                   aTTCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgc            >  1:972447/1‑148 (MQ=255)
                                                                                                                                    ttCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgcg           <  2:2399563/148‑1 (MQ=255)
                                                                                                                                    ttCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgcg           <  1:1873745/148‑1 (MQ=255)
                                                                                                                                    ttCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCcgcg           <  2:2399645/148‑1 (MQ=255)
                                                                                                                                     tCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCgg          >  1:2718163/1‑148 (MQ=255)
                                                                                                                                     tCCTGCTGACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCgg          >  1:2125696/1‑148 (MQ=255)
                                                                                                                                            gACCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGt   >  2:1432910/1‑148 (MQ=255)
                                                                                                                                             aCCGAGTTTAAGCAACAGATTCACTTTCTGCCCTTCCAGAATCTGCGCCTGACCGTTGATCTTGATCCCCAAACCACCGTCGGAGAAATCCTGAACGGTACACGAGAAGAGGTGACCATCTTCGCGGGCAATTGCCGCGGGCATCGTc  >  1:767328/1‑148 (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: