Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,332,486 G→T 13.0% R800L (CGA→CTA)  acnA → aconitate hydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,332,4860GT13.0% 849.9 / 69.7 245R800L (CGA→CTA) acnAaconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (95/118);  new base T (23/9);  total (118/127)
Fisher's exact test for biased strand distribution p-value = 4.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 247 total aligned reads displayed.

TTATGATGCTGCGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGA  >  CP009273/1332337‑1332635
                                                                                                                                                     |                                                                                                                                                     
  atgatgCTGCGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGaa                                                                                                                                                     >  2:3625325/1‑150 (MQ=255)
  atgatgCTGCGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGaa                                                                                                                                                     >  2:3625341/1‑150 (MQ=255)
   tgatgCTGCGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAAt                                                                                                                                                    >  1:1852894/1‑150 (MQ=255)
   tgatgCTGCGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAAt                                                                                                                                                    >  1:1614011/1‑150 (MQ=255)
       gctgcGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCAc                                                                                                                                                <  1:2699938/150‑1 (MQ=255)
         tgcGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGAGATTGCCGAATCGTTTGAACGAATTCACCg                                                                                                                                              >  1:1685137/1‑150 (MQ=255)
          gcGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGt                                                                                                                                             >  1:3048906/1‑150 (MQ=255)
               gcgcTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGaa                                                                                                                                        <  1:3678249/150‑1 (MQ=255)
                cgcTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAAt                                                                                                                                       >  1:4437759/1‑150 (MQ=255)
                 gcTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAAtt                                                                                                                                      >  1:3254508/1‑150 (MQ=255)
                 gcTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAAtt                                                                                                                                      >  2:1582500/1‑150 (MQ=255)
                   tataAGCTGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTa                                                                                                                                    <  1:3880562/150‑1 (MQ=255)
                   tataAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTa                                                                                                                                    <  2:384299/150‑1 (MQ=255)
                     taAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAAt                                                                                                                                  <  2:4337612/150‑1 (MQ=255)
                       aGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCTGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTg                                                                                                                                >  2:1154322/1‑150 (MQ=255)
                       aGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTg                                                                                                                                >  2:1425147/1‑150 (MQ=255)
                        gCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTgg                                                                                                                               <  2:1790237/150‑1 (MQ=255)
                           ggAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCAt                                                                                                                            <  2:729205/150‑1 (MQ=255)
                             aGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATgg                                                                                                                          <  1:651799/150‑1 (MQ=255)
                             aGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATgg                                                                                                                          <  2:805490/150‑1 (MQ=255)
                              gCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATggg                                                                                                                         >  2:4241425/1‑150 (MQ=255)
                               cAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGc                                                                                                                        >  2:919721/1‑150 (MQ=255)
                                aaaCGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCa                                                                                                                       >  1:2473511/1‑150 (MQ=255)
                                 aaCGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCAt                                                                                                                      <  2:2021679/150‑1 (MQ=255)
                                    gccgcTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCt                                                                                                                   <  1:2238508/150‑1 (MQ=255)
                                     ccgcTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTg                                                                                                                  <  2:570103/150‑1 (MQ=255)
                                       gcTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGcc                                                                                                                >  2:3416503/1‑150 (MQ=255)
                                             ggTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGa                                                                                                          <  2:2658385/150‑1 (MQ=255)
                                               tGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAAt                                                                                                        >  1:2466953/1‑150 (MQ=255)
                                               tGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAAt                                                                                                        >  1:2106380/1‑150 (MQ=255)
                                                  ttGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTc                                                                                                     >  2:631093/1‑150 (MQ=255)
                                                   tGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCATTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTcc                                                                                                    <  1:767339/150‑1 (MQ=255)
                                                    gCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCg                                                                                                   <  2:2342204/150‑1 (MQ=255)
                                                        ggAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAg                                                                                               <  1:2528756/150‑1 (MQ=255)
                                                        ggAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAg                                                                                               <  2:647594/150‑1 (MQ=255)
                                                         gacagagTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAgg                                                                                              <  2:1913140/147‑1 (MQ=255)
                                                         gAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGTACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAgg                                                                                              <  2:247029/150‑1 (MQ=255)
                                                           aaGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCg                                                                                            <  1:4566155/150‑1 (MQ=255)
                                                            agagTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGt                                                                                           <  1:2968685/150‑1 (MQ=255)
                                                              agTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTaa                                                                                         >  1:3818393/1‑150 (MQ=255)
                                                                tATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAAcg                                                                                       >  1:134461/1‑150 (MQ=255)
                                                                tATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAAcg                                                                                       >  2:2819947/1‑150 (MQ=255)
                                                                tATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAAcg                                                                                       <  2:3478119/150‑1 (MQ=255)
                                                                tATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAAcg                                                                                       <  2:608982/150‑1 (MQ=255)
                                                                tATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAAcg                                                                                       <  2:3478128/150‑1 (MQ=255)
                                                                tATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAAcg                                                                                       <  1:1060296/150‑1 (MQ=255)
                                                                   ggATCAGGCTACAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGt                                                                                    >  2:1853042/1‑150 (MQ=255)
                                                                    gATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTa                                                                                   >  2:1182771/1‑150 (MQ=255)
                                                                    gATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTa                                                                                   >  1:3572206/1‑150 (MQ=255)
                                                                    gATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTa                                                                                   <  1:2003709/150‑1 (MQ=255)
                                                                        aGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGATTTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAAc                                                                               <  1:2684585/150‑1 (MQ=255)
                                                                          gCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGt                                                                             >  1:1158170/1‑150 (MQ=255)
                                                                           cTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGtt                                                                            <  1:1833690/150‑1 (MQ=255)
                                                                             ccAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAg                                                                          <  1:2724083/150‑1 (MQ=255)
                                                                                  cGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTa                                                                     <  2:599336/150‑1 (MQ=255)
                                                                                     gACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAAcc                                                                  <  1:4489622/150‑1 (MQ=255)
                                                                                       cTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCgg                                                                >  1:2275559/1‑150 (MQ=255)
                                                                                           gcggcAAAATGTCCGCGTCTGCTTGGTATTCGTGTTGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaa                                                            <  1:2446576/150‑1 (MQ=255)
                                                                                           gcggcAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaa                                                            <  1:4022747/150‑1 (MQ=255)
                                                                                           gcggcAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaa                                                            >  1:1516729/1‑150 (MQ=255)
                                                                                            cggcAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaag                                                           <  2:2511297/150‑1 (MQ=255)
                                                                                              gcAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagag                                                         <  1:429596/150‑1 (MQ=255)
                                                                                              gcAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagag                                                         <  2:3008719/150‑1 (MQ=255)
                                                                                                aaaaGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagagaa                                                       <  1:1178987/150‑1 (MQ=255)
                                                                                                aaaaGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagagaa                                                       >  1:1920414/1‑150 (MQ=255)
                                                                                                 aaaGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagagaag                                                      >  2:1566103/1‑150 (MQ=255)
                                                                                                 aaaGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagagaag                                                      <  2:2639792/150‑1 (MQ=255)
                                                                                                  aaGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagagaaga                                                     >  1:343041/1‑150 (MQ=255)
                                                                                                  aaGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGgaagagaaga                                                     <  1:919721/150‑1 (MQ=255)
                                                                                                    ggTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGAtt                                                   <  1:3452870/150‑1 (MQ=255)
                                                                                                          cgTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATAtt                                             >  1:37162/1‑150 (MQ=255)
                                                                                                          cgTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATAtt                                             <  1:1136035/150‑1 (MQ=255)
                                                                                                           gTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTg                                            <  1:891003/150‑1 (MQ=255)
                                                                                                           gTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTg                                            <  2:4521408/150‑1 (MQ=255)
                                                                                                             cTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGc                                          >  2:4250563/1‑150 (MQ=255)
                                                                                                                cTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGAt                                       <  1:4586481/150‑1 (MQ=255)
                                                                                                                  tGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCt                                     >  2:2081016/1‑150 (MQ=255)
                                                                                                                  tGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCt                                     >  1:870899/1‑150 (MQ=255)
                                                                                                                  tGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCt                                     >  1:871007/1‑150 (MQ=255)
                                                                                                                   ggTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTg                                    >  2:2720320/1‑150 (MQ=255)
                                                                                                                    gTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGTGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACTTTATGGCTAACCGCGGAAGAGAAGATTGATATTGGCGATCTGc                                   >  1:2814786/1‑150 (MQ=255)
                                                                                                                          gtgtGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAAcc                             >  1:2541267/1‑150 (MQ=255)
                                                                                                                          gtgtGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAAcc                             >  2:1330430/1‑150 (MQ=255)
                                                                                                                           tgtgGTGATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCt                            >  2:4019949/1‑150 (MQ=255)
                                                                                                                           tgtgGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCt                            <  2:3364930/150‑1 (MQ=255)
                                                                                                                            gtggtgATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTa                           <  2:4300815/150‑1 (MQ=255)
                                                                                                                             tggtgATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTAc                          <  1:478194/150‑1 (MQ=255)
                                                                                                                              ggtgATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACa                         >  2:3119127/1‑150 (MQ=255)
                                                                                                                               gtgATTGCCGAATCGTTTGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACaa                        >  2:3920168/1‑150 (MQ=255)
                                                                                                                                tgATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAAc                       >  2:4575129/1‑150 (MQ=255)
                                                                                                                                   ttGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCg                    <  2:3238995/150‑1 (MQ=255)
                                                                                                                                   ttGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCg                    >  2:1039172/1‑150 (MQ=255)
                                                                                                                                       cGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGgcgc                >  2:1931903/1‑150 (MQ=255)
                                                                                                                                              tttGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGTCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGtt         <  2:4430989/150‑1 (MQ=255)
                                                                                                                                              tttGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGtt         <  2:3750061/150‑1 (MQ=255)
                                                                                                                                              tttGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGtt         <  2:3749939/150‑1 (MQ=255)
                                                                                                                                                tGAACTAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTcc       <  2:1106773/150‑1 (MQ=255)
                                                                                                                                                 gAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCg      >  2:450399/1‑150 (MQ=255)
                                                                                                                                                 gAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCg      <  1:4380436/150‑1 (MQ=255)
                                                                                                                                                     gAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGa  <  2:233667/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATGATGCTGCGATGCGCTATAAGCAGGAGCAAACGCCGCTGGCGGTGATTGCCGGGAAAGAGTATGGATCAGGCTCCAGTCGTGACTGGGCGGCAAAAGGTCCGCGTCTGCTTGGTATTCGTGTGGTGATTGCCGAATCGTTTGAACGAATTCACCGTTCGAATTTAATTGGCATGGGCATCCTGCCGCTGGAATTTCCGCAAGGCGTAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGA  >  CP009273/1332337‑1332635

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: