Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 943,573 G→T 100% A297S (GCC→TCC)  ycaM → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273943,5730GT97.5% 808.1 / 9.7 236A297S (GCC→TCC) ycaMputative transporter
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (106/124);  total (107/129)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 236 total aligned reads displayed.

GTTTGTGCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATGGCTAT  >  CP009273/943426‑943719
                                                                                                                                                   |                                                                                                                                                  
gTTTGTGCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTcc                                                                                                                                                  >  2:1913529/1‑150 (MQ=255)
  ttGTGCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCAt                                                                                                                                                >  1:891189/1‑150 (MQ=255)
     tgCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGc                                                                                                                                             <  2:4557/150‑1 (MQ=255)
       ccATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGa                                                                                                                                           <  2:2794963/150‑1 (MQ=255)
        cATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGaa                                                                                                                                          >  2:4049710/1‑150 (MQ=255)
         aTTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAAt                                                                                                                                         >  1:1865656/1‑150 (MQ=255)
         aTTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAAt                                                                                                                                         <  2:3007138/150‑1 (MQ=255)
         aTTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAAt                                                                                                                                         <  2:3230423/150‑1 (MQ=255)
         aTTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAAt                                                                                                                                         <  2:3230657/150‑1 (MQ=255)
            cTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATAcc                                                                                                                                      >  2:486111/1‑150 (MQ=255)
             tGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATAccc                                                                                                                                     >  1:268621/1‑150 (MQ=255)
             tGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATAccc                                                                                                                                     >  1:1488345/1‑150 (MQ=255)
              gggCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCt                                                                                                                                    >  2:2523635/1‑150 (MQ=255)
              gggCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCt                                                                                                                                    >  2:2523627/1‑150 (MQ=255)
                  tCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGa                                                                                                                                >  2:4108666/1‑150 (MQ=255)
                    gCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACa                                                                                                                              <  2:2399392/150‑1 (MQ=255)
                       ggCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGt                                                                                                                           <  2:1664694/150‑1 (MQ=255)
                        gCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTa                                                                                                                          <  1:1116074/150‑1 (MQ=255)
                        gCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTa                                                                                                                          <  1:537679/150‑1 (MQ=255)
                           aTGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTATTGACCTTCGGTCAGTATTTCGCCTTTCAGAAGCTGGGCGAGTTTTACACCATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATTCCCTGGGACAAGTAgcg                                                                                                                       <  2:2657829/150‑1 (MQ=255)
                           aTGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcg                                                                                                                       >  1:323987/1‑150 (MQ=255)
                            tGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcgg                                                                                                                      <  2:1526708/150‑1 (MQ=255)
                             ggggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggc                                                                                                                     <  1:3373927/150‑1 (MQ=255)
                             ggggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggc                                                                                                                     <  1:868844/150‑1 (MQ=255)
                             ggggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggc                                                                                                                     <  2:1162902/150‑1 (MQ=255)
                             ggggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggc                                                                                                                     <  1:438398/150‑1 (MQ=255)
                              gggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggcg                                                                                                                    <  2:1853903/150‑1 (MQ=255)
                              gggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggcg                                                                                                                    <  2:3952966/150‑1 (MQ=255)
                              gggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAgcggcg                                                                                                                    <  2:860040/150‑1 (MQ=255)
                               ggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGgcgc                                                                                                                   >  1:2620775/1‑150 (MQ=255)
                                 atgatgTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTg                                                                                                                 <  2:1073436/150‑1 (MQ=255)
                                  tgttgTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTgg                                                                                                                <  1:3049313/147‑1 (MQ=255)
                                  tgatgTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTgg                                                                                                                <  1:2447250/150‑1 (MQ=255)
                                     tgTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTAt                                                                                                             <  1:3354819/150‑1 (MQ=255)
                                       tttGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATtc                                                                                                           <  2:234925/150‑1 (MQ=255)
                                       tttGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATtc                                                                                                           >  2:1413982/1‑150 (MQ=255)
                                        ttGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATtct                                                                                                          <  1:3800203/150‑1 (MQ=255)
                                        ttGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATtct                                                                                                          <  1:1701254/150‑1 (MQ=255)
                                         tGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATtctc                                                                                                         <  2:446691/150‑1 (MQ=255)
                                          gATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCg                                                                                                        >  1:2785605/1‑150 (MQ=255)
                                                 gTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATg                                                                                                 >  1:1181276/1‑150 (MQ=255)
                                                   aaTATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGcc                                                                                               >  1:402515/1‑150 (MQ=255)
                                                    atatCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGccc                                                                                              >  1:1952301/1‑150 (MQ=255)
                                                     tatCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGcccc                                                                                             <  1:1816508/150‑1 (MQ=255)
                                                       tCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGc                                                                                           >  1:2461274/1‑150 (MQ=255)
                                                        cccGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCt                                                                                          >  2:363974/1‑150 (MQ=255)
                                                         ccGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCtt                                                                                         >  1:667165/1‑150 (MQ=255)
                                                               gACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTg                                                                                   >  1:397100/1‑150 (MQ=255)
                                                               gACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTg                                                                                   >  1:583152/1‑150 (MQ=255)
                                                                aCTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGc                                                                                  >  1:220519/1‑150 (MQ=255)
                                                                aCTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGc                                                                                  >  1:988951/1‑150 (MQ=255)
                                                                 cTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCt                                                                                 >  2:2104719/1‑150 (MQ=255)
                                                                   tAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTAt                                                                               >  1:2373275/1‑150 (MQ=255)
                                                                      tGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAg                                                                            <  1:2081738/150‑1 (MQ=255)
                                                                       gACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAgg                                                                           <  1:1105145/150‑1 (MQ=255)
                                                                           aaCGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGAt                                                                       >  2:58227/1‑150 (MQ=255)
                                                                           aaCGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGAt                                                                       <  1:1260638/150‑1 (MQ=255)
                                                                           aaCGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGAt                                                                       >  1:1776696/1‑150 (MQ=255)
                                                                           aaCGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGAt                                                                       <  1:1904653/150‑1 (MQ=255)
                                                                            aCGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATg                                                                      <  1:1664438/150‑1 (MQ=255)
                                                                            aCGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATg                                                                      >  1:3233654/1‑150 (MQ=255)
                                                                                tCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGa                                                                  >  1:981295/1‑150 (MQ=255)
                                                                                 cAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGAc                                                                 <  2:509097/150‑1 (MQ=255)
                                                                                 cAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGAc                                                                 <  2:509001/150‑1 (MQ=255)
                                                                                 cAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGAc                                                                 <  1:3483681/150‑1 (MQ=255)
                                                                                  aGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACa                                                                <  2:617945/150‑1 (MQ=255)
                                                                                           ccTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAatata                                                       <  2:881693/150‑1 (MQ=255)
                                                                                           ccTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAatata                                                       >  2:695088/1‑150 (MQ=255)
                                                                                           ccTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAatata                                                       >  2:694818/1‑150 (MQ=255)
                                                                                             tttCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATAtt                                                     <  1:3779942/150‑1 (MQ=255)
                                                                                              ttCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTc                                                    <  1:10684/150‑1 (MQ=255)
                                                                                              ttCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTc                                                    <  1:1744610/150‑1 (MQ=255)
                                                                                                cAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTcca                                                  <  1:906069/150‑1 (MQ=255)
                                                                                                 aGAAGCTGGGCGAGTATTACAACTTGGTTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTccag                                                 <  1:3943572/150‑1 (MQ=255)
                                                                                                 aGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTccag                                                 >  1:2971479/1‑150 (MQ=255)
                                                                                                 aGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTccag                                                 >  1:3625054/1‑150 (MQ=255)
                                                                                                  gAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTccagc                                                <  2:3035739/150‑1 (MQ=255)
                                                                                                   aaGCTGTGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTccagcc                                               <  2:3777264/150‑1 (MQ=255)
                                                                                                       tGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGtt                                           >  2:2224211/1‑150 (MQ=255)
                                                                                                       tGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGtt                                           >  2:3959940/1‑150 (MQ=255)
                                                                                                        gggCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGttt                                          <  1:4082316/150‑1 (MQ=255)
                                                                                                        gggCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGttt                                          <  2:3092267/150‑1 (MQ=255)
                                                                                                           cGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATg                                       <  1:403126/150‑1 (MQ=255)
                                                                                                                aTTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTa                                  >  1:1413850/1‑150 (MQ=255)
                                                                                                                   acaacaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCa                               >  2:3581503/1‑150 (MQ=255)
                                                                                                                    caacaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCaa                              >  2:2468653/1‑150 (MQ=255)
                                                                                                                      acaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACg                            >  1:318177/1‑150 (MQ=255)
                                                                                                                       caTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGc                           >  1:2674287/1‑150 (MQ=255)
                                                                                                                          gggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTc                        <  1:2099230/150‑1 (MQ=255)
                                                                                                                          gggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTc                        >  2:3692828/1‑150 (MQ=255)
                                                                                                                             tAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGTGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTgg                     <  1:2628531/150‑1 (MQ=255)
                                                                                                                                tACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTAc                  <  1:1691501/150‑1 (MQ=255)
                                                                                                                                 aCTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACg                 >  1:3385943/1‑150 (MQ=255)
                                                                                                                                 aCTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACg                 >  2:217475/1‑150 (MQ=255)
                                                                                                                                      aaTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGt            >  1:1537937/1‑150 (MQ=255)
                                                                                                                                      aaTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGt            >  1:3635183/1‑150 (MQ=255)
                                                                                                                                      aaTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGt            >  1:302612/1‑150 (MQ=255)
                                                                                                                                      aaTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGt            >  2:319732/1‑150 (MQ=255)
                                                                                                                                        tGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTa          <  1:3804905/150‑1 (MQ=255)
                                                                                                                                          gTGATTTACTCCATTGCGAATACCCTGGTACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAAt        <  1:3450111/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTTGTGCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATGGCTAT  >  CP009273/943426‑943719

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: