Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,590,906 C→A 100% A5A (GCG→GCT livK ← leucine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,590,9060CA96.9% 787.8 / 13.9 255A5A (GCG→GCTlivKleucine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base A (105/142);  total (107/148)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 255 total aligned reads displayed.

TTGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAG  >  CP009273/3590764‑3591052
                                                                                                                                              |                                                                                                                                                  
ttGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCAtt                                                                                                                                               >  2:3024889/1‑148 (MQ=255)
    cTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGt                                                                                                                                           <  2:1016417/148‑1 (MQ=255)
     tGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGtt                                                                                                                                          >  2:275720/1‑148 (MQ=255)
      gACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGGGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGttt                                                                                                                                         <  1:1461700/148‑1 (MQ=255)
       aCGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTc                                                                                                                                        >  2:1601499/1‑148 (MQ=255)
       aCGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTc                                                                                                                                        <  2:1329109/148‑1 (MQ=255)
       aCGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTACAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTc                                                                                                                                        >  2:2891595/1‑148 (MQ=255)
          cgcgcCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATa                                                                                                                                     >  2:2026160/1‑148 (MQ=255)
           gcgcCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAg                                                                                                                                    <  1:1471141/148‑1 (MQ=255)
            cgcCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGt                                                                                                                                   >  2:1987010/1‑148 (MQ=255)
              ccGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCa                                                                                                                                 >  1:1729328/1‑148 (MQ=255)
              ccGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCa                                                                                                                                 >  1:828286/1‑148 (MQ=255)
                 ttAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaaa                                                                                                                              >  1:2524106/1‑148 (MQ=255)
                  tAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaaaa                                                                                                                             <  2:678447/148‑1 (MQ=255)
                   aaaTTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAAt                                                                                                                            <  1:2634152/148‑1 (MQ=255)
                     aTTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATcc                                                                                                                          >  1:1000561/1‑148 (MQ=255)
                       tCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATcccc                                                                                                                        <  1:2697918/148‑1 (MQ=255)
                        ccATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCa                                                                                                                       <  1:500659/148‑1 (MQ=255)
                            atCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCg                                                                                                                   <  2:3377020/148‑1 (MQ=255)
                            atCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCg                                                                                                                   >  1:1475524/1‑148 (MQ=255)
                                 cccACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtg                                                                                                              <  2:2582308/148‑1 (MQ=255)
                                 cccACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtg                                                                                                              <  2:840131/148‑1 (MQ=255)
                                    aCTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttg                                                                                                           <  1:2820543/148‑1 (MQ=255)
                                     cTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgt                                                                                                          >  1:1079706/1‑148 (MQ=255)
                                       gggCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgt                                                                                                        <  2:3125267/148‑1 (MQ=255)
                                       gggCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgt                                                                                                        <  2:3088486/148‑1 (MQ=255)
                                           aaTCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCt                                                                                                    <  1:1934760/148‑1 (MQ=255)
                                                ggCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGtt                                                                                               >  1:3171033/1‑148 (MQ=255)
                                                  ccGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGtttt                                                                                             <  2:2795849/148‑1 (MQ=255)
                                                    ggACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTa                                                                                           <  2:3533664/148‑1 (MQ=255)
                                                    ggACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTa                                                                                           <  1:3224724/148‑1 (MQ=255)
                                                     gACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTAt                                                                                          >  2:2284323/1‑148 (MQ=255)
                                                      aCATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTAtg                                                                                         >  1:1520624/1‑148 (MQ=255)
                                                       cATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTAtgt                                                                                        <  2:2403079/148‑1 (MQ=255)
                                                         tCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATgtgt                                                                                      <  2:2727600/148‑1 (MQ=255)
                                                          cgcgCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGtt                                                                                     >  1:613093/1‑148 (MQ=255)
                                                               cGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACaa                                                                                <  1:214944/148‑1 (MQ=255)
                                                               cGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACaa                                                                                <  2:1193010/148‑1 (MQ=255)
                                                               cGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACaa                                                                                <  2:1589346/148‑1 (MQ=255)
                                                               cGACAACGGCGACTTTAAGATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGGTGTGTTGCTTTGTTTTTATGTGTTAACaa                                                                                <  2:2911816/148‑1 (MQ=255)
                                                                     cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                          >  2:2229337/1‑148 (MQ=255)
                                                                     cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                          <  2:1146948/148‑1 (MQ=255)
                                                                      ggCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGAc                                                                         <  2:1729328/148‑1 (MQ=255)
                                                                        cGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTg                                                                       <  1:953698/148‑1 (MQ=255)
                                                                          aCTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGtt                                                                     <  1:2705207/148‑1 (MQ=255)
                                                                          aCTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGtt                                                                     <  2:44254/148‑1 (MQ=255)
                                                                             ttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtt                                                                  <  2:2628042/148‑1 (MQ=255)
                                                                             ttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtt                                                                  <  1:153735/148‑1 (MQ=255)
                                                                             ttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtt                                                                  >  2:2195569/1‑148 (MQ=255)
                                                                             ttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtt                                                                  <  1:2174312/148‑1 (MQ=255)
                                                                             ttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtt                                                                  <  2:1968785/148‑1 (MQ=255)
                                                                               aaTATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtttt                                                                <  1:3422769/148‑1 (MQ=255)
                                                                                 tatCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtttttt                                                              <  1:737345/148‑1 (MQ=255)
                                                                                 tatCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtttttt                                                              >  1:2261558/1‑148 (MQ=255)
                                                                                  atCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttttttt                                                             <  1:1563491/148‑1 (MQ=255)
                                                                                  atCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttttttt                                                             <  1:3325697/148‑1 (MQ=255)
                                                                                   tCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTta                                                            >  1:779393/1‑148 (MQ=255)
                                                                                   tCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTta                                                            >  1:678952/1‑148 (MQ=255)
                                                                                      tCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATAc                                                         <  1:1137575/148‑1 (MQ=255)
                                                                                      tCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATAc                                                         >  1:2945447/1‑148 (MQ=255)
                                                                                      tCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATAc                                                         >  2:1818663/1‑148 (MQ=255)
                                                                                       cAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACt                                                        >  2:2964211/1‑148 (MQ=255)
                                                                                         gCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGc                                                      <  2:1527231/148‑1 (MQ=255)
                                                                                         gCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGc                                                      >  2:3488697/1‑148 (MQ=255)
                                                                                         gCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGc                                                      >  2:2098963/1‑148 (MQ=255)
                                                                                          ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                     <  2:2318615/148‑1 (MQ=255)
                                                                                          ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                     >  2:2111022/1‑148 (MQ=255)
                                                                                          ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                     <  1:1600616/148‑1 (MQ=255)
                                                                                           cATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCAc                                                    >  2:2441852/1‑148 (MQ=255)
                                                                                           cATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCAc                                                    >  2:1281194/1‑148 (MQ=255)
                                                                                            aTAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACt                                                   >  2:910257/1‑148 (MQ=255)
                                                                                             tAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTg                                                  >  1:3483150/1‑148 (MQ=255)
                                                                                                 ggTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGtttt                                              <  2:1950151/148‑1 (MQ=255)
                                                                                                 ggTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGtttt                                              <  1:2188541/148‑1 (MQ=255)
                                                                                                   tgtgtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTg                                            >  2:129326/1‑148 (MQ=255)
                                                                                                      gtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCt                                         <  1:2104680/148‑1 (MQ=255)
                                                                                                       tgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTg                                        <  1:1987010/148‑1 (MQ=255)
                                                                                                         aaaTTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTc                                      >  1:1089250/1‑148 (MQ=255)
                                                                                                            ttGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGa                                   <  1:916527/148‑1 (MQ=255)
                                                                                                             tGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGAt                                  >  1:3149737/1‑148 (MQ=255)
                                                                                                               ccAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTa                                >  1:1262884/1‑148 (MQ=255)
                                                                                                                 aGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAg                              <  1:2577903/148‑1 (MQ=255)
                                                                                                                   tGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAggg                            >  2:2628591/1‑148 (MQ=255)
                                                                                                                   tGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGAGTAAggg                            >  1:3119160/1‑148 (MQ=255)
                                                                                                                    gCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAgggg                           <  1:2284323/148‑1 (MQ=255)
                                                                                                                         catcCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAgc                      >  1:200115/1‑148 (MQ=255)
                                                                                                                             ccTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAg                  >  1:1218898/1‑148 (MQ=255)
                                                                                                                                gCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAAGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTTCACTGTTTTTGACTGTCTGATTAAGGGGTTAGCTCAGTAt               >  1:2305402/1‑148 (MQ=255)
                                                                                                                                  gatgatAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAttt             <  1:1232596/148‑1 (MQ=255)
                                                                                                                                    tgatAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCTTGATGTTTTGTTGATTTGTTTTTATGTGTTAACAAATCAGACTGTTCTGTTTTTATACTTCACTTTTTTTTAATGTATGATTAAGGGGTTAGCGCAGTATTTtg           >  1:367947/1‑148 (MQ=255)
                                                                                                                                    tgatAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTtg           <  2:1079706/148‑1 (MQ=255)
                                                                                                                                     gatAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTtgt          <  1:1049611/148‑1 (MQ=255)
                                                                                                                                     gatAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTtgt          <  1:2441852/148‑1 (MQ=255)
                                                                                                                                      atAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCCTTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTtgtg         <  2:2514346/148‑1 (MQ=255)
                                                                                                                                       tAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGa        >  2:1538007/1‑148 (MQ=255)
                                                                                                                                       tAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGa        <  2:1093890/148‑1 (MQ=255)
                                                                                                                                       tAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGa        >  1:1654801/1‑148 (MQ=255)
                                                                                                                                        aGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTACTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGat       >  2:3329785/1‑148 (MQ=255)
                                                                                                                                           tttAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataat    <  1:2783694/148‑1 (MQ=255)
                                                                                                                                           tttAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataat    <  2:948866/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TTGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAG  >  CP009273/3590764‑3591052

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: