Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,245,537 G→T 100% 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,5370GT96.0% 555.9 / 12.3 188N286K (AAC→AAA) plsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base T (80/99);  total (83/105)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 189 total aligned reads displayed.

GGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAT  >  CP009273/4245392‑4245680
                                                                                                                                                 |                                                                                                                                               
ggTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                               <  1:2091381/148‑1 (MQ=255)
 gTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTa                                                                                                                                              <  1:3052703/148‑1 (MQ=255)
  ttGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTAt                                                                                                                                             <  1:616567/148‑1 (MQ=255)
  ttGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTAt                                                                                                                                             <  1:732898/148‑1 (MQ=255)
  ttGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTAt                                                                                                                                             <  1:1012010/148‑1 (MQ=255)
   tGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATg                                                                                                                                            >  2:3007876/1‑148 (MQ=255)
   tGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATg                                                                                                                                            >  2:1107225/1‑148 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                           <  2:1213423/148‑1 (MQ=255)
       cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACg                                                                                                                                        >  2:1458253/1‑148 (MQ=255)
        ccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGATCGCGCTCAGCTTTATGGACGt                                                                                                                                       <  1:2827332/148‑1 (MQ=255)
         cggcggCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGtt                                                                                                                                      <  2:1032695/148‑1 (MQ=255)
          ggcggcGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTg                                                                                                                                     <  1:122095/148‑1 (MQ=255)
           gcggcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGa                                                                                                                                    >  2:1507044/1‑148 (MQ=255)
               cgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGcc                                                                                                                                <  1:601759/148‑1 (MQ=255)
                gATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGccc                                                                                                                               <  1:1127503/148‑1 (MQ=255)
                gATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGccc                                                                                                                               <  1:78402/148‑1 (MQ=255)
                    tGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGt                                                                                                                           <  2:1827642/148‑1 (MQ=255)
                    tGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGt                                                                                                                           <  2:168250/148‑1 (MQ=255)
                              ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                 <  1:2121770/148‑1 (MQ=255)
                              ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                 >  1:3156653/1‑148 (MQ=255)
                              ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                 >  1:2751009/1‑148 (MQ=255)
                              ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                 >  1:1515125/1‑148 (MQ=255)
                               cAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTc                                                                                                                >  1:220370/1‑148 (MQ=255)
                               cAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTc                                                                                                                <  1:2972291/148‑1 (MQ=255)
                               cAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTc                                                                                                                <  1:2887046/148‑1 (MQ=255)
                                   cccTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAgg                                                                                                            <  1:2349216/148‑1 (MQ=255)
                                     cTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTg                                                                                                          >  2:2686187/1‑148 (MQ=255)
                                      tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGa                                                                                                         <  2:172998/148‑1 (MQ=255)
                                      tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGa                                                                                                         >  1:1763679/1‑148 (MQ=255)
                                      tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGa                                                                                                         >  2:762650/1‑148 (MQ=255)
                                       ggTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGaa                                                                                                        >  2:3193179/1‑148 (MQ=255)
                                           aTACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGccc                                                                                                    >  2:21364/1‑148 (MQ=255)
                                             aCAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAg                                                                                                  >  2:503178/1‑148 (MQ=255)
                                                 cACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATa                                                                                              >  1:1627453/1‑148 (MQ=255)
                                                     tAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgt                                                                                          >  2:2916803/1‑148 (MQ=255)
                                                     tAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgt                                                                                          >  2:3159486/1‑148 (MQ=255)
                                                         aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagt                                                                                      >  2:207941/1‑148 (MQ=255)
                                                          aaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtc                                                                                     >  2:2340459/1‑148 (MQ=255)
                                                          aaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtc                                                                                     >  2:2340481/1‑148 (MQ=255)
                                                           agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                    >  1:3058221/1‑148 (MQ=255)
                                                           agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                    >  1:789554/1‑148 (MQ=255)
                                                           agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                    >  2:2913488/1‑148 (MQ=255)
                                                            gcagcaGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGtcagtcag                                                                                   >  1:2596729/1‑148 (MQ=255)
                                                            gcagcaGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGtcagtcag                                                                                   >  2:2505652/1‑148 (MQ=255)
                                                                cagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGa                                                                               >  1:1263180/1‑148 (MQ=255)
                                                                cagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGa                                                                               <  1:2922048/148‑1 (MQ=255)
                                                                   gTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatc                                                                            >  1:1227907/1‑148 (MQ=255)
                                                                     aGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatcat                                                                          <  1:528325/148‑1 (MQ=255)
                                                                     aGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatcat                                                                          <  1:528339/148‑1 (MQ=255)
                                                                      gTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatcatc                                                                         >  1:792462/1‑148 (MQ=255)
                                                                       tCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTCAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtct                                                                        >  2:2021561/1‑148 (MQ=255)
                                                                       tCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtct                                                                        <  2:1786314/148‑1 (MQ=255)
                                                                        ccATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtctc                                                                       <  2:1809972/148‑1 (MQ=255)
                                                                        ccATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtctc                                                                       <  2:3060045/148‑1 (MQ=255)
                                                                            gtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTaa                                                                   >  1:3193196/1‑148 (MQ=255)
                                                                            gtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTaa                                                                   >  1:2883127/1‑148 (MQ=255)
                                                                            gtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTaa                                                                   >  1:1160704/1‑148 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                              <  2:1075148/148‑1 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                              <  1:1063042/148‑1 (MQ=255)
                                                                                  gCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaaa                                                             <  1:1170134/148‑1 (MQ=255)
                                                                                   cGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAt                                                            <  2:1820660/148‑1 (MQ=255)
                                                                                      tGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCTAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATTCGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTct                                                         <  2:616335/148‑2 (MQ=255)
                                                                                      tGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcg                                                         <  2:2278664/148‑1 (MQ=255)
                                                                                       ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                                        >  1:769206/1‑148 (MQ=255)
                                                                                         cAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgccg                                                      >  1:3150202/1‑148 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                  <  2:881486/148‑1 (MQ=255)
                                                                                              cacaTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtc                                                 <  1:2913436/148‑1 (MQ=255)
                                                                                              cacaTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtc                                                 <  1:2913668/148‑1 (MQ=255)
                                                                                               acaTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtct                                                <  1:2835512/148‑1 (MQ=255)
                                                                                                 atatACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATctct                                              >  1:3224062/1‑148 (MQ=255)
                                                                                                 atatACCAGCTCATGACAGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATctct                                              >  1:1739712/1‑148 (MQ=255)
                                                                                                   ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                            <  2:353313/148‑1 (MQ=255)
                                                                                                     aCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCGTCCa                                          <  1:2848794/148‑1 (MQ=255)
                                                                                                      ccAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCAt                                         >  1:2909504/1‑148 (MQ=255)
                                                                                                       cAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                                        >  2:1457689/1‑148 (MQ=255)
                                                                                                        aGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                                       >  1:891330/1‑148 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGTTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                 <  1:1107225/148‑1 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                 <  1:8339/148‑1 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                 <  1:1685549/148‑1 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                 <  2:3242776/148‑1 (MQ=255)
                                                                                                                  cgtcgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgc                             <  2:2564230/148‑1 (MQ=255)
                                                                                                                       tGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTTGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTct                        >  1:2133542/1‑148 (MQ=255)
                                                                                                                       tGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTTGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTct                        >  1:2135401/1‑148 (MQ=255)
                                                                                                                        gggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctg                       <  2:170322/148‑1 (MQ=255)
                                                                                                                         ggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                      <  1:3007876/148‑1 (MQ=255)
                                                                                                                          gCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgct                     >  2:1291434/1‑148 (MQ=255)
                                                                                                                           ccAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgctg                    <  2:2901815/148‑1 (MQ=255)
                                                                                                                           ccAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgctg                    <  2:3088890/148‑1 (MQ=255)
                                                                                                                             aGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcg                  >  2:25976/1‑148 (MQ=255)
                                                                                                                              gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                 <  1:902434/148‑1 (MQ=255)
                                                                                                                                 ggCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttt              <  2:2766471/148‑1 (MQ=255)
                                                                                                                                  gCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtttt             <  1:1458253/148‑1 (MQ=255)
                                                                                                                                    gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc           <  1:2653224/148‑1 (MQ=255)
                                                                                                                                    gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc           <  1:2931060/148‑1 (MQ=255)
                                                                                                                                         cgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATggg      <  1:2517975/148‑1 (MQ=255)
                                                                                                                                         cgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATggg      <  1:1473354/148‑1 (MQ=255)
                                                                                                                                           cTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgag    <  2:1644628/148‑1 (MQ=255)
                                                                                                                                           cTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgag    <  2:1417047/148‑1 (MQ=255)
                                                                                                                                            tCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgaga   <  2:1687924/148‑1 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt  <  1:433050/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAT  >  CP009273/4245392‑4245680

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