Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 166,662 C→A 100% 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,6620CA98.2% 702.2 / 3.1 222S133* (TCG→TAG) fhuCiron‑hydroxamate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (99/119);  total (102/120)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 223 total aligned reads displayed.

CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGC  >  CP009273/166517‑166806
                                                                                                                                                 |                                                                                                                                                
cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGACGCTATCTAGc                                                                                                                                                >  1:3075535/1‑148 (MQ=255)
cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                                <  2:1191713/148‑1 (MQ=255)
cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                                >  2:1609486/1‑148 (MQ=255)
cGCAAAGTGGCTTATTTGCCGCAGCAGCCTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                                >  2:2656602/1‑148 (MQ=255)
 gCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCt                                                                                                                                               >  2:2187382/1‑148 (MQ=255)
  cAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTg                                                                                                                                              <  2:1411350/148‑1 (MQ=255)
  cAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTg                                                                                                                                              >  1:2251085/1‑148 (MQ=255)
   aaaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTgg                                                                                                                                             <  1:2707489/148‑1 (MQ=255)
       tGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTgg                                                                                                                                         <  2:514250/148‑1 (MQ=255)
        ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                        <  2:1614115/148‑1 (MQ=255)
        ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                        <  2:2624577/148‑1 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCt                                                                                                                                       <  1:973139/148‑1 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCt                                                                                                                                       >  1:2608791/1‑148 (MQ=255)
            tATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaa                                                                                                                                    <  2:633705/148‑1 (MQ=255)
              tttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaaa                                                                                                                                  <  1:2014510/148‑1 (MQ=255)
              tttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaaa                                                                                                                                  >  2:214792/1‑148 (MQ=255)
              tttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCGGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaaa                                                                                                                                  <  2:2329489/148‑1 (MQ=255)
                   cGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCt                                                                                                                             <  2:1382914/148‑1 (MQ=255)
                        cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGgcgc                                                                                                                        >  1:1084380/1‑148 (MQ=255)
                           cTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                     >  2:1323563/1‑148 (MQ=255)
                           cTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                     >  1:2886581/1‑148 (MQ=255)
                           cTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                     >  1:2962005/1‑148 (MQ=255)
                           cTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                     >  1:506053/1‑148 (MQ=255)
                            ttCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCg                                                                                                                    >  1:3070725/1‑148 (MQ=255)
                            ttCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCg                                                                                                                    >  1:1893583/1‑148 (MQ=255)
                                  cGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGt                                                                                                              >  2:2622850/1‑148 (MQ=255)
                                     cAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGa                                                                                                           <  2:437957/148‑1 (MQ=255)
                                       gAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                                         >  1:441188/1‑148 (MQ=255)
                                         aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                                       >  1:761012/1‑148 (MQ=255)
                                         aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                                       >  1:2695886/1‑148 (MQ=255)
                                           ggATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtct                                                                                                     >  2:3193047/1‑148 (MQ=255)
                                           ggATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtct                                                                                                     >  1:2922782/1‑148 (MQ=255)
                                                  cgtgcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTggc                                                                                              >  2:3055587/1‑148 (MQ=255)
                                                      cgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGgcggcg                                                                                          <  2:652423/148‑1 (MQ=255)
                                                        tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctc                                                                                                    <  1:3281371/136‑1 (MQ=255)
                                                         gAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAAc                                                                                       >  1:3159927/1‑148 (MQ=255)
                                                         gAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAAc                                                                                       >  2:743157/1‑148 (MQ=255)
                                                           aCTGGTGGCGATTGGTCGTTACCCGGGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAACAGGCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGt                                                                                     >  2:3299864/1‑148 (MQ=255)
                                                            cTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                                    <  2:1422177/148‑1 (MQ=255)
                                                              ggtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGATAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGAGGCGAACGTCAg                                                                                  <  1:2235350/148‑1 (MQ=255)
                                                               gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGc                                                                                 >  2:2003119/1‑148 (MQ=255)
                                                                tggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCg                                                                                >  2:2389817/1‑148 (MQ=255)
                                                                 ggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCgg                                                                               <  1:1594692/148‑1 (MQ=255)
                                                                  gCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCggg                                                                              >  1:73139/1‑148 (MQ=255)
                                                                  gCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCggg                                                                              >  1:1847821/1‑148 (MQ=255)
                                                                       tGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                         >  1:1697810/1‑148 (MQ=255)
                                                                       tGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                         <  2:3237681/148‑1 (MQ=255)
                                                                        ggTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGa                                                                        <  1:148212/148‑1 (MQ=255)
                                                                          tCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATc                                                                      >  2:2350083/1‑148 (MQ=255)
                                                                           cGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCg                                                                     >  2:870080/1‑148 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                                  <  2:1197889/148‑1 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                                  >  2:168070/1‑148 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCGGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                                  <  2:1781231/148‑1 (MQ=255)
                                                                               aCCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                 <  1:2851320/148‑1 (MQ=255)
                                                                               aCCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                 >  2:1717529/1‑148 (MQ=255)
                                                                               aCCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                 >  2:643127/1‑148 (MQ=255)
                                                                               aCCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                 >  2:1922710/1‑148 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                                <  2:2362128/148‑1 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                                <  1:3475291/148‑1 (MQ=255)
                                                                                cccGTGGCATGGCGCGCGGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATg                                                                <  2:599853/148‑1 (MQ=255)
                                                                                 ccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGc                                                               >  2:206303/1‑148 (MQ=255)
                                                                                   gTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtg                                                             <  1:1795834/148‑1 (MQ=255)
                                                                                       cATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                                         <  2:161875/148‑1 (MQ=255)
                                                                                             gcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAgg                                                   <  1:335943/148‑1 (MQ=255)
                                                                                             gcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAgg                                                   <  2:1634360/148‑1 (MQ=255)
                                                                                             gcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAgg                                                   <  1:2104629/148‑1 (MQ=255)
                                                                                              cgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGa                                                  <  1:2638934/148‑1 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                               <  2:493669/148‑1 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                               <  2:493978/148‑1 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                               <  2:2757257/148‑1 (MQ=255)
                                                                                                  ggggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGc                                              <  2:2705830/148‑1 (MQ=255)
                                                                                                     gcgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGt                                           <  2:2023236/148‑1 (MQ=255)
                                                                                                       gcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTg                                         <  2:2237014/148‑1 (MQ=255)
                                                                                                         tttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCGCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTc                                       <  1:2711152/148‑1 (MQ=255)
                                                                                                          ttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCt                                      >  2:3252569/1‑148 (MQ=255)
                                                                                                           tGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTg                                     <  1:2517869/148‑1 (MQ=255)
                                                                                                            ggggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGt                                    <  1:467154/148‑1 (MQ=255)
                                                                                                            ggggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGt                                    <  1:2906359/148‑1 (MQ=255)
                                                                                                             gggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGtt                                   >  1:3311622/1‑148 (MQ=255)
                                                                                                                cggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                <  1:1290673/148‑1 (MQ=255)
                                                                                                                cggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                >  1:2786253/1‑148 (MQ=255)
                                                                                                                cggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                <  1:3186917/148‑1 (MQ=255)
                                                                                                                  gcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcg                              <  2:2889280/148‑1 (MQ=255)
                                                                                                                   cAGATCGCGAAAAAGTCGAGGACGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                             >  2:125936/1‑148 (MQ=255)
                                                                                                                   cAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTTGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                             <  2:1017289/148‑1 (MQ=255)
                                                                                                                     gATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGGTGTCTGTTGCTcgacg                           >  2:932255/1‑148 (MQ=255)
                                                                                                                      aTCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgacga                          >  2:1946055/1‑148 (MQ=255)
                                                                                                                       tCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGaa                         <  2:1566167/148‑1 (MQ=255)
                                                                                                                           gAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGa                     <  1:2752952/148‑1 (MQ=255)
                                                                                                                             aaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAcc                   <  2:550243/148‑1 (MQ=255)
                                                                                                                                   cGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcg             <  2:1881917/148‑1 (MQ=255)
                                                                                                                                    gAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcgc            <  1:2187382/148‑1 (MQ=255)
                                                                                                                                      ggAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTg          >  1:3021649/1‑148 (MQ=255)
                                                                                                                                      ggAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTg          >  2:2586164/1‑148 (MQ=255)
                                                                                                                                      ggAAGCTATCGAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTg          <  2:2365494/148‑1 (MQ=255)
                                                                                                                                          gCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata      >  2:932115/1‑148 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGatat     <  2:1448171/148‑1 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGatat     <  2:2489887/148‑1 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGatat     >  1:1487697/1‑148 (MQ=255)
                                                                                                                                             aTCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCg   >  1:550615/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGC  >  CP009273/166517‑166806

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