Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 166,662 C→A 89.9% S133* (TCG→TAG)  fhuC → iron‑hydroxamate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273166,6620CA89.9% 401.2 / 40.8 158S133* (TCG→TAG) fhuCiron‑hydroxamate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (5/11);  new base A (68/74);  total (73/85)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 158 total aligned reads displayed.

CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGAT  >  CP009273/166517‑166801
                                                                                                                                                 |                                                                                                                                           
cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                           >  1:2054578/1‑148 (MQ=255)
  cAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTg                                                                                                                                         <  1:1213875/148‑1 (MQ=255)
    aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGt                                                                                                                                       >  1:817445/1‑148 (MQ=255)
    aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGt                                                                                                                                       >  1:817487/1‑148 (MQ=255)
        ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                   <  1:1358700/148‑1 (MQ=255)
        ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                   <  2:1553492/148‑1 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCt                                                                                                                                  <  1:214641/148‑1 (MQ=255)
          cTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCtt                                                                                                                                 <  1:402004/148‑1 (MQ=255)
            tATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaa                                                                                                                               >  2:1830572/1‑148 (MQ=255)
             aTTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaa                                                                                                                              <  1:1986631/148‑1 (MQ=255)
             aTTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaa                                                                                                                              <  1:2783851/148‑1 (MQ=255)
               ttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAAc                                                                                                                            >  2:2287909/1‑148 (MQ=255)
                 gCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCg                                                                                                                          <  1:246441/148‑1 (MQ=255)
                   cGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCt                                                                                                                        <  1:1316304/148‑1 (MQ=255)
                        cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGgcgc                                                                                                                   >  2:878194/1‑148 (MQ=255)
                          gcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCAt                                                                                                                 <  2:787227/148‑1 (MQ=255)
                            ttCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCg                                                                                                               >  1:2440717/1‑148 (MQ=255)
                            ttCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCg                                                                                                               >  1:544083/1‑148 (MQ=255)
                             tcctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCgg                                                                                                              >  2:512865/1‑148 (MQ=255)
                             tcctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCgg                                                                                                              <  2:2195695/148‑1 (MQ=255)
                                tccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                           <  1:26396/148‑1 (MQ=255)
                                tccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                           >  2:3159263/1‑148 (MQ=255)
                                 ccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                          >  1:278261/1‑148 (MQ=255)
                                 ccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                          <  1:1441497/148‑1 (MQ=255)
                                     cAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGa                                                                                                      <  1:786442/148‑1 (MQ=255)
                                      aGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGAt                                                                                                     <  2:3017575/148‑1 (MQ=255)
                                      aGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGAt                                                                                                     <  2:129322/148‑1 (MQ=255)
                                        aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAg                                                                                                   >  2:539444/1‑148 (MQ=255)
                                        aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAg                                                                                                   >  1:1784998/1‑148 (MQ=255)
                                        aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAg                                                                                                   >  2:1735297/1‑148 (MQ=255)
                                        aaGGGAAGACCGTGCGGGAACTGGTGGCGAATGGTCGTTACCCGTGGCATGGCGAGCTGGGGCGCTTTGGGGAGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTATCAGGTTGGCATAAATCCGATGGCGCATCACCTCGTAGATAg                                                                                                   >  2:1655219/1‑148 (MQ=255)
                                             aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctct                                                                                              <  2:2457888/148‑1 (MQ=255)
                                             aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctct                                                                                              <  2:2457960/148‑1 (MQ=255)
                                             aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctct                                                                                              >  1:289325/1‑148 (MQ=255)
                                              tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                             <  1:1175533/148‑1 (MQ=255)
                                              tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                             <  1:1175515/148‑1 (MQ=255)
                                              tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                             >  2:2666904/1‑148 (MQ=255)
                                              tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                             >  2:2666910/1‑148 (MQ=255)
                                               gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTctctct                                                                                            >  1:3042435/1‑148 (MQ=255)
                                                         gAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAAc                                                                                  >  1:1670838/1‑148 (MQ=255)
                                                                 ggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCgg                                                                          <  2:1834279/148‑1 (MQ=255)
                                                                     aTTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGt                                                                      <  1:1452020/148‑1 (MQ=255)
                                                                       tGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                    >  2:284569/1‑148 (MQ=255)
                                                                       tGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                    <  1:2285252/148‑1 (MQ=255)
                                                                       tGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                    <  1:2285283/148‑1 (MQ=255)
                                                                        ggTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGa                                                                   <  2:923332/148‑1 (MQ=255)
                                                                        ggTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGa                                                                   <  2:1528222/148‑1 (MQ=255)
                                                                        ggTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGa                                                                   <  2:2713622/148‑1 (MQ=255)
                                                                             ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGcc                                                              <  1:1736352/148‑1 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGACGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                             >  1:1379406/1‑148 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGACCGTCAGCGGGCGTGGATCGCCa                                                             >  1:1379437/1‑148 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                             >  1:305829/1‑148 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                             >  1:306252/1‑148 (MQ=255)
                                                                               aCCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                            >  1:1572111/1‑148 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                           >  2:2230672/1‑148 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                           >  2:2230653/1‑148 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                           >  2:2227697/1‑148 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                           >  2:1590430/1‑148 (MQ=255)
                                                                                 ccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGc                                                          >  1:3182518/1‑148 (MQ=255)
                                                                                   gTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtg                                                        <  1:1073570/148‑1 (MQ=255)
                                                                                    tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtgg                                                       <  1:296173/148‑1 (MQ=255)
                                                                                     ggCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggt                                                      <  2:2617580/148‑1 (MQ=255)
                                                                                       cATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                                    >  1:3109864/1‑148 (MQ=255)
                                                                                       cATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                                    >  1:3109499/1‑148 (MQ=255)
                                                                                         tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                                  <  2:1662141/148‑1 (MQ=255)
                                                                                               gcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGAt                                            <  2:39652/148‑1 (MQ=255)
                                                                                               gcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGACCGCCATGCTGGTGGCGCAGGAt                                            <  1:3159263/148‑1 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                          <  2:817487/148‑1 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                          <  1:3300457/148‑1 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                          <  2:817445/148‑1 (MQ=255)
                                                                                                   gggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGcc                                        >  1:1131308/1‑148 (MQ=255)
                                                                                                    ggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCg                                       >  2:1602223/1‑148 (MQ=255)
                                                                                                     gcgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGt                                      <  2:2922927/148‑1 (MQ=255)
                                                                                                     gcgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGt                                      >  1:7209/1‑148 (MQ=255)
                                                                                                       gcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTg                                    >  2:1167321/1‑148 (MQ=255)
                                                                                                        cTTTGGGGCGGCAGATCGCGAAACAGTCGAGGACGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGt                                   <  2:2964303/148‑1 (MQ=255)
                                                                                                          ttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCt                                 <  1:928913/148‑1 (MQ=255)
                                                                                                             gggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGtt                              >  1:1234217/1‑148 (MQ=255)
                                                                                                              ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                             <  2:71709/148‑1 (MQ=255)
                                                                                                              ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                             <  2:693267/148‑1 (MQ=255)
                                                                                                              ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                             >  1:119632/1‑148 (MQ=255)
                                                                                                              ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                             <  1:1923118/148‑1 (MQ=255)
                                                                                                                cggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                           <  2:1449764/148‑1 (MQ=255)
                                                                                                                  gcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcg                         <  1:3026082/148‑1 (MQ=255)
                                                                                                                   cAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                        >  1:273767/1‑148 (MQ=255)
                                                                                                                   cAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                        >  1:274396/1‑148 (MQ=255)
                                                                                                                     gATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgacg                      >  1:458323/1‑148 (MQ=255)
                                                                                                                      aTCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgacga                     <  1:3264795/148‑1 (MQ=255)
                                                                                                                      aTCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgacga                     >  1:2033934/1‑148 (MQ=255)
                                                                                                                        cgcgAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAAc                   >  2:1176144/1‑148 (MQ=255)
                                                                                                                           gAAAAAGTCGAGGCAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGa                <  1:180828/148‑1 (MQ=255)
                                                                                                                           gAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGa                <  2:3041424/148‑1 (MQ=255)
                                                                                                                            aaaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAc               >  2:1009870/1‑148 (MQ=255)
                                                                                                                            aaaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAc               >  2:1187094/1‑148 (MQ=255)
                                                                                                                             aaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAcc              <  1:2351477/148‑1 (MQ=255)
                                                                                                                              aaaGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCt             <  1:2513210/148‑1 (MQ=255)
                                                                                                                                 gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGTGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCgg          <  2:2871136/148‑1 (MQ=255)
                                                                                                                                     aGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCt      <  1:941637/148‑1 (MQ=255)
                                                                                                                                       gAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTgg    <  2:1254150/148‑1 (MQ=255)
                                                                                                                                         aGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGat  >  2:2591545/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGAT  >  CP009273/166517‑166801

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: