Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,500,728 C→A 23.8% E183* (GAA→TAA)  fecE ← iron‑dicitrate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,500,7280CA23.8% 850.8 / 234.9 344E183* (GAA→TAA) fecEiron‑dicitrate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (155/107);  new base A (53/29);  total (208/136)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 353 total aligned reads displayed.

TCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAA  >  CP009273/4500579‑4500877
                                                                                                                                                     |                                                                                                                                                     
tctcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTc                                                                                                                                                       >  1:2734558/1‑150 (MQ=255)
tctcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTc                                                                                                                                                       <  1:4483511/150‑1 (MQ=255)
tctcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTc                                                                                                                                                       <  1:128793/150‑1 (MQ=255)
tctcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTc                                                                                                                                                       <  1:4671139/150‑1 (MQ=255)
 ctcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCg                                                                                                                                                      <  2:3323786/150‑1 (MQ=255)
 ctcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCg                                                                                                                                                      >  2:4067343/1‑150 (MQ=255)
  tcAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGc                                                                                                                                                     >  2:3569780/1‑150 (MQ=255)
   cAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGcc                                                                                                                                                    <  1:771561/150‑1 (MQ=255)
      cAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCAt                                                                                                                                                 >  1:2434828/1‑150 (MQ=255)
       aaTCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATc                                                                                                                                                >  1:1620641/1‑150 (MQ=255)
          ccTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAAc                                                                                                                                             >  1:3076234/1‑150 (MQ=255)
            tGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCg                                                                                                                                           >  1:4343222/1‑150 (MQ=255)
               ggTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCAt                                                                                                                                        <  2:4102251/150‑1 (MQ=255)
                 tcatcaCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCa                                                                                                                                      >  1:4086451/1‑150 (MQ=255)
                  catcaCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAg                                                                                                                                     >  2:3133205/1‑150 (MQ=255)
                  catcaCCTCTTCTGGTGTGACTTGCGCCATAACATGTCCGTTTGCAATTACCACCAGTTGATCGCAGTACCGGCTAGCAGGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCCACCGCATCAg                                                                                                                                     >  2:1394297/1‑150 (MQ=255)
                   atcaCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAgg                                                                                                                                    >  2:2049443/1‑150 (MQ=255)
                     caCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTc                                                                                                                                  <  2:4432732/150‑1 (MQ=255)
                      aCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTcc                                                                                                                                 <  1:329825/150‑1 (MQ=255)
                       ccTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCa                                                                                                                                >  2:159122/1‑150 (MQ=255)
                             ctGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGt                                                                                                                          >  1:3143762/1‑150 (MQ=255)
                             ctGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGt                                                                                                                          >  2:2496632/1‑150 (MQ=255)
                             ctGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGt                                                                                                                          >  1:829972/1‑150 (MQ=255)
                             ctGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGt                                                                                                                          <  1:584329/150‑1 (MQ=255)
                                 tgtgCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGtgatt                                                                                                                      >  1:2489608/1‑150 (MQ=255)
                                   tgCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGtgattga                                                                                                                    <  1:2209809/150‑1 (MQ=255)
                                      cTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGatat                                                                                                                 <  2:4009807/150‑1 (MQ=255)
                                      cTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGATAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGatat                                                                                                                 >  2:658256/1‑150 (MQ=255)
                                          cgcCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAg                                                                                                             >  1:4622850/1‑150 (MQ=255)
                                           gcCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGa                                                                                                            <  2:3834169/150‑1 (MQ=255)
                                             cATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCGAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATa                                                                                                          <  2:4503262/150‑1 (MQ=255)
                                              aTAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAg                                                                                                         <  1:3395193/150‑1 (MQ=255)
                                               tAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAgg                                                                                                        <  1:4061095/150‑1 (MQ=255)
                                                aaCATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAggt                                                                                                       >  1:3995715/1‑150 (MQ=255)
                                                aaCATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAggt                                                                                                       >  1:2321467/1‑150 (MQ=255)
                                                 tcATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAggtg                                                                                                      >  2:1887873/2‑150 (MQ=255)
                                                   aTGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAggtggt                                                                                                    >  1:4369916/1‑150 (MQ=255)
                                                    tGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGtt                                                                                                   >  1:1086024/1‑150 (MQ=255)
                                                    tGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGtt                                                                                                   <  1:1189708/150‑1 (MQ=255)
                                                     gTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTg                                                                                                  <  1:980242/150‑1 (MQ=255)
                                                      tCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTgg                                                                                                 >  1:1039911/1‑150 (MQ=255)
                                                        cGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCt                                                                                               <  2:783532/150‑1 (MQ=255)
                                                        cGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCt                                                                                               >  1:3169975/1‑150 (MQ=255)
                                                          tttGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtca                                                                                             <  1:149481/150‑1 (MQ=255)
                                                          tttGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtca                                                                                             <  1:940473/150‑1 (MQ=255)
                                                           ttGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcat                                                                                            <  2:2827591/150‑1 (MQ=255)
                                                           ttGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcat                                                                                            <  1:373551/150‑1 (MQ=255)
                                                            tGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcatc                                                                                           <  2:872877/150‑1 (MQ=255)
                                                               cATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCGACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAg                                                                                        >  1:2861344/1‑150 (MQ=255)
                                                                aTTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGt                                                                                       <  2:3494278/150‑1 (MQ=255)
                                                                     caccaGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATAc                                                                                  >  1:2029298/1‑150 (MQ=255)
                                                                      accaGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACa                                                                                 >  1:1457379/1‑150 (MQ=255)
                                                                      accaGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACa                                                                                 <  1:1243654/150‑1 (MQ=255)
                                                                      accaGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACa                                                                                 <  1:3758071/150‑1 (MQ=255)
                                                                        caGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAAc                                                                               <  2:2208333/150‑1 (MQ=255)
                                                                         aGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACg                                                                              >  2:2860315/1‑150 (MQ=255)
                                                                         aGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACg                                                                              >  1:1157480/1‑150 (MQ=255)
                                                                             gATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCg                                                                          >  2:3064730/1‑150 (MQ=255)
                                                                              aTCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGt                                                                         <  1:1958425/150‑1 (MQ=255)
                                                                               tCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACATCGACCACCGTTTTCCCCAGAGTCCTGAGTTAGCCCAACAACCGCATCAGGTCCACCAGGTGATAGATATCAAGATAGGAGGTTGGCTCATCAATTAATACAACGGGCgaa                                                                        >  1:3207074/1‑148 (MQ=255)
                                                                                cGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTAt                                                                       >  1:3255675/1‑150 (MQ=255)
                                                                                cGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTAt                                                                       <  2:782100/150‑1 (MQ=255)
                                                                                 gCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTAtt                                                                      <  2:3249607/150‑1 (MQ=255)
                                                                                 gCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTAtt                                                                      <  1:4755407/150‑1 (MQ=255)
                                                                                  cAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTc                                                                     >  2:3414976/1‑150 (MQ=255)
                                                                                    gTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTg                                                                   <  2:2905024/150‑1 (MQ=255)
                                                                                         ggCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCa                                                              <  1:3801808/150‑1 (MQ=255)
                                                                                           cTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAgg                                                            <  2:1727571/150‑1 (MQ=255)
                                                                                            tAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGa                                                           >  2:2445970/1‑150 (MQ=255)
                                                                                             aGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGAc                                                          >  1:4588756/1‑150 (MQ=255)
                                                                                               ccTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCa                                                        >  1:171078/1‑150 (MQ=255)
                                                                                                 tGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCAtt                                                      >  2:2293319/1‑150 (MQ=255)
                                                                                                 tGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCAtt                                                      <  2:1160743/150‑1 (MQ=255)
                                                                                                     tAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGcta                                                  <  2:2114593/150‑3 (MQ=255)
                                                                                                     tAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCa                                                  <  2:2664131/150‑1 (MQ=255)
                                                                                                       aGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGa                                                >  1:3023998/1‑150 (MQ=255)
                                                                                                          tCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAAt                                             <  1:4421444/150‑1 (MQ=255)
                                                                                                             tGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATgcg                                          >  2:1346928/1‑150 (MQ=255)
                                                                                                                agcaCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATgcgcgc                                       >  2:4598017/1‑150 (MQ=255)
                                                                                                                agcaCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAAAAATgcgcgc                                       >  2:275846/1‑150 (MQ=255)
                                                                                                                agcaCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAAAAATgcgcgc                                       >  2:2711952/1‑150 (MQ=255)
                                                                                                                   acaGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctgg                                    <  1:661418/150‑1 (MQ=255)
                                                                                                                    caGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctggc                                   <  2:4315263/150‑1 (MQ=255)
                                                                                                                    caGCGACCACCGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctggc                                   >  1:2127350/1‑150 (MQ=255)
                                                                                                                     aGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAGATGCgcgctggcg                                  <  1:3481739/150‑1 (MQ=255)
                                                                                                                     aGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctggcg                                  >  2:464397/1‑150 (MQ=255)
                                                                                                                        gACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctggcgctg                               <  2:39874/150‑1 (MQ=255)
                                                                                                                         accaccGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGa                              <  1:345110/150‑1 (MQ=255)
                                                                                                                              cGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGAccgcc                         >  1:510809/1‑150 (MQ=255)
                                                                                                                              cGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGAccgcc                         >  2:4869446/1‑150 (MQ=255)
                                                                                                                              cGTTTTCCCCTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGAccgcc                         >  2:2052515/1‑150 (MQ=255)
                                                                                                                                ttttCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCgg                       >  2:4245007/1‑150 (MQ=255)
                                                                                                                                    ccccTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAg                   <  1:1883330/150‑1 (MQ=255)
                                                                                                                                      ccTGAGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCt                 >  2:1211521/1‑150 (MQ=255)
                                                                                                                                          aGTCCGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGt             >  1:424488/1‑150 (MQ=255)
                                                                                                                                              cGGAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTCAc         >  2:1244902/1‑150 (MQ=255)
                                                                                                                                               ggAGTTAGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAAcc        >  1:2905082/1‑150 (MQ=255)
                                                                                                                                                gAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACcg       >  2:1086917/1‑150 (MQ=255)
                                                                                                                                                 aGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACcga      >  2:1779908/1‑150 (MQ=255)
                                                                                                                                                     cGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGaa  <  2:2434828/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAA  >  CP009273/4500579‑4500877

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