Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,169,719 C→A 18.2% E28* (GAA→TAA)  gatZ ← D‑tagatose 1,6‑bisphosphate aldolase 2, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,169,7190CA18.2% 628.3 / 103.5 218E28* (GAA→TAA) gatZD‑tagatose 1,6‑bisphosphate aldolase 2, subunit
Reads supporting (aligned to +/- strand):  ref base C (88/90);  new base A (16/24);  total (104/114)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

AAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACTTTG  >  CP009273/2169570‑2169868
                                                                                                                                                     |                                                                                                                                                     
aaaCCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTc                                                                                                                                                       <  2:4179902/150‑1 (MQ=255)
  aCCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGa                                                                                                                                                     >  2:2415945/1‑150 (MQ=255)
   cccAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGAt                                                                                                                                                    <  2:2750583/150‑1 (MQ=255)
    ccAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATa                                                                                                                                                   <  1:522374/150‑1 (MQ=255)
     cAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATaa                                                                                                                                                  >  1:3570476/1‑150 (MQ=255)
      aaCTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAAc                                                                                                                                                 <  2:1322919/150‑1 (MQ=255)
       aCTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAAcc                                                                                                                                                >  2:3321011/1‑150 (MQ=255)
        cTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCa                                                                                                                                               <  2:4203897/150‑1 (MQ=255)
        cTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCa                                                                                                                                               <  2:3463166/150‑1 (MQ=255)
        cTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCa                                                                                                                                               >  1:3731995/1‑150 (MQ=255)
          ttATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAAc                                                                                                                                             <  1:4231153/150‑1 (MQ=255)
           tATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACg                                                                                                                                            >  2:2719228/1‑150 (MQ=255)
            aTCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACgg                                                                                                                                           >  1:813736/1‑150 (MQ=255)
               ggCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATg                                                                                                                                        <  1:2109763/150‑1 (MQ=255)
               ggCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATg                                                                                                                                        >  2:700286/1‑150 (MQ=255)
               ggCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCCACGGATg                                                                                                                                        <  2:3243073/150‑1 (MQ=255)
                 cAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATggg                                                                                                                                      >  1:1205040/1‑150 (MQ=255)
                  aaTCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGc                                                                                                                                     >  1:271973/1‑150 (MQ=255)
                  aaTCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGc                                                                                                                                     >  1:272358/1‑150 (MQ=255)
                   aTCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCa                                                                                                                                    >  2:3626918/1‑150 (MQ=255)
                   aTCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCAGCTTCGATAACCAACGGATGGGCa                                                                                                                                    <  2:2513819/150‑1 (MQ=255)
                    tCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAg                                                                                                                                   <  1:718510/150‑1 (MQ=255)
                    tCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAg                                                                                                                                   >  1:3723323/1‑150 (MQ=255)
                    tCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAg                                                                                                                                   >  1:4708216/1‑150 (MQ=255)
                      gTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGaa                                                                                                                                 <  1:2195369/150‑1 (MQ=255)
                           aaCAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGAc                                                                                                                            <  2:16355/150‑1 (MQ=255)
                           aaCAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGAc                                                                                                                            >  2:4636523/1‑150 (MQ=255)
                            aCAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACt                                                                                                                           <  2:3187584/150‑1 (MQ=255)
                               aaTTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGaa                                                                                                                        >  1:1019686/1‑150 (MQ=255)
                                aTTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAAc                                                                                                                       >  2:581877/1‑150 (MQ=255)
                                  tCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACat                                                                                                                     >  1:1249684/1‑150 (MQ=255)
                                  tCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACat                                                                                                                     >  1:1249987/1‑150 (MQ=255)
                                   cgcgAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACata                                                                                                                    >  1:861376/1‑150 (MQ=255)
                                      gAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGc                                                                                                                 >  2:4525262/1‑150 (MQ=255)
                                        aaGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCg                                                                                                               <  2:3499930/150‑1 (MQ=255)
                                        aaGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCg                                                                                                               <  2:2058351/150‑1 (MQ=255)
                                        aaGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCg                                                                                                               >  1:4647438/1‑150 (MQ=255)
                                            cTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGatat                                                                                                           <  1:4401012/150‑1 (MQ=255)
                                            cTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGatat                                                                                                           >  1:2800107/1‑150 (MQ=255)
                                              gCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGt                                                                                                         <  2:604945/150‑1 (MQ=255)
                                                 ggTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCa                                                                                                      >  1:2197620/1‑150 (MQ=255)
                                                 ggTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCa                                                                                                      >  1:357466/1‑150 (MQ=255)
                                                     tCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAg                                                                                                  >  1:704132/1‑150 (MQ=255)
                                                       aTTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCt                                                                                                >  2:4445514/1‑150 (MQ=255)
                                                        ttCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCtt                                                                                               <  1:3615533/150‑1 (MQ=255)
                                                         tCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCATCAGCttt                                                                                              >  2:2121469/1‑150 (MQ=255)
                                                         tCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCttt                                                                                              >  1:3118873/1‑150 (MQ=255)
                                                           cGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTAt                                                                                            <  1:3045684/150‑1 (MQ=255)
                                                           cGGGATCACCGCCAAATTGATTGACCTGGGTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTAt                                                                                            <  2:3856350/150‑1 (MQ=255)
                                                            ggTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCAATATGTTCACCAGCTTTATg                                                                                           >  1:2424274/1‑150 (MQ=255)
                                                               ataACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCg                                                                                        <  1:591723/150‑1 (MQ=255)
                                                                   ccGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCa                                                                                    <  1:2900423/150‑1 (MQ=255)
                                                                     gCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAAt                                                                                  <  2:2401383/150‑1 (MQ=255)
                                                                      ccAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAAtt                                                                                 <  2:1738720/150‑1 (MQ=255)
                                                                      ccAAAGTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAAtt                                                                                 <  2:1714549/150‑1 (MQ=255)
                                                                       cAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTa                                                                                <  2:1509770/150‑1 (MQ=255)
                                                                        aaaTTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGGTCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTaa                                                                               <  2:3471493/150‑1 (MQ=255)
                                                                           ttgattgaCCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGt                                                                            >  2:1772368/1‑150 (MQ=255)
                                                                            tgattgaCCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGtt                                                                           <  2:1213050/150‑1 (MQ=255)
                                                                             gattgaCCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGttt                                                                          <  1:2073022/150‑1 (MQ=255)
                                                                                tgaCCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCa                                                                       <  2:4572860/150‑1 (MQ=255)
                                                                                 gaCCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCAt                                                                      <  1:2076641/150‑1 (MQ=255)
                                                                                 gaCCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCAt                                                                      >  1:3281325/1‑150 (MQ=255)
                                                                                    cTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGtt                                                                   >  1:2544902/1‑150 (MQ=255)
                                                                                     tGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGttt                                                                  <  1:2367719/150‑1 (MQ=255)
                                                                                     tGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGttt                                                                  <  2:1228278/150‑1 (MQ=255)
                                                                                      ggTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTAACCAGCTTTATGCCGGGCAATTAACGTGTTCATGtttt                                                                 >  1:4075330/1‑150 (MQ=255)
                                                                                       gTTTGACGTTGCTTCAATCAGCACTGTGCGCGTGCTGTAGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTCTGCCGGGCAATTAACGTTTTCATGttttt                                                                <  2:435644/150‑1 (MQ=255)
                                                                                           gACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCt                                                            >  1:682684/1‑150 (MQ=255)
                                                                                               ttGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGtt                                                        >  2:3813368/1‑150 (MQ=255)
                                                                                               ttGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGtt                                                        <  1:1062345/150‑1 (MQ=255)
                                                                                                tGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTa                                                       <  2:694489/150‑1 (MQ=255)
                                                                                                 gCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTaa                                                      >  2:1899658/1‑150 (MQ=255)
                                                                                                   ttCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAAt                                                    <  2:3509760/150‑1 (MQ=255)
                                                                                                       aTCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCa                                                >  1:1705534/1‑150 (MQ=255)
                                                                                                        tCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAg                                               <  2:3375016/150‑1 (MQ=255)
                                                                                                        tCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAg                                               <  2:4404581/150‑1 (MQ=255)
                                                                                                         cAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAgg                                              <  1:3469249/150‑1 (MQ=255)
                                                                                                          aGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGt                                             <  1:769258/150‑1 (MQ=255)
                                                                                                                ttGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACATACTGAACATATGCCGATATGTTAACCAGCTTTATGCCGGGCAAGTAACGTTTTCATGTTTTTTCCTTGTTAAATGTCAGGTGCGTTa                                       >  2:2017877/1‑150 (MQ=255)
                                                                                                                       tGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTAtgccctg                                <  2:3709430/150‑1 (MQ=255)
                                                                                                                           gTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATgccctgccct                            <  1:2485497/150‑1 (MQ=255)
                                                                                                                            ttGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTc                           <  2:2900435/150‑1 (MQ=255)
                                                                                                                            ttGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTc                           <  1:658815/150‑1 (MQ=255)
                                                                                                                             tGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAAGTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTCTGCCCTGCCCTCg                          >  2:3122704/1‑150 (MQ=255)
                                                                                                                               cGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTAACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTACTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCa                        >  1:836714/1‑150 (MQ=255)
                                                                                                                                 aTCAAATGCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGc                      >  1:629487/1‑150 (MQ=255)
                                                                                                                                   cAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCca                    >  1:3909458/1‑150 (MQ=255)
                                                                                                                                   cAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCca                    >  2:1496853/1‑150 (MQ=255)
                                                                                                                                        gCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATc               <  2:1124101/150‑1 (MQ=255)
                                                                                                                                        gCCAGCGCCGCTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATc               >  2:4875591/1‑150 (MQ=255)
                                                                                                                                             cgccgcTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAAGTCACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCaa          >  2:2666192/1‑150 (MQ=255)
                                                                                                                                              gccgcTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAAt         >  2:1004813/1‑150 (MQ=255)
                                                                                                                                                cgcTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCa       <  1:1449481/150‑1 (MQ=255)
                                                                                                                                                 gcTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCAc      <  1:4112609/150‑1 (MQ=255)
                                                                                                                                                  cTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACt     >  1:4423308/1‑150 (MQ=255)
                                                                                                                                                  cTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACt     >  1:4388811/1‑150 (MQ=255)
                                                                                                                                                  cTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACt     >  2:2369027/1‑150 (MQ=255)
                                                                                                                                                  cTTAGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACt     >  1:2600242/1‑150 (MQ=255)
                                                                                                                                                    tCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACttt   <  2:4733281/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGGCAATTAACGTTTTCATGTTTTTTCCTTGTTAAATGGCAGGTGCGTTATGCCCTGCCCTCGCAGCCACAATCGGCAATCACTTTG  >  CP009273/2169570‑2169868

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: ‑