Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 249,629 G→A 8.3% R145Q (CGA→CAA)  yafP → putative acyl‑CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273249,6290GA 8.3% 694.2 / 38.4 220R145Q (CGA→CAA) yafPputative acyl‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (112/84);  new base A (10/8);  total (122/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 224 total aligned reads displayed.

GTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTAATTCT  >  CP009273/249483‑249777
                                                                                                                                                  |                                                                                                                                                    
 ttAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGata                                                                                                                                                   >  1:2553739/1‑149 (MQ=255)
  tAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGatat                                                                                                                                                  <  1:2165262/149‑1 (MQ=255)
   aaaaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGAtata                                                                                                                                                 >  1:2828281/1‑149 (MQ=255)
   aaaaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGAtata                                                                                                                                                 >  1:202159/1‑149 (MQ=255)
   aaaaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCAAtata                                                                                                                                                 >  2:2086047/1‑149 (MQ=255)
    aaaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATaa                                                                                                                                                >  1:590281/1‑149 (MQ=255)
     aaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATaaa                                                                                                                                               <  1:5111862/149‑1 (MQ=255)
     aaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATaaa                                                                                                                                               <  1:3099541/149‑1 (MQ=255)
     aaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCAATATaaa                                                                                                                                               >  2:4999740/1‑149 (MQ=255)
               tAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACAtt                                                                                                                                     <  2:1608539/149‑1 (MQ=255)
                 aGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaa                                                                                                                                   >  1:183439/1‑149 (MQ=255)
                 aGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaa                                                                                                                                   >  2:4688242/1‑149 (MQ=255)
                 aGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaa                                                                                                                                   >  2:2363917/1‑149 (MQ=255)
                  gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaaa                                                                                                                                  >  1:284201/1‑149 (MQ=255)
                  gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaaa                                                                                                                                  >  1:489693/1‑149 (MQ=255)
                  gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaaa                                                                                                                                  >  2:1845758/1‑149 (MQ=255)
                   tCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAAt                                                                                                                                 <  1:3695809/149‑1 (MQ=255)
                    cTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATc                                                                                                                                <  2:1412585/149‑1 (MQ=255)
                     tGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATcc                                                                                                                               >  1:4552244/1‑149 (MQ=255)
                      gAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCa                                                                                                                              >  1:4871518/1‑149 (MQ=255)
                       aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAg                                                                                                                             >  1:3486037/1‑149 (MQ=255)
                          ccGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCtt                                                                                                                          <  1:844515/149‑1 (MQ=255)
                             aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCa                                                                                                                       >  1:879295/1‑149 (MQ=255)
                              aCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCaa                                                                                                                      >  1:209688/1‑149 (MQ=255)
                               cTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAAt                                                                                                                     >  2:3067275/1‑149 (MQ=255)
                               cTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAAt                                                                                                                     <  1:2648748/149‑1 (MQ=255)
                                  aCGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGaa                                                                                                                  <  1:1737852/149‑1 (MQ=255)
                                  aCGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGaa                                                                                                                  <  1:975081/149‑1 (MQ=255)
                                    ggTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGaaaa                                                                                                                <  1:1106489/149‑1 (MQ=255)
                                      tGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAGTa                                                                                                              >  1:3101270/1‑149 (MQ=255)
                                      tGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATa                                                                                                              >  1:2250326/1‑149 (MQ=255)
                                      tGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATa                                                                                                              >  1:3263511/1‑149 (MQ=255)
                                       ggACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATaa                                                                                                             <  2:2305696/149‑1 (MQ=255)
                                       ggACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATaa                                                                                                             <  2:1386445/149‑1 (MQ=255)
                                        gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAc                                                                                                            >  2:4408758/1‑149 (MQ=255)
                                           gCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCAATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                                                                                                         >  2:3457591/1‑149 (MQ=255)
                                                cATAACCGCAAAACCCTTTTTTGAAGGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCt                                                                                                    <  2:2658248/149‑1 (MQ=255)
                                                cATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCt                                                                                                    <  2:58958/149‑1 (MQ=255)
                                                cATAACCGCAAAACCATTTTTTGAACTTTATGGTTTTCATAAAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCt                                                                                                    <  2:1449799/149‑1 (MQ=255)
                                                 aTAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTg                                                                                                   >  1:4711598/1‑149 (MQ=255)
                                                  tAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTgg                                                                                                  >  1:3516042/1‑149 (MQ=255)
                                                   aaCCGCAAAACCCTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTa                                                                                                >  1:4373891/1‑149 (MQ=255)
                                                         aaaaCCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgc                                                                                           >  1:4121680/1‑149 (MQ=255)
                                                          aaaCCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcg                                                                                          >  1:854317/1‑149 (MQ=255)
                                                          aaaCCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcg                                                                                          >  2:768488/1‑149 (MQ=255)
                                                          aaaCCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcg                                                                                          >  1:4844941/1‑149 (MQ=255)
                                                          aaaCCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcg                                                                                          >  1:1671011/1‑149 (MQ=255)
                                                               cTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGttt                                                                                     <  2:1677483/149‑1 (MQ=255)
                                                               cTTTTTTGAACGTGATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGttt                                                                                     <  2:4294635/149‑1 (MQ=255)
                                                                  ttttGAACGTTATGGTTTTCAGACATTTAAGCAGCAGCGCGTTGAATGCAGGGGAGCGTTGTTTACTAATATTTATATGCGATAAAAACCGCAACATAAAATCCAGCTTGCAAAGAAAATAACTACAGCCAGGAATGTGCGCGTTATCt                                                                                  >  2:643360/1‑149 (MQ=255)
                                                                    ttGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGa                                                                                >  1:3038253/1‑149 (MQ=255)
                                                                         cGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCca                                                                           <  2:3900782/149‑1 (MQ=255)
                                                                           ttATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCcaca                                                                         >  2:4038721/1‑149 (MQ=255)
                                                                            tATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACaa                                                                        >  2:2031660/1‑149 (MQ=255)
                                                                             aTGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACaaa                                                                       >  2:3577063/1‑149 (MQ=255)
                                                                             aTGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACaaa                                                                       >  1:2115873/1‑149 (MQ=255)
                                                                              tGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAAc                                                                      >  1:1426996/1‑149 (MQ=255)
                                                                                gTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTg                                                                    >  1:4020568/1‑149 (MQ=255)
                                                                                gTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTg                                                                    >  1:1317008/1‑149 (MQ=255)
                                                                                gTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGAAATGAAAATAACGCCCGCATGGTATGTGCGCGTTTTCTGATTCCACAAACTg                                                                    >  1:1400424/1‑149 (MQ=255)
                                                                                            ttAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTaaa                                                        >  2:4154467/1‑149 (MQ=255)
                                                                                              aaGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATc                                                      <  2:712023/149‑1 (MQ=255)
                                                                                                 cagcagcGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATg                                                   <  2:2264816/149‑1 (MQ=255)
                                                                                                 cagcagcGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATg                                                   <  2:1282512/149‑1 (MQ=255)
                                                                                                   gcagcGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGAc                                                 >  1:757999/1‑149 (MQ=255)
                                                                                                    cagcGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACa                                                <  2:1096568/149‑1 (MQ=255)
                                                                                                      gcgcGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACaaa                                              >  2:4848236/1‑149 (MQ=255)
                                                                                                      gcgcGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACaaa                                              >  2:2231966/1‑149 (MQ=255)
                                                                                                           ttGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCtt                                         >  1:2036822/1‑149 (MQ=255)
                                                                                                              aaTGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTAt                                      <  1:2952455/149‑1 (MQ=255)
                                                                                                                 gCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCaa                                   <  2:2968997/149‑1 (MQ=255)
                                                                                                                   cGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAAAATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATg                                 <  2:1191960/149‑1 (MQ=255)
                                                                                                                    ggggAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGa                                >  1:2765184/1‑149 (MQ=255)
                                                                                                                    ggggAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGa                                <  2:4708447/149‑1 (MQ=255)
                                                                                                                    ggggAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGa                                >  1:4329281/1‑149 (MQ=255)
                                                                                                                     gggAGCGTGTTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGAc                               <  1:5199243/149‑1 (MQ=255)
                                                                                                                      ggAGCGTGTTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGAct                              <  1:3178507/149‑1 (MQ=255)
                                                                                                                      ggAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGAct                              >  1:1635421/1‑149 (MQ=255)
                                                                                                                        aGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGActct                            <  2:3038253/149‑1 (MQ=255)
                                                                                                                        aGCGTGGTTTAATAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGActct                            <  1:1550547/149‑1 (MQ=255)
                                                                                                                         gCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCtt                           >  1:3318371/1‑149 (MQ=255)
                                                                                                                         gCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCtt                           <  2:2517549/149‑1 (MQ=255)
                                                                                                                             ggTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCa                       <  1:3902072/149‑1 (MQ=255)
                                                                                                                             ggTTTACTAATTTTTATATGCAATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCa                       <  2:1938456/149‑1 (MQ=255)
                                                                                                                                 tACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGAcc                   >  1:4199551/1‑149 (MQ=255)
                                                                                                                                  aCTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCt                  <  2:1699017/149‑1 (MQ=255)
                                                                                                                                  aCTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCt                  <  2:1699237/149‑1 (MQ=255)
                                                                                                                                     aaTTTTTATATGCTATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTc               <  2:3167476/149‑1 (MQ=255)
                                                                                                                                     aaTTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTc               >  1:1744867/1‑149 (MQ=255)
                                                                                                                                       tttttATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCa             >  2:4702279/1‑149 (MQ=255)
                                                                                                                                        ttttATATGCAATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAg            >  1:2341186/1‑149 (MQ=255)
                                                                                                                                          ttATATGAAATATAAACCGCAACATTAAATCCAGCTTGAAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGa          <  1:4999740/149‑1 (MQ=255)
                                                                                                                                           tataTGCTATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCTGACCTTTCCAGGAt         <  2:4488832/149‑1 (MQ=255)
                                                                                                                                           tataTGCGATATAAACCGAAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGAt         <  1:4126234/149‑1 (MQ=255)
                                                                                                                                           tataTGCAATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGAt         <  1:3457591/149‑1 (MQ=255)
                                                                                                                                             tttGCTATATAAACCGCAACATTAAATCCAGCTTGCAATTAAAAAAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTa       <  1:2663038/147‑1 (MQ=255)
                                                                                                                                              atGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTaa      >  1:4956606/1‑149 (MQ=255)
                                                                                                                                                gCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTAAtt    >  2:3510232/1‑149 (MQ=255)
                                                                                                                                                  gATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTAATTCt  >  2:4319964/1‑149 (MQ=255)
                                                                                                                                                  gATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTAATTCt  >  2:1703603/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACAAACTGCAAGGAGGTAAATCATGACAAATCCTTTATCAATGACTCTTTGCAGACCTTTCCAGGATTAATTCT  >  CP009273/249483‑249777

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: