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,6620CA97.3% 1024.4 / 15.2 330S133* (TCG→TAG) fhuCiron‑hydroxamate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base A (140/181);  total (145/185)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 331 total aligned reads displayed.

CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGCCC  >  CP009273/166517‑166808
                                                                                                                                                 |                                                                                                                                                  
cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTg                                                                                                                                                <  1:2199546/150‑1 (MQ=255)
    aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTg                                                                                                                                            <  1:1012540/150‑1 (MQ=255)
     aGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTgg                                                                                                                                           <  2:2429007/150‑1 (MQ=255)
       tGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCt                                                                                                                                         <  2:5141014/150‑1 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTa                                                                                                                                       <  2:3345908/150‑1 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTa                                                                                                                                       >  2:3093243/1‑150 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTa                                                                                                                                       <  1:5107593/150‑1 (MQ=255)
           ttATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaa                                                                                                                                     >  1:4849933/1‑150 (MQ=255)
            tATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaaa                                                                                                                                    <  2:1218342/150‑1 (MQ=255)
              tttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTGAAAAcc                                                                                                                                  >  1:4780036/1‑150 (MQ=255)
                tGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGc                                                                                                                                >  1:603933/1‑150 (MQ=255)
                  ccGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTg                                                                                                                              <  1:3513766/150‑1 (MQ=255)
                   cGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGGGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTgg                                                                                                                             <  2:4856608/150‑1 (MQ=255)
                   cGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTgg                                                                                                                             >  1:1166760/1‑150 (MQ=255)
                        cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCAt                                                                                                                        <  2:4677786/150‑1 (MQ=255)
                         agcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                       >  2:1381479/1‑150 (MQ=255)
                           cTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCgg                                                                                                                     >  1:1998350/1‑150 (MQ=255)
                              cctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                                  >  2:1914915/1‑150 (MQ=255)
                              cctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                                  >  2:697506/1‑150 (MQ=255)
                               ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                                 >  1:1069297/1‑150 (MQ=255)
                               ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                                 <  1:3842228/150‑1 (MQ=255)
                               ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                                 <  1:3252055/150‑1 (MQ=255)
                                  cGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCg                                                                                                              <  1:4382136/150‑1 (MQ=255)
                                  cGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCg                                                                                                              >  1:252531/1‑150 (MQ=255)
                                    gCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGAt                                                                                                            <  2:4152292/150‑1 (MQ=255)
                                     cAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                                           <  2:4758836/150‑1 (MQ=255)
                                     cAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTAACCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCTGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                                           <  2:1719165/150‑1 (MQ=255)
                                      aGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAg                                                                                                          >  2:5006559/1‑150 (MQ=255)
                                      aGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCGTTGGGGCGGCAGATCGCGAACAAGTCGAGGAAGCTAGCTAGCTGGGTGGCTTCTAACCGCTGGCGCATCGGCGGGGCGagag                                                                                                          >  2:3031338/1‑150 (MQ=255)
                                       gAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                                         <  1:1668635/150‑1 (MQ=255)
                                        aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtc                                                                                                        >  2:3127242/1‑150 (MQ=255)
                                         aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACACGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtct                                                                                                       <  1:3978330/150‑1 (MQ=255)
                                             aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAATGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTctctct                                                                                                   >  2:4101594/1‑150 (MQ=255)
                                             aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTctctct                                                                                                   >  2:1256274/1‑150 (MQ=255)
                                              tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTg                                                                                                  <  1:4593778/150‑1 (MQ=255)
                                                 ccGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTggcg                                                                                               >  1:4449881/1‑150 (MQ=255)
                                                  cgtgcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTggcgg                                                                                              >  2:2539508/1‑150 (MQ=255)
                                                     gcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGa                                                                                           <  2:2998037/150‑1 (MQ=255)
                                                      cgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGaa                                                                                          >  2:3934346/1‑150 (MQ=255)
                                                       gtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAAc                                                                                         >  1:2859893/1‑150 (MQ=255)
                                                        tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACg                                                                                        <  2:3630129/150‑1 (MQ=255)
                                                          aaCTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                                      <  1:4919463/150‑1 (MQ=255)
                                                             tggtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGc                                                                                   <  2:2800714/150‑1 (MQ=255)
                                                             tggtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGc                                                                                   <  2:1299877/150‑1 (MQ=255)
                                                             tggtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGc                                                                                   <  1:1178805/150‑1 (MQ=255)
                                                               gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCgg                                                                                 >  1:2958754/1‑150 (MQ=255)
                                                                tggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTCTCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCggg                                                                                >  2:2565870/1‑150 (MQ=255)
                                                                 ggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGc                                                                               <  1:3724786/150‑1 (MQ=255)
                                                                  gCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCg                                                                              <  2:3466125/150‑1 (MQ=255)
                                                                    gATTGGTCGATACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTg                                                                            <  1:4199753/150‑1 (MQ=255)
                                                                     aTTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTgg                                                                           >  2:2229418/1‑150 (MQ=255)
                                                                             ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                   >  1:1390/1‑150 (MQ=255)
                                                                               aCCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGc                                                                 <  2:3811054/150‑1 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCt                                                                <  1:531161/150‑1 (MQ=255)
                                                                                cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCt                                                                <  1:3786489/150‑1 (MQ=255)
                                                                                     ggCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                                           <  1:1950322/150‑1 (MQ=255)
                                                                                      gCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgc                                                          <  2:2362087/150‑1 (MQ=255)
                                                                                       cATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                                         <  1:3641782/150‑1 (MQ=255)
                                                                                        aTGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcgc                                                        >  2:4285949/1‑150 (MQ=255)
                                                                                        aTGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcgc                                                        >  2:2030207/1‑150 (MQ=255)
                                                                                         tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCa                                                       >  1:2529401/1‑150 (MQ=255)
                                                                                         tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCa                                                       <  2:4034562/150‑1 (MQ=255)
                                                                                          ggCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAg                                                      >  2:4122065/1‑150 (MQ=255)
                                                                                             gcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGAt                                                   <  1:3118851/150‑1 (MQ=255)
                                                                                               gcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGACGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAg                                                 <  1:935986/150‑1 (MQ=255)
                                                                                                cTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGc                                                >  1:3662705/1‑150 (MQ=255)
                                                                                                 tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGcc                                               >  1:1430846/1‑150 (MQ=255)
                                                                                                    ggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGtt                                            <  2:3423078/150‑1 (MQ=255)
                                                                                                     gcgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTg                                           <  1:1564212/150‑1 (MQ=255)
                                                                                                      cgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGt                                          <  2:4018080/150‑1 (MQ=255)
                                                                                                         tttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTg                                       <  2:268598/150‑1 (MQ=255)
                                                                                                         tttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTg                                       >  2:4920584/1‑150 (MQ=255)
                                                                                                            ggggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                                    <  1:2320718/150‑1 (MQ=255)
                                                                                                             gggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGc                                   <  1:1914915/150‑1 (MQ=255)
                                                                                                              ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                  <  1:1303424/150‑1 (MQ=255)
                                                                                                              ggcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                  <  2:3733291/150‑1 (MQ=255)
                                                                                                               gcggcAGATCGCGAAAAAGTCGATGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTc                                 <  1:2467052/150‑1 (MQ=255)
                                                                                                                cggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcg                                <  1:3060924/150‑1 (MQ=255)
                                                                                                                cggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcg                                <  1:3356544/150‑1 (MQ=255)
                                                                                                                 ggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                               >  1:2682793/1‑150 (MQ=255)
                                                                                                                   cAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgacg                             >  2:1467685/1‑150 (MQ=255)
                                                                                                                       tCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAAcc                         <  2:283962/150‑1 (MQ=255)
                                                                                                                        cgcgAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCg                        <  1:4871459/150‑1 (MQ=255)
                                                                                                                        cgcgAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCg                        >  1:2552668/1‑150 (MQ=255)
                                                                                                                             aaaaGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTc                   <  1:2047635/150‑1 (MQ=255)
                                                                                                                                aGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgc                <  1:1493805/150‑1 (MQ=255)
                                                                                                                                 gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcg               <  2:863718/150‑1 (MQ=255)
                                                                                                                                 gTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCCCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcg               <  2:612853/150‑1 (MQ=255)
                                                                                                                                    gAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTg            >  2:1541812/1‑150 (MQ=255)
                                                                                                                                    gAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTg            >  1:3175179/1‑150 (MQ=255)
                                                                                                                                      ggAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGa          >  2:2212929/1‑150 (MQ=255)
                                                                                                                                        aaGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata        >  1:4832554/1‑150 (MQ=255)
                                                                                                                                        aaGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata        >  1:4183017/1‑150 (MQ=255)
                                                                                                                                        aaGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata        <  2:4573201/150‑1 (MQ=255)
                                                                                                                                        aaGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata        <  2:4544900/150‑1 (MQ=255)
                                                                                                                                          gCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATc      >  2:4802962/1‑150 (MQ=255)
                                                                                                                                          gCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATc      <  2:1166594/150‑1 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCg     >  1:4523366/1‑150 (MQ=255)
                                                                                                                                            tATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGc    >  1:1853833/1‑150 (MQ=255)
                                                                                                                                            tATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGc    >  1:2733802/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGCCC  >  CP009273/166517‑166808

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: