Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,224,391 G→T 47.2% V236F (GTT→TTT)  pgi → glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,224,3910GT47.2% 156.6 / 306.2 244V236F (GTT→TTT) pgiglucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (57/72);  new base T (57/58);  total (114/130)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 244 total aligned reads displayed.

TCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCC  >  CP009273/4224251‑4224537
                                                                                                                                            |                                                                                                                                                  
  gcgTAAGTGCTTAAAAAATTAAACCCGGAACCCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAac                                                                                                                                           <  1:732281/148‑1 (MQ=255)
     gAAGTGCTTAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACt                                                                                                                                        <  2:1353874/148‑1 (MQ=255)
      aaGTGCTGAAAAAATTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACtt                                                                                                                                       <  1:1081762/148‑1 (MQ=255)
      aaGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACtt                                                                                                                                       >  1:1846980/1‑148 (MQ=255)
             gAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTgcggcg                                                                                                                                <  2:3168719/148‑1 (MQ=255)
              aaaaaaGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGgcgc                                                                                                                               >  2:666680/1‑148 (MQ=255)
               aaaaaTAAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACAATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCt                                                                                                                              <  1:590727/148‑1 (MQ=255)
               aaaaaGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCt                                                                                                                              <  1:1562832/148‑1 (MQ=255)
                 aaaGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCATAACACTTTGCGGCGCttt                                                                                                                            <  2:103419/148‑1 (MQ=255)
                  aaaaaaaCCCTCAAACCACGTTTTTCTTGGTAGAAAAAAAAACCTACACCACACCGGAAAATATCCCCAACGCCCATAGCGCGCGTGACTGTTACCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTc                                                                                                                           <  1:988843/144‑1 (MQ=255)
                  aaGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTc                                                                                                                           >  2:3072980/1‑148 (MQ=255)
                     tAAACCCTGAATCCACGCTGTTCTTGGTAGAATCTAAAACCTTCACCCCTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAAAACTTTGCAAAACACTTTGCGGCGCTTTccac                                                                                                                        <  2:3508828/148‑1 (MQ=255)
                     tAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTccac                                                                                                                        >  1:2532623/1‑148 (MQ=255)
                     tAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTccac                                                                                                                        >  1:3748876/1‑148 (MQ=255)
                     tAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTGCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTGTGCAAAACACTTTGCGGCGCTTTccac                                                                                                                        <  2:2480672/148‑1 (MQ=255)
                        aCCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCaa                                                                                                                     >  2:1066007/1‑148 (MQ=255)
                        aCCCGCAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCaa                                                                                                                     <  1:81417/148‑1 (MQ=255)
                         cccGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAAt                                                                                                                    <  2:3328363/148‑1 (MQ=255)
                           cGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGc                                                                                                                  >  2:1253875/1‑148 (MQ=255)
                               aaCCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCaaa                                                                                                              >  1:665164/1‑148 (MQ=255)
                               aaCCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCaaa                                                                                                              >  1:2972302/1‑148 (MQ=255)
                                   aCGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCg                                                                                                          >  1:42294/1‑148 (MQ=255)
                                   aCGCTGTTCTTGGTAGCATCAAAAACCTTCCCCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCg                                                                                                          <  1:1643376/148‑1 (MQ=255)
                                       tGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTgg                                                                                                      >  1:942322/1‑148 (MQ=255)
                                        gTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGc                                                                                                     >  2:475639/1‑148 (MQ=255)
                                         ttCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCg                                                                                                    <  2:1521870/148‑1 (MQ=255)
                                          tCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGa                                                                                                   <  1:2329900/148‑1 (MQ=255)
                                          tCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGa                                                                                                   >  2:456054/1‑148 (MQ=255)
                                           cTTGGTAGCATCTAAAACCTTCACCAATCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAg                                                                                                  <  1:1180167/148‑1 (MQ=255)
                                            tttgTAGCATCTAACACCTTCACCACTACGTAAACTATGACCAACGCCCATAGCGCGCGGGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGt                                                                                                 <  1:72466/145‑1 (MQ=255)
                                            ttGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGt                                                                                                 <  1:3116867/148‑1 (MQ=255)
                                             tGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGtt                                                                                                <  1:2716658/148‑1 (MQ=255)
                                                tAGCATCTAAAACCTTCACCTCTCAGGAAACTATGACCAACGCCCATAGCGCGCTTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTgg                                                                                             <  2:2983127/148‑1 (MQ=255)
                                                tAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTgg                                                                                             >  2:18447/1‑148 (MQ=255)
                                                 aGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGt                                                                                            >  1:207934/1‑148 (MQ=255)
                                                 aGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGt                                                                                            >  2:1934032/1‑148 (MQ=255)
                                                   cATCTAAAACCTTCACCACTCATGAAACTATGACCAACGCCCATAGCGCGCGGGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTAt                                                                                          <  1:2879578/148‑1 (MQ=255)
                                                    aTCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTAtt                                                                                         >  1:3737255/1‑148 (MQ=255)
                                                        aaaaCCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATa                                                                                     >  1:2159673/1‑148 (MQ=255)
                                                        aaaaCCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATa                                                                                     <  1:779385/148‑1 (MQ=255)
                                                         aaaCCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTGGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATAc                                                                                    >  2:1599711/1‑148 (MQ=255)
                                                         aaaCCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATAc                                                                                    >  1:1981855/1‑148 (MQ=255)
                                                          aaCCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACt                                                                                   <  2:1185988/148‑1 (MQ=255)
                                                           aCCTTCACCACTCAGGAAACTATGACCAACGCCCATCGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTg                                                                                  >  2:3589680/1‑148 (MQ=255)
                                                           aCCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTg                                                                                  >  2:2708293/1‑148 (MQ=255)
                                                               tCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCaa                                                                              >  1:3183990/1‑148 (MQ=255)
                                                                caccacTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAAc                                                                             <  2:821203/148‑1 (MQ=255)
                                                                   cacTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATg                                                                          >  2:3781993/1‑148 (MQ=255)
                                                                    acTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGt                                                                         <  2:1317012/148‑1 (MQ=255)
                                                                    acTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGt                                                                         <  1:3661582/148‑1 (MQ=255)
                                                                       cAGGAAACTAGGACCCACGCCCATCGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCg                                                                      <  1:2606661/148‑1 (MQ=255)
                                                                              cTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTg                                                               >  1:2762062/1‑148 (MQ=255)
                                                                                aTGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTggg                                                             >  1:2421626/1‑148 (MQ=255)
                                                                                aTGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTggg                                                             >  1:500220/1‑148 (MQ=255)
                                                                                  gACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCTAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGAc                                                           <  2:1107182/148‑1 (MQ=255)
                                                                                   aCCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGTTTATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACt                                                          <  2:2287883/148‑1 (MQ=255)
                                                                                   aCCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACt                                                          >  1:950026/1‑148 (MQ=255)
                                                                                   aCCAACGCCCAGAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGAGACTGCCAACATGTTCGAGTTCTGGGACt                                                          <  2:170350/148‑1 (MQ=255)
                                                                                    ccAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTg                                                         <  2:665164/148‑1 (MQ=255)
                                                                                    ccAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTg                                                         >  1:1058967/1‑148 (MQ=255)
                                                                                      aaCGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTggg                                                       <  2:2668040/148‑1 (MQ=255)
                                                                                       aCGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGt                                                      <  1:1253980/148‑1 (MQ=255)
                                                                                        cGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGtt                                                     <  2:2016427/148‑1 (MQ=255)
                                                                                         gCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTg                                                    <  1:117511/148‑1 (MQ=255)
                                                                                          cccATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTgg                                                   <  2:1480941/148‑1 (MQ=255)
                                                                                            cATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTggcg                                                 >  1:2112548/1‑148 (MQ=255)
                                                                                            cATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTggcg                                                 >  1:1834370/1‑148 (MQ=255)
                                                                                              tAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTggcggc                                               >  1:3486832/1‑148 (MQ=255)
                                                                                              tAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTggcggc                                               >  1:2584069/1‑148 (MQ=255)
                                                                                                gcgcgcGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGCATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCg                                             <  2:42294/148‑1 (MQ=255)
                                                                                                gcgcgcGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCg                                             <  1:1337941/148‑1 (MQ=255)
                                                                                                gcgcgcGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCg                                             <  2:1283645/148‑1 (MQ=255)
                                                                                                 cgcgcgTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGt                                            <  2:2025161/148‑1 (MQ=255)
                                                                                                   cgcgTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTa                                          <  2:2785936/148‑1 (MQ=255)
                                                                                                   cgcgTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTa                                          <  1:2362044/148‑1 (MQ=255)
                                                                                                   cgcgTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTa                                          <  2:2531451/148‑1 (MQ=255)
                                                                                                      gTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTActc                                       <  1:3806710/148‑1 (MQ=255)
                                                                                                       tGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTActct                                      >  2:2961431/1‑148 (MQ=255)
                                                                                                         aCTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCttt                                    <  2:3603644/148‑1 (MQ=255)
                                                                                                          cTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTtg                                   <  2:2153750/148‑1 (MQ=255)
                                                                                                            ggTTCCTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTtgtg                                 >  2:587942/1‑148 (MQ=255)
                                                                                                            ggTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTtgtg                                 <  1:2027669/148‑1 (MQ=255)
                                                                                                                ccTGAAAGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCa                             <  1:2046675/148‑1 (MQ=255)
                                                                                                                    aaaGCGGCAGGTGATGAAAAACACGGTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGcta                         <  2:785149/148‑3 (MQ=255)
                                                                                                                     aaGCGGCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGAt                        <  2:2393203/148‑1 (MQ=255)
                                                                                                                     aaGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCACTGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGAt                        <  2:3488897/148‑1 (MQ=255)
                                                                                                                      cgCGGCAGGTGATGAAACACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGAtt                       <  2:3260916/147‑1 (MQ=255)
                                                                                                                      aGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGAtt                       <  2:3109328/148‑1 (MQ=255)
                                                                                                                         ggCAGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGc                    <  1:3589680/148‑1 (MQ=255)
                                                                                                                            aGGTGATGAAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTg                 <  2:1744822/148‑1 (MQ=255)
                                                                                                                            aGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTg                 <  2:1353232/148‑1 (MQ=255)
                                                                                                                              gTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCTATTGGCCTGTc               >  2:291125/1‑148 (MQ=255)
                                                                                                                              gTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTc               >  2:1879520/1‑148 (MQ=255)
                                                                                                                                gatgaAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTCTTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGCTTGGCCTGTCgc             >  1:2809927/1‑147 (MQ=255)
                                                                                                                                gatgaAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGa             <  1:1142160/148‑1 (MQ=255)
                                                                                                                                 atgaAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGAt            >  1:1924429/1‑148 (MQ=255)
                                                                                                                                 atgaAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGAt            <  2:1087607/148‑1 (MQ=255)
                                                                                                                                  tgaAAAACACTTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGCtt           >  2:518987/1‑148 (MQ=255)
                                                                                                                                      aaaCACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACGCTTTGTGGTGAGCTATTTGCCAGTCGATttttg       >  1:1423388/1‑147 (MQ=255)
                                                                                                                                           cGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTcc  >  2:411423/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
TCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTCACCACTCAGGAAACTATGACCAACGCCCATAGCGCGCGTGACTGGTTCCTGAAAGCGGCAGGTGATGAAAAACACGTTGCAAAACACTTTGCGGCGCTTTCCACCAATGCCAAAGCCGTTGGCGAGTTTGGTATTGATACTGCCAACATGTTCGAGTTCTGGGACTGGGTTGGCGGCCGTTACTCTTTGTGGTCAGCGATTGGCCTGTCGATTGTTCTCTCC  >  CP009273/4224251‑4224537

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: