Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 166,662 C→A 93.8% 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,6620CA93.8% 512.1 / 21.7 181S133* (TCG→TAG) fhuCiron‑hydroxamate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (3/9);  new base A (86/83);  total (89/92)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

CCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATAT  >  CP009273/166516‑166803
                                                                                                                                                  |                                                                                                                                             
 cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGc                                                                                                                                             <  2:522764/148‑1 (MQ=255)
 cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                             >  1:3213081/1‑148 (MQ=255)
 cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                             >  2:3515733/1‑148 (MQ=255)
 cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGc                                                                                                                                             <  1:864793/148‑1 (MQ=255)
 cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATATCGc                                                                                                                                             >  1:1219015/1‑148 (MQ=255)
   cAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTg                                                                                                                                           <  1:1515826/148‑1 (MQ=255)
     aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTATCTGGt                                                                                                                                         >  2:337268/1‑148 (MQ=255)
     aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGt                                                                                                                                         <  2:1357707/148‑1 (MQ=255)
     aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGt                                                                                                                                         >  2:3062955/1‑148 (MQ=255)
      aGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGtt                                                                                                                                        >  1:2482456/1‑148 (MQ=255)
       gTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTg                                                                                                                                       <  2:2909366/148‑1 (MQ=255)
         ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                     <  2:626671/148‑1 (MQ=255)
         ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                     <  1:3163788/148‑1 (MQ=255)
         ggCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                     <  2:2690520/148‑1 (MQ=255)
          gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCt                                                                                                                                    <  2:3043826/148‑1 (MQ=255)
          gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCt                                                                                                                                    >  2:1705736/1‑148 (MQ=255)
            ttATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGATCTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTTGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGTTTGGCTTa                                                                                                                                  <  2:749509/148‑1 (MQ=255)
             tATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaa                                                                                                                                 >  2:2066390/1‑148 (MQ=255)
             tATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCGTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaa                                                                                                                                 >  2:2066387/1‑148 (MQ=255)
                ttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAAc                                                                                                                              <  2:1128782/148‑1 (MQ=255)
                 tGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAAcc                                                                                                                             >  1:3265593/1‑148 (MQ=255)
                  gCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCg                                                                                                                            <  1:585667/148‑1 (MQ=255)
                   ccGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTGACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGc                                                                                                                           <  2:3220138/148‑1 (MQ=255)
                      cagcagCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCATTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTgg                                                                                                                        <  2:545815/148‑1 (MQ=255)
                         cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGgcgc                                                                                                                     >  2:2628136/1‑148 (MQ=255)
                         cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGgcgc                                                                                                                     >  1:1141394/1‑148 (MQ=255)
                         cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGgcgc                                                                                                                     >  2:2995817/1‑148 (MQ=255)
                               cctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGc                                                                                                               >  1:3485922/1‑148 (MQ=255)
                                ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCt                                                                                                              <  1:3426505/148‑1 (MQ=255)
                                ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCt                                                                                                              <  1:3344193/148‑1 (MQ=255)
                                  ccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                            <  2:1242316/148‑1 (MQ=255)
                                  ccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                            <  1:163238/148‑1 (MQ=255)
                                          aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                                    >  2:1519101/1‑148 (MQ=255)
                                          aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                                    <  1:232190/148‑1 (MQ=255)
                                          aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                                    <  1:2027216/148‑1 (MQ=255)
                                            ggATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtct                                                                                                  <  2:2857291/148‑1 (MQ=255)
                                             gATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctc                                                                                                 >  1:1696656/1‑148 (MQ=255)
                                              aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctct                                                                                                >  1:2058743/1‑148 (MQ=255)
                                               tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                               <  2:1102284/148‑1 (MQ=255)
                                               tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                               <  2:167484/148‑1 (MQ=255)
                                               tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                               <  2:2017483/148‑1 (MQ=255)
                                               tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                               <  2:167518/148‑1 (MQ=255)
                                               tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                               <  1:3526386/148‑1 (MQ=255)
                                                gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTctctct                                                                                              <  1:842936/148‑1 (MQ=255)
                                                 aCCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTg                                                                                             >  2:2420720/1‑148 (MQ=255)
                                                      gcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTggcggc                                                                                        >  1:1436014/1‑148 (MQ=255)
                                                      gcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTggcggc                                                                                        >  1:1436186/1‑148 (MQ=255)
                                                         tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGaa                                                                                     <  2:1063149/148‑1 (MQ=255)
                                                         tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTAGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTATCTCTGGCGGCGaa                                                                                     >  1:2492295/1‑148 (MQ=255)
                                                             cTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                                 >  1:565565/1‑148 (MQ=255)
                                                             cTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                                 >  1:2351072/1‑148 (MQ=255)
                                                                 tggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCg                                                                             >  2:1612807/1‑148 (MQ=255)
                                                                 tggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCg                                                                             >  1:610440/1‑148 (MQ=255)
                                                                 tggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCg                                                                             <  2:3020652/148‑1 (MQ=255)
                                                                      aTTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGACCGTCAGCGGGCGt                                                                        >  1:445254/1‑148 (MQ=255)
                                                                      aTTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGt                                                                        <  2:1272884/148‑1 (MQ=255)
                                                                      aTTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGt                                                                        <  2:3401403/148‑1 (MQ=255)
                                                                      aTTGGTCGTGACCCGTGGCATGGCGCGCGGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGt                                                                        <  2:1272893/148‑1 (MQ=255)
                                                                        tGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                      <  2:2680115/148‑1 (MQ=255)
                                                                          gTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGAt                                                                    <  1:637842/148‑1 (MQ=255)
                                                                              ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGcc                                                                <  2:1831164/148‑1 (MQ=255)
                                                                              ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGcc                                                                >  2:481358/1‑148 (MQ=255)
                                                                               tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                               >  1:1181709/1‑148 (MQ=255)
                                                                               tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                               <  1:1948786/148‑1 (MQ=255)
                                                                                   cGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCt                                                           >  1:3328852/1‑148 (MQ=255)
                                                                                    gTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtg                                                          <  1:2200668/148‑1 (MQ=255)
                                                                                     tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtgg                                                         <  1:1685661/148‑1 (MQ=255)
                                                                                      ggCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggt                                                        <  1:3496048/148‑1 (MQ=255)
                                                                                          tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                                    >  1:2266902/1‑148 (MQ=255)
                                                                                          tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                                    >  1:3001777/1‑148 (MQ=255)
                                                                                          tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                                    >  1:478849/1‑148 (MQ=255)
                                                                                            gcgcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCa                                                  <  1:1256147/148‑1 (MQ=255)
                                                                                                    gggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGcc                                          <  2:268597/148‑1 (MQ=255)
                                                                                                      gcgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGt                                        <  1:3249791/148‑1 (MQ=255)
                                                                                                         cTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGt                                     <  1:947203/148‑1 (MQ=255)
                                                                                                           ttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCt                                   >  1:1323161/1‑148 (MQ=255)
                                                                                                           ttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATAGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCt                                   >  2:3170452/1‑148 (MQ=255)
                                                                                                            tGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTg                                  <  2:2172420/148‑1 (MQ=255)
                                                                                                               ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                               <  1:51872/148‑1 (MQ=255)
                                                                                                               ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                               <  1:2946641/148‑1 (MQ=255)
                                                                                                                gcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGc                              <  2:2845215/148‑1 (MQ=255)
                                                                                                                    cAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                          <  1:2401256/148‑1 (MQ=255)
                                                                                                                      gATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgacg                        >  1:1602570/1‑148 (MQ=255)
                                                                                                                        tCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGaa                      >  2:3049266/1‑148 (MQ=255)
                                                                                                                        tCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGaa                      >  2:3322225/1‑148 (MQ=255)
                                                                                                                             aaaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAc                 >  2:437210/1‑148 (MQ=255)
                                                                                                                              aaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAcc                >  1:2747903/1‑148 (MQ=255)
                                                                                                                               aaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCt               <  1:2793151/148‑1 (MQ=255)
                                                                                                                                 aGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCg             >  1:1830213/1‑148 (MQ=255)
                                                                                                                                  gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCgg            >  2:1144888/1‑148 (MQ=255)
                                                                                                                                  gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCgg            >  2:2897152/1‑148 (MQ=255)
                                                                                                                                  gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCgg            >  2:2897136/1‑148 (MQ=255)
                                                                                                                                  gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCgg            >  2:1122481/1‑148 (MQ=255)
                                                                                                                                  gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGGTGTCTGTTGCTCGACGAACCGACCTCgg            >  2:1122421/1‑148 (MQ=255)
                                                                                                                                     gAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcgc         <  1:3062955/148‑1 (MQ=255)
                                                                                                                                       ggAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTg       <  2:868401/148‑1 (MQ=255)
                                                                                                                                        gAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTgg      >  1:264428/1‑148 (MQ=255)
                                                                                                                                        gAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTgg      >  2:84840/1‑148 (MQ=255)
                                                                                                                                           gCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata   >  1:568380/1‑148 (MQ=255)
                                                                                                                                            cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGatat  <  2:1575420/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATAT  >  CP009273/166516‑166803

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: