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,5370GT98.4% 766.6 / 3.0 243N286K (AAC→AAAplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (92/147);  total (95/148)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01
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 249 total aligned reads displayed.

GGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245392‑4245681
                                                                                                                                                 |                                                                                                                                                
ggTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                                <  1:979526/148‑1 (MQ=255)
ggTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                                >  1:599328/1‑148 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  1:3320087/148‑1 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  2:2802231/148‑1 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  2:2803073/148‑1 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  1:3219685/148‑1 (MQ=255)
    gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATgg                                                                                                                                            <  2:3273305/148‑1 (MQ=255)
     aTCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGa                                                                                                                                           <  2:1132969/148‑1 (MQ=255)
     aTCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGa                                                                                                                                           <  2:3186712/148‑1 (MQ=255)
     aTCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGa                                                                                                                                           <  1:1364027/148‑1 (MQ=255)
      tCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGAc                                                                                                                                          >  1:87672/1‑148 (MQ=255)
           gcggcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGa                                                                                                                                     <  1:814152/148‑1 (MQ=255)
              gcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGc                                                                                                                                  <  1:1218847/148‑1 (MQ=255)
                 atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               >  1:2771968/1‑148 (MQ=255)
                 atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               >  1:2964424/1‑148 (MQ=255)
                  tatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGTAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTg                                                                                                                              <  1:2678591/148‑1 (MQ=255)
                  tatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTg                                                                                                                              >  1:712014/1‑148 (MQ=255)
                  tatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTg                                                                                                                              >  1:712009/1‑148 (MQ=255)
                  tatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTg                                                                                                                              <  2:2491123/148‑1 (MQ=255)
                       ggcggcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTaaa                                                                                                                         <  1:2386861/148‑1 (MQ=255)
                        gcggcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAg                                                                                                                        <  1:2513333/148‑1 (MQ=255)
                         cggcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGt                                                                                                                       <  2:2216942/148‑1 (MQ=255)
                            cACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCgg                                                                                                                    <  2:2330072/148‑1 (MQ=255)
                             aCCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGt                                                                                                                   >  2:1727014/1‑148 (MQ=255)
                              ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                  >  2:1785154/1‑148 (MQ=255)
                                aGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTcc                                                                                                                <  2:1592286/148‑1 (MQ=255)
                                aGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTcc                                                                                                                <  2:1592591/148‑1 (MQ=255)
                                 gCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCa                                                                                                               <  2:2545088/148‑1 (MQ=255)
                                           aTACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGccc                                                                                                     >  2:701525/1‑148 (MQ=255)
                                               aGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGaa                                                                                                 <  1:2599806/148‑1 (MQ=255)
                                                  aCGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATAc                                                                                              <  2:194213/148‑1 (MQ=255)
                                                   cGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACg                                                                                             >  2:1483148/1‑148 (MQ=255)
                                                    gTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACgg                                                                                            <  2:2450049/148‑1 (MQ=255)
                                                     tAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgt                                                                                           >  2:2811679/1‑148 (MQ=255)
                                                         aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagt                                                                                       >  1:2884261/1‑148 (MQ=255)
                                                          aaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtc                                                                                      <  1:1164332/148‑1 (MQ=255)
                                                          aaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtc                                                                                      >  2:1675753/1‑148 (MQ=255)
                                                          aaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACCgtcagtc                                                                                      <  1:109561/148‑1 (MQ=255)
                                                           agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     >  1:100047/1‑148 (MQ=255)
                                                           agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     >  1:3132677/1‑148 (MQ=255)
                                                            gcagcaGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGtcagtcag                                                                                    <  1:3499027/148‑1 (MQ=255)
                                                             cagcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAgg                                                                                   <  2:3378599/148‑1 (MQ=255)
                                                             cagcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAgg                                                                                   <  2:1640795/148‑1 (MQ=255)
                                                                 agGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGaa                                                                               <  2:2627889/148‑1 (MQ=255)
                                                                 agGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGaa                                                                               <  2:539537/148‑1 (MQ=255)
                                                                  ggTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAat                                                                              >  2:3130742/1‑148 (MQ=255)
                                                                    tAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatca                                                                            >  1:1827790/1‑148 (MQ=255)
                                                                    tAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatca                                                                            >  1:160414/1‑148 (MQ=255)
                                                                    tAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatca                                                                            >  1:3414228/1‑148 (MQ=255)
                                                                           tgtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGGTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTa                                                                     >  2:1032078/1‑148 (MQ=255)
                                                                                cTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAgaga                                                                >  1:1025606/1‑148 (MQ=255)
                                                                                cTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGAAGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAgaga                                                                <  1:2929944/148‑1 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               >  1:2409925/1‑148 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               >  1:2409751/1‑148 (MQ=255)
                                                                                 tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               <  1:1980813/148‑1 (MQ=255)
                                                                                  gCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaaa                                                              <  1:3043037/148‑1 (MQ=255)
                                                                                    ggTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAtt                                                            <  1:3121550/148‑1 (MQ=255)
                                                                                       ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                                         <  2:1517354/148‑1 (MQ=255)
                                                                                        gCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                                        >  2:79294/1‑148 (MQ=255)
                                                                                          aaGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgccgc                                                      >  1:2273514/1‑148 (MQ=255)
                                                                                            ggCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCaa                                                    <  2:3127303/148‑1 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   >  2:2836977/1‑148 (MQ=255)
                                                                                             gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   <  1:1861749/148‑1 (MQ=255)
                                                                                               acaTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtct                                                 <  2:1130257/148‑1 (MQ=255)
                                                                                               acaTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtct                                                 <  1:2461962/148‑1 (MQ=255)
                                                                                                caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtctc                                                <  1:248552/148‑1 (MQ=255)
                                                                                                caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtctc                                                >  2:380350/1‑148 (MQ=255)
                                                                                                 atatACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATctct                                               <  2:2752718/148‑1 (MQ=255)
                                                                                                  tataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCtt                                              <  2:812249/148‑1 (MQ=255)
                                                                                                   ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                             >  1:634532/1‑148 (MQ=255)
                                                                                                   ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                             >  1:2437259/1‑148 (MQ=255)
                                                                                                   ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTAGATGCCATGGTAAAGTCGGTTCCAGGTCAAGCCCAGAATACCGTCAGTCATGCGAATCATCTCGTAAGAGAAATTCGCCGCAATATATTc                                             >  1:2214920/1‑148 (MQ=255)
                                                                                                            cATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGc                                    >  2:1436758/1‑148 (MQ=255)
                                                                                                              tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                  >  2:3338881/1‑148 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  2:1182787/148‑1 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  1:969678/148‑1 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  1:2623181/148‑1 (MQ=255)
                                                                                                               gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAAt                                 <  2:2933234/148‑1 (MQ=255)
                                                                                                                 ccGTCGTGGGCCATCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcg                               <  1:535935/148‑1 (MQ=255)
                                                                                                                 ccGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcg                               <  2:2436752/148‑1 (MQ=255)
                                                                                                                     cgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                           >  1:2425298/1‑148 (MQ=255)
                                                                                                                        gggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCTCGTTctg                        <  2:883156/148‑1 (MQ=255)
                                                                                                                        gggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctg                        <  2:2964424/148‑1 (MQ=255)
                                                                                                                         ggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                       >  2:2687643/1‑148 (MQ=255)
                                                                                                                         ggCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                       >  2:3256969/1‑148 (MQ=255)
                                                                                                                             aGCTTGAGAACGCGCTCAGCTTTATGGACTTTTATGCCCTGGTAAATTCGGTTCCAGGTGTTGACCAGAATACGGTCAGTCAGGCTATTCTTCTCGTATGAGATTTTCGCCGCTATCTCTTACATCAGTGCATTCGCGTTCTGCTtcg                   <  2:1298014/148‑1 (MQ=255)
                                                                                                                              gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  >  2:3171250/1‑148 (MQ=255)
                                                                                                                              gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  <  2:499384/148‑1 (MQ=255)
                                                                                                                              gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  >  1:2147517/1‑148 (MQ=255)
                                                                                                                               cTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCt                 <  1:1507510/148‑1 (MQ=255)
                                                                                                                                tGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  2:444295/148‑1 (MQ=255)
                                                                                                                                 ggCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttt               <  2:1841591/148‑1 (MQ=255)
                                                                                                                                 ggCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttt               <  2:347872/148‑1 (MQ=255)
                                                                                                                                  gCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtttt              >  1:2467513/1‑148 (MQ=255)
                                                                                                                                    gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc            <  2:1540741/148‑1 (MQ=255)
                                                                                                                                          gcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGga      <  2:2634482/148‑1 (MQ=255)
                                                                                                                                          gcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGCCAGGCGANTCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGga      <  1:1785154/148‑1 (MQ=255)
                                                                                                                                            tCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgaga    >  2:422219/1‑148 (MQ=255)
                                                                                                                                             cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   <  1:3379116/148‑1 (MQ=255)
                                                                                                                                              aGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtt  <  1:1652400/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: