Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,890,151 G→T 17.9% A6A (GCG→GCT yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,890,1510GT17.9% 654.3 / 85.7 233A6A (GCG→GCTyieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (87/104);  new base T (20/22);  total (107/126)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 237 total aligned reads displayed.

CTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAATTAT  >  CP009273/3890004‑3890297
                                                                                                                                                   |                                                                                                                                                  
cTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCg                                                                                                                                                    >  2:1630867/1‑148 (MQ=255)
   tCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTa                                                                                                                                                 >  1:2034244/1‑148 (MQ=255)
   tCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTa                                                                                                                                                 <  2:2019479/148‑1 (MQ=255)
   tCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTa                                                                                                                                                 >  1:3087171/1‑148 (MQ=255)
      tataCCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTAttt                                                                                                                                              <  2:811895/148‑1 (MQ=255)
        taCCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTAttttt                                                                                                                                            <  2:26160/148‑1 (MQ=255)
          ccAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCg                                                                                                                                          <  2:1054200/148‑1 (MQ=255)
          ccAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCg                                                                                                                                          <  2:3551728/148‑1 (MQ=255)
           cAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGa                                                                                                                                         <  1:1550995/148‑1 (MQ=255)
           cAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGa                                                                                                                                         <  1:2767771/148‑1 (MQ=255)
            aaaaTGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGAc                                                                                                                                        <  2:3436087/148‑1 (MQ=255)
            aaaaTGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGAc                                                                                                                                        <  2:1391594/148‑1 (MQ=255)
              aaTGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTg                                                                                                                                      >  2:3288865/1‑148 (MQ=255)
               aTGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGc                                                                                                                                     >  1:2969299/1‑148 (MQ=255)
                 gCGACACTAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGa                                                                                                                                   >  2:2034865/1‑148 (MQ=255)
                     cAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGt                                                                                                                               >  2:915447/1‑148 (MQ=255)
                       aTAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTAc                                                                                                                             >  1:2619470/1‑148 (MQ=255)
                            cgcgTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTgg                                                                                                                        >  1:1005088/1‑148 (MQ=255)
                             gcgTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGt                                                                                                                       <  2:455539/148‑1 (MQ=255)
                             gcgTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGt                                                                                                                       <  1:427732/148‑1 (MQ=255)
                                 ttGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGAGACGGTACGCTGGTCGAc                                                                                                                   >  2:2239523/1‑148 (MQ=255)
                                   gtgAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAg                                                                                                                 <  2:1353000/148‑1 (MQ=255)
                                   gtgAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAg                                                                                                                 <  1:3499795/148‑1 (MQ=255)
                                   gtgAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAg                                                                                                                 <  2:2877203/148‑1 (MQ=255)
                                       gagaCTTTTTTATTGGATTTAATTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGaa                                                                                                             <  1:1606828/148‑1 (MQ=255)
                                       gagaCTTTTTTATGGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGaa                                                                                                             <  1:1606704/148‑1 (MQ=255)
                                        agaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAg                                                                                                            <  2:243695/148‑1 (MQ=255)
                                        agaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAg                                                                                                            >  1:3302842/1‑148 (MQ=255)
                                        agaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAg                                                                                                            >  2:1689674/1‑148 (MQ=255)
                                        agaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAg                                                                                                            <  2:1305693/148‑1 (MQ=255)
                                          aCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTc                                                                                                          <  2:891367/148‑1 (MQ=255)
                                          aCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTc                                                                                                          <  2:20884/148‑1 (MQ=255)
                                            ttttttATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCAt                                                                                                        >  2:1671345/1‑148 (MQ=255)
                                              ttttATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCAttt                                                                                                      <  2:1172332/148‑1 (MQ=255)
                                                ttATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGc                                                                                                    <  1:1326636/148‑1 (MQ=255)
                                                   ttGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGACGTCATTTGctct                                                                                                 >  2:904785/1‑148 (MQ=255)
                                                          tACTTATACGGCAACGTTTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAATGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCat                                                                                          <  1:2397175/148‑1 (MQ=255)
                                                           aCTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCata                                                                                         >  1:3154622/1‑148 (MQ=255)
                                                              tataCGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGt                                                                                      >  2:191617/1‑148 (MQ=255)
                                                              tataCGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGt                                                                                      >  1:1021385/1‑148 (MQ=255)
                                                                  cGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACg                                                                                  >  1:3389674/1‑148 (MQ=255)
                                                                  cGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACg                                                                                  <  1:1630867/148‑1 (MQ=255)
                                                                   ggCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGa                                                                                 >  1:2882522/1‑148 (MQ=255)
                                                                    gCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGAt                                                                                <  2:3275243/148‑1 (MQ=255)
                                                                    gCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGAt                                                                                <  2:2737196/148‑1 (MQ=255)
                                                                    gCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGAt                                                                                >  2:3402803/1‑148 (MQ=255)
                                                                     cAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATg                                                                               >  2:2403776/1‑148 (MQ=255)
                                                                     cAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATg                                                                               >  2:3064164/1‑148 (MQ=255)
                                                                        cGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGttt                                                                            >  1:248629/1‑148 (MQ=255)
                                                                        cGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGttt                                                                            >  1:1532406/1‑148 (MQ=255)
                                                                                ttgcgCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAAtt                                                                    >  1:227341/1‑148 (MQ=255)
                                                                                 tgcgCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAAttt                                                                   <  1:1957294/148‑1 (MQ=255)
                                                                                   cgcAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTgg                                                                 <  1:2866695/148‑1 (MQ=255)
                                                                                    gcAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGt                                                                <  2:1284068/148‑1 (MQ=255)
                                                                                      aaaTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTAt                                                              <  1:1817917/148‑1 (MQ=255)
                                                                                        aTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTa                                                            <  1:1671345/148‑1 (MQ=255)
                                                                                          cGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACg                                                          <  2:3466945/148‑1 (MQ=255)
                                                                                           ggCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGc                                                         <  1:2403808/148‑1 (MQ=255)
                                                                                                 aCGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATc                                                   <  1:3288865/148‑1 (MQ=255)
                                                                                                 aCGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTGGGTATTACGCTCGATc                                                   <  2:56468/148‑1 (MQ=255)
                                                                                                  cGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATcc                                                  <  1:1436252/148‑1 (MQ=255)
                                                                                                   gAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCt                                                 <  2:1170756/148‑1 (MQ=255)
                                                                                                   gAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCt                                                 <  2:2948865/148‑1 (MQ=255)
                                                                                                     aaaCGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCATGa                                               <  2:2882522/148‑1 (MQ=255)
                                                                                                       aCGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAag                                             >  2:2766901/1‑148 (MQ=255)
                                                                                                       aCGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAag                                             <  1:915447/148‑1 (MQ=255)
                                                                                                         gTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAagag                                           <  2:1533822/148‑1 (MQ=255)
                                                                                                          ttAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAgg                                          >  1:784687/1‑148 (MQ=255)
                                                                                                          ttAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAgg                                          >  2:679651/1‑148 (MQ=255)
                                                                                                             aaCTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTAt                                       >  1:2214097/1‑148 (MQ=255)
                                                                                                              aCTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTAtt                                      >  1:134159/1‑148 (MQ=255)
                                                                                                              aCTGATGAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTAtt                                      <  2:2458656/148‑1 (MQ=255)
                                                                                                                 gATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCaa                                   >  1:1374385/1‑148 (MQ=255)
                                                                                                                  aTTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCaaa                                  <  2:833531/148‑1 (MQ=255)
                                                                                                                  aTTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCaaa                                  >  1:1079995/1‑148 (MQ=255)
                                                                                                                    tAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACg                                <  2:181712/148‑1 (MQ=255)
                                                                                                                     aaaaaGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGt                               <  1:630393/148‑1 (MQ=255)
                                                                                                                     aaaaaGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGt                               >  1:2559322/1‑148 (MQ=255)
                                                                                                                      aaaaGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGAGGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGtt                              >  2:2121301/1‑148 (MQ=255)
                                                                                                                       aaaggaaaggCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGttt                             <  2:2105691/148‑1 (MQ=255)
                                                                                                                       aaaggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGttt                             <  2:1896511/148‑1 (MQ=255)
                                                                                                                       aaaggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGttt                             >  1:1724919/1‑148 (MQ=255)
                                                                                                                       aaaggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGttt                             <  2:2076602/148‑1 (MQ=255)
                                                                                                                          ggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCaa                          >  1:3249684/1‑148 (MQ=255)
                                                                                                                           gaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCaaa                         <  1:3553085/148‑1 (MQ=255)
                                                                                                                            aaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAg                        >  1:2160210/1‑148 (MQ=255)
                                                                                                                             aaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGAATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAgg                       <  2:285351/148‑1 (MQ=255)
                                                                                                                               ggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGgtg                     >  1:3378263/1‑148 (MQ=255)
                                                                                                                                      cccGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACt              <  2:425946/148‑1 (MQ=255)
                                                                                                                                        cGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGt            >  1:3453165/1‑148 (MQ=255)
                                                                                                                                        cGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGt            <  2:3324687/148‑1 (MQ=255)
                                                                                                                                        cGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGt            <  1:3118890/148‑1 (MQ=255)
                                                                                                                                         ggATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTa           <  1:3181478/148‑1 (MQ=255)
                                                                                                                                         ggATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTa           <  2:1251194/148‑1 (MQ=255)
                                                                                                                                           aTAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACg         >  1:973402/1‑148 (MQ=255)
                                                                                                                                           aTAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACg         >  1:79725/1‑148 (MQ=255)
                                                                                                                                           aTAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACg         >  1:72647/1‑148 (MQ=255)
                                                                                                                                            tAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGa        <  1:191617/148‑1 (MQ=255)
                                                                                                                                             aGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGaa       >  1:717662/1‑148 (MQ=255)
                                                                                                                                                 gCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAATTa   <  1:536740/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAATTAT  >  CP009273/3890004‑3890297

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: