Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,245,537 G→T 91.4% N286K (AAC→AAA plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,245,5370GT91.4% 526.3 / 43.3 199N286K (AAC→AAAplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base T (67/115);  total (76/123)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 203 total aligned reads displayed.

GGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245392‑4245681
                                                                                                                                                 |                                                                                                                                                
ggTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTTATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                                <  2:810620/148‑1 (MQ=255)
ggTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                                >  2:2887703/1‑148 (MQ=255)
  ttGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTAt                                                                                                                                              <  2:2433909/148‑1 (MQ=255)
   tGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATg                                                                                                                                             >  1:386502/1‑148 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  2:2345945/148‑1 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  1:84356/148‑1 (MQ=255)
     aTCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGa                                                                                                                                           <  1:3427421/148‑1 (MQ=255)
     aTCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGa                                                                                                                                           <  2:3414938/148‑1 (MQ=255)
          ggcggcGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTg                                                                                                                                      <  1:3324191/148‑1 (MQ=255)
              gcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGc                                                                                                                                  >  1:391646/1‑148 (MQ=255)
               cgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGcg                                                                                                                                 <  1:2841373/148‑2 (MQ=255)
                 atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               <  2:2614201/148‑1 (MQ=255)
                 atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               <  2:2812572/148‑1 (MQ=255)
                 atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               >  1:2548359/1‑148 (MQ=255)
                  tatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTg                                                                                                                              <  2:1288123/148‑1 (MQ=255)
                      cggcggCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCGAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTaa                                                                                                                          <  2:2255764/148‑1 (MQ=255)
                       ggcggcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTaaa                                                                                                                         <  2:1908844/148‑1 (MQ=255)
                           gcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCg                                                                                                                     <  1:3089811/148‑1 (MQ=255)
                            cACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCgg                                                                                                                    >  2:1621402/1‑148 (MQ=255)
                            cACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCAAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCgg                                                                                                                    <  2:2619646/148‑1 (MQ=255)
                                aGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTcc                                                                                                                <  1:1918855/148‑1 (MQ=255)
                                 gCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCa                                                                                                               <  2:1391383/148‑1 (MQ=255)
                                 gCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCa                                                                                                               <  2:1091622/148‑1 (MQ=255)
                                  ccccTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAg                                                                                                              >  1:844427/1‑148 (MQ=255)
                                     cTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTg                                                                                                           <  1:3151878/148‑1 (MQ=255)
                                      tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCatgctc                                                                                                          >  1:1405953/1‑143 (MQ=255)
                                      tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGa                                                                                                          <  2:482829/148‑1 (MQ=255)
                                        gTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAg                                                                                                        >  1:1251317/1‑148 (MQ=255)
                                          gATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGcc                                                                                                      >  1:1690977/1‑148 (MQ=255)
                                          gATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGcc                                                                                                      >  1:2440289/1‑148 (MQ=255)
                                           aTACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGccc                                                                                                     <  1:3139749/148‑1 (MQ=255)
                                             aCAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAg                                                                                                   >  1:1146494/1‑148 (MQ=255)
                                              cAGCACGTAAGAAAGCAGCAGGTAGTCCTTGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGa                                                                                                  <  2:408489/148‑1 (MQ=255)
                                               aGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGaa                                                                                                 <  1:2906917/148‑1 (MQ=255)
                                                 cACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATa                                                                                               <  1:3480144/148‑1 (MQ=255)
                                                  aCGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATAc                                                                                              <  2:1865999/148‑1 (MQ=255)
                                                      aaGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtc                                                                                          <  1:1911216/148‑1 (MQ=255)
                                                        gAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcag                                                                                        >  1:157941/1‑148 (MQ=255)
                                                             cagcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAgg                                                                                   <  2:2556383/148‑1 (MQ=255)
                                                                  ggTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAat                                                                              <  2:1545129/148‑1 (MQ=255)
                                                                    tAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatca                                                                            <  2:1636588/148‑1 (MQ=255)
                                                                    tAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatca                                                                            <  1:3455937/148‑1 (MQ=255)
                                                                    tAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatca                                                                            <  1:2216409/148‑1 (MQ=255)
                                                                     aGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatcat                                                                           <  1:682203/148‑1 (MQ=255)
                                                                       tCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtct                                                                         >  1:361054/1‑148 (MQ=255)
                                                                          aTGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                                      >  2:3503901/1‑148 (MQ=255)
                                                                           tgtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTa                                                                     >  2:3297615/1‑148 (MQ=255)
                                                                               aCTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAagag                                                                 <  1:268217/148‑1 (MQ=255)
                                                                                cTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAgaga                                                                >  2:2999515/1‑148 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               <  1:2755382/148‑1 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               <  1:2755606/148‑1 (MQ=255)
                                                                                   cGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAt                                                             <  2:232193/148‑1 (MQ=255)
                                                                                   cGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAt                                                             <  1:2314756/148‑1 (MQ=255)
                                                                                    ggTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAtt                                                            <  2:14221/148‑1 (MQ=255)
                                                                                      tGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcg                                                          <  1:2439474/148‑1 (MQ=255)
                                                                                           aGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCa                                                     <  2:2291888/148‑1 (MQ=255)
                                                                                           aGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCa                                                     <  2:2273364/148‑1 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   >  1:1470253/1‑148 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   <  2:2979796/148‑1 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   >  1:87106/1‑148 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   >  2:1386623/1‑148 (MQ=255)
                                                                                                caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtctc                                                <  1:2660528/148‑1 (MQ=255)
                                                                                                caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtctc                                                >  1:2174976/1‑148 (MQ=255)
                                                                                                  tataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCtt                                              >  2:2896613/1‑148 (MQ=255)
                                                                                                   ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                             <  2:1323597/148‑1 (MQ=255)
                                                                                                   ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                             >  1:2656652/1‑148 (MQ=255)
                                                                                                       cAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                                         <  1:3244571/148‑1 (MQ=255)
                                                                                                        aGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                                        <  2:1212891/148‑1 (MQ=255)
                                                                                                         gCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAg                                       <  2:386502/148‑1 (MQ=255)
                                                                                                         gCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAg                                       <  2:380895/148‑1 (MQ=255)
                                                                                                          cTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGt                                      >  1:2625308/1‑148 (MQ=255)
                                                                                                           tCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAGGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTg                                     >  2:3060780/1‑148 (MQ=255)
                                                                                                           tCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTg                                     <  1:722902/148‑1 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                  <  2:2382929/148‑1 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                  <  2:2832776/148‑1 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                  <  1:2533927/148‑1 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 >  1:1320157/1‑148 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  1:2225081/148‑1 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  1:2942971/148‑1 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  2:2987825/148‑1 (MQ=255)
                                                                                                                   gtcgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                             <  2:3382091/148‑1 (MQ=255)
                                                                                                                    tcgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGt                            <  1:229028/148‑1 (MQ=255)
                                                                                                                     cgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                           >  1:3001574/1‑148 (MQ=255)
                                                                                                                         ggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                       <  1:3073506/148‑1 (MQ=255)
                                                                                                                               cTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCt                 >  1:3367048/1‑148 (MQ=255)
                                                                                                                                    gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc            <  2:2332641/148‑1 (MQ=255)
                                                                                                                                     aaCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCa           <  2:2637004/148‑1 (MQ=255)
                                                                                                                                      aCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCGGCGCTTTTTCAt          >  1:2325072/1‑148 (MQ=255)
                                                                                                                                       cgcgCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATg         >  1:1085694/1‑148 (MQ=255)
                                                                                                                                         cgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATggg       >  2:1957381/1‑148 (MQ=255)
                                                                                                                                         cgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATggg       <  1:3422642/148‑1 (MQ=255)
                                                                                                                                           cTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgag     >  2:94562/1‑148 (MQ=255)
                                                                                                                                            tCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgaga    >  1:492629/1‑148 (MQ=255)
                                                                                                                                            tCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgaga    >  1:3060941/1‑148 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   <  1:1870234/148‑1 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   >  2:2122962/1‑148 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   >  2:2121634/1‑148 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   <  2:1718010/148‑1 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   <  1:804477/148‑1 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   <  2:943216/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245392‑4245681

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: