Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,242,509 C→T 100% L26F (CTC→TTC)  ubiC → chorismate‑‑pyruvate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,242,5090CT96.7% 786.9 / 12.3 244L26F (CTC→TTC) ubiCchorismate‑‑pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base T (110/125);  total (114/130)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 244 total aligned reads displayed.

TTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTT  >  CP009273/4242364‑4242653
                                                                                                                                                 |                                                                                                                                                
ttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTc                                                                                                                                                >  1:1614780/1‑148 (MQ=255)
  gCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGa                                                                                                                                              >  1:2058616/1‑148 (MQ=255)
  gCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGa                                                                                                                                              <  2:1744759/148‑1 (MQ=255)
   cTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGAc                                                                                                                                             >  2:954546/1‑148 (MQ=255)
    tttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACt                                                                                                                                            >  2:1403351/1‑148 (MQ=255)
    tttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACt                                                                                                                                            <  1:998668/148‑1 (MQ=255)
    tttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACt                                                                                                                                            >  2:2009874/1‑148 (MQ=255)
     ttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                           >  1:1001/1‑148 (MQ=255)
     ttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGCAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTCCGAGCGATGAGCTATAGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                           >  1:732/1‑148 (MQ=255)
      tACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTgg                                                                                                                                          <  1:2342239/148‑1 (MQ=255)
      tACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTgg                                                                                                                                          <  2:2737631/148‑1 (MQ=255)
      tACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTgg                                                                                                                                          <  1:2978450/148‑1 (MQ=255)
          tctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGt                                                                                                                                      >  1:1110413/1‑148 (MQ=255)
          tctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGt                                                                                                                                      >  1:3259404/1‑148 (MQ=255)
            tcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTg                                                                                                                                    <  2:1643421/148‑1 (MQ=255)
             cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGc                                                                                                                                   >  2:243776/1‑148 (MQ=255)
             cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGc                                                                                                                                   <  2:2533294/148‑1 (MQ=255)
             cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGc                                                                                                                                   >  2:3027464/1‑148 (MQ=255)
                  aaaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggag                                                                                                                              >  2:1360894/1‑148 (MQ=255)
                   aaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagg                                                                                                                             >  1:2649845/1‑148 (MQ=255)
                   aaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagg                                                                                                                             <  2:472446/148‑1 (MQ=255)
                      acTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAtt                                                                                                                          <  2:3524886/148‑1 (MQ=255)
                        ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTcc                                                                                                                        <  1:668130/148‑1 (MQ=255)
                        ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTcc                                                                                                                        <  1:1895310/148‑1 (MQ=255)
                         tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                       <  2:2037696/148‑1 (MQ=255)
                         tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                       <  2:2413346/148‑1 (MQ=255)
                          ttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                      <  2:1224339/148‑1 (MQ=255)
                            cTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                    >  1:3295968/1‑148 (MQ=255)
                               cGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaa                                                                                                                 <  2:2826895/148‑1 (MQ=255)
                                gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaa                                                                                                                <  1:1374199/148‑1 (MQ=255)
                                gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaa                                                                                                                <  2:2982753/148‑1 (MQ=255)
                                 aTACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaaa                                                                                                               >  2:26352/1‑148 (MQ=255)
                                 aTACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaaa                                                                                                               >  1:1044487/1‑148 (MQ=255)
                                    cAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                    <  1:708532/140‑1 (MQ=255)
                                      aTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGttt                                                                                                          <  1:3371755/148‑1 (MQ=255)
                                       tGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGAGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGtttt                                                                                                         >  2:471613/1‑148 (MQ=255)
                                        gCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTg                                                                                                        >  2:1825418/1‑148 (MQ=255)
                                        gCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACAACCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTg                                                                                                        <  1:3081496/148‑1 (MQ=255)
                                             tACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcag                                                                                                   <  2:3016842/148‑1 (MQ=255)
                                             tACGGTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcag                                                                                                   <  2:899751/148‑1 (MQ=255)
                                              aCGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagc                                                                                                  <  2:1558273/148‑1 (MQ=255)
                                                 ttATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAgg                                                                                               >  1:2229739/1‑148 (MQ=255)
                                                 ttATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTAGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAgg                                                                                               >  2:749468/1‑148 (MQ=255)
                                                  tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGATGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAggg                                                                                              >  2:957620/1‑148 (MQ=255)
                                                   aTGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGa                                                                                             <  1:1765037/148‑1 (MQ=255)
                                                     gTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaaa                                                                                           >  1:945481/1‑148 (MQ=255)
                                                      tAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaaaa                                                                                          >  2:1184568/1‑148 (MQ=255)
                                                         cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                       >  1:1919206/1‑148 (MQ=255)
                                                         cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTATTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                       >  2:646052/1‑148 (MQ=255)
                                                                 tCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGt                                                                               <  2:3135080/148‑1 (MQ=255)
                                                                      aTGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                          <  2:2651309/148‑1 (MQ=255)
                                                                      aTGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                          <  1:1036819/148‑1 (MQ=255)
                                                                       tGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgat                                                                         <  1:386155/148‑1 (MQ=255)
                                                                        gTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatg                                                                        >  1:339783/1‑148 (MQ=255)
                                                                          cacaCCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatgat                                                                      <  1:2051665/148‑1 (MQ=255)
                                                                          cacaCCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatgat                                                                      <  1:1461384/148‑1 (MQ=255)
                                                                             acCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcg                                                                   >  1:1096132/1‑148 (MQ=255)
                                                                                ccGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGa                                                                >  1:3103717/1‑148 (MQ=255)
                                                                                ccGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGa                                                                >  2:2940171/1‑148 (MQ=255)
                                                                                 cgcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGaa                                                               >  1:1347393/1‑148 (MQ=255)
                                                                                    gTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAggg                                                            <  2:3259404/148‑1 (MQ=255)
                                                                                     ttAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACAGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAgtct                                                           >  2:2889202/1‑145 (MQ=255)
                                                                                        aCGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTg                                                        <  2:1919206/148‑1 (MQ=255)
                                                                                        aCGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTg                                                        <  1:2813578/148‑1 (MQ=255)
                                                                                         cGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGt                                                       >  1:384627/1‑148 (MQ=255)
                                                                                         cGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGt                                                       >  1:1084479/1‑148 (MQ=255)
                                                                                           cAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCg                                                     <  2:1940439/148‑1 (MQ=255)
                                                                                              cTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGc                                                  >  1:2244661/1‑148 (MQ=255)
                                                                                                    gcgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATg                                            <  1:2670960/148‑1 (MQ=255)
                                                                                                      gctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaa                                          <  2:2357522/148‑1 (MQ=255)
                                                                                                      gctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaa                                          <  1:883191/148‑1 (MQ=255)
                                                                                                       ctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                         >  2:3319255/1‑148 (MQ=255)
                                                                                                        tgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAAt                                        <  2:2102891/148‑1 (MQ=255)
                                                                                                         gcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATc                                       <  1:1915693/148‑1 (MQ=255)
                                                                                                              aTTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCga                                  >  2:1753713/1‑148 (MQ=255)
                                                                                                                tGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaag                                >  2:1609678/1‑148 (MQ=255)
                                                                                                                 gTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaaga                               <  1:444931/148‑1 (MQ=255)
                                                                                                                 gTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaaga                               <  1:1999921/148‑1 (MQ=255)
                                                                                                                     agagATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTg                           <  2:520874/148‑1 (MQ=255)
                                                                                                                      gagaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGc                          >  1:3118791/1‑148 (MQ=255)
                                                                                                                      gagaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGc                          >  1:3239114/1‑148 (MQ=255)
                                                                                                                      gagaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGc                          <  2:3196411/148‑1 (MQ=255)
                                                                                                                        gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCg                        <  1:2953080/148‑1 (MQ=255)
                                                                                                                         aTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgc                       <  2:2699761/148‑1 (MQ=255)
                                                                                                                          tCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCANGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                      >  1:1802801/1‑148 (MQ=255)
                                                                                                                           ccctgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgctg                     >  1:1167570/1‑148 (MQ=255)
                                                                                                                            cctgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgctgc                    <  2:2153527/148‑1 (MQ=255)
                                                                                                                                ccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGcc                >  1:59916/1‑148 (MQ=255)
                                                                                                                                   tgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAATCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaa             >  1:2490334/1‑148 (MQ=255)
                                                                                                                                     gATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAag           <  1:1850386/148‑1 (MQ=255)
                                                                                                                                     gATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAag           <  1:2926974/148‑1 (MQ=255)
                                                                                                                                      aTCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAaga          >  1:3554304/1‑148 (MQ=255)
                                                                                                                                        ccGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGt        >  2:3376642/1‑148 (MQ=255)
                                                                                                                                        ccGCAACTGTTCGACTGGCTGGTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGt        <  2:2541227/148‑1 (MQ=255)
                                                                                                                                            aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCg    <  1:3319255/148‑1 (MQ=255)
                                                                                                                                            aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCg    >  1:2365515/1‑148 (MQ=255)
                                                                                                                                            aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCg    <  1:472017/148‑1 (MQ=255)
                                                                                                                                             aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGt   <  1:1939632/148‑1 (MQ=255)
                                                                                                                                              cTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt  >  2:2750838/1‑148 (MQ=255)
                                                                                                                                              cTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt  >  2:2369472/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTT  >  CP009273/4242364‑4242653

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: