Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,460,036 T→A 5.2% N17Y (AAT→TAT)  bfr ← bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,460,0360TA 5.2% 1069.6 / 37.2 307N17Y (AAT→TAT) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base T (133/156);  new base A (8/8);  total (141/166)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 313 total aligned reads displayed.

TTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTG  >  CP009273/3459887‑3460185
                                                                                                                                                     |                                                                                                                                                     
ttCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCAtt                                                                                                                                                       <  1:1177532/150‑1 (MQ=255)
ttCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCAtt                                                                                                                                                       >  1:3786911/1‑150 (MQ=255)
  ccAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTc                                                                                                                                                     <  1:4938003/150‑1 (MQ=255)
    aGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTccc                                                                                                                                                   <  2:699471/150‑1 (MQ=255)
     gAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCACTGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCa                                                                                                                                                  <  2:104721/150‑1 (MQ=255)
      aaaaaGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCaa                                                                                                                                                 >  1:2328533/1‑150 (MQ=255)
       aaaaGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAAc                                                                                                                                                >  1:3038039/1‑150 (MQ=255)
       aaaaGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAAc                                                                                                                                                <  1:4881169/150‑1 (MQ=255)
        aaaGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATATCCCAACa                                                                                                                                               >  2:3076475/1‑150 (MQ=255)
         aaGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAg                                                                                                                                              <  2:23494/150‑1 (MQ=255)
            aaTGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGttt                                                                                                                                           <  1:3250060/150‑1 (MQ=255)
             aTGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTttg                                                                                                                                          >  2:3370183/1‑150 (MQ=255)
                 gcTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTga                                                                                                                                      >  2:1316418/1‑150 (MQ=255)
                    cAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGAt                                                                                                                                   >  1:1772498/1‑150 (MQ=255)
                    cAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGAt                                                                                                                                   <  1:539606/150‑1 (MQ=255)
                        ataACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAAtt                                                                                                                               >  1:1982298/1‑150 (MQ=255)
                          aaCGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTta                                                                                                                             <  1:1664500/150‑1 (MQ=255)
                            cGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTtata                                                                                                                           >  2:1120964/1‑150 (MQ=255)
                            cGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTtata                                                                                                                           <  1:3692597/150‑1 (MQ=255)
                               tCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACt                                                                                                                        <  2:1451548/150‑1 (MQ=255)
                               tCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACt                                                                                                                        <  1:2404551/150‑1 (MQ=255)
                                cGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACtt                                                                                                                       <  2:4374184/150‑1 (MQ=255)
                                 ggCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACttt                                                                                                                      >  2:642998/1‑150 (MQ=255)
                                 ggCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACttt                                                                                                                      >  2:3264959/1‑150 (MQ=255)
                                  gCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATATCCCAACAGTTTGTTGAGATAATTTATAACTTTa                                                                                                                     <  1:794792/150‑1 (MQ=255)
                                      gtTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTAt                                                                                                                 <  1:4666207/150‑1 (MQ=255)
                                          catctcatcAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCAcc                                                                                                             <  1:4017678/150‑1 (MQ=255)
                                          catctcatcAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTATGTTGAGATAATTTATAACTTTAGTATAAcc                                                                                                             >  1:3811178/1‑150 (MQ=255)
                                            tctcatcAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCtt                                                                                                           <  2:205791/150‑1 (MQ=255)
                                              tcatcaATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTc                                                                                                         >  1:3314772/1‑150 (MQ=255)
                                                atcaATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCAt                                                                                                       <  1:3568687/150‑1 (MQ=255)
                                                   aaTGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCAtttt                                                                                                    <  2:4157658/150‑1 (MQ=255)
                                                      ggATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTg                                                                                                 >  1:4177572/1‑150 (MQ=255)
                                                      ggATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTg                                                                                                 >  1:1338454/1‑150 (MQ=255)
                                                        aTTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGAc                                                                                               >  2:1603375/1‑150 (MQ=255)
                                                           cATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTcc                                                                                            >  1:4851845/1‑150 (MQ=255)
                                                           cATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTcc                                                                                            <  2:1534716/150‑1 (MQ=255)
                                                              gATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTc                                                                                         <  2:2723765/150‑1 (MQ=255)
                                                              gATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTc                                                                                         <  2:2209010/150‑1 (MQ=255)
                                                              gATACTACACATCATTGAGACGTTTGAGACCCCAGTTATTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTc                                                                                         <  2:3803834/150‑1 (MQ=255)
                                                                  cTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCtt                                                                                     >  2:1820755/1‑150 (MQ=255)
                                                                    ccACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTcc                                                                                   <  2:4486480/150‑1 (MQ=255)
                                                                     cacaTCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCa                                                                                  <  2:5604141/150‑1 (MQ=255)
                                                                      acaTCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTATTGACTCCTTCGCTTCCAc                                                                                 >  2:4700750/1‑150 (MQ=255)
                                                                       catcatTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCAct                                                                                >  2:1376962/1‑150 (MQ=255)
                                                                        atcatTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCActc                                                                               >  1:2318844/1‑150 (MQ=255)
                                                                        atcatTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCActc                                                                               >  1:4570435/1‑150 (MQ=255)
                                                                          catTGAGACGTTTGAGACCCCAGTTTTTAAACATTCTGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCtt                                                                             >  1:5386573/1‑150 (MQ=255)
                                                                            ttGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTg                                                                           <  1:2368613/150‑1 (MQ=255)
                                                                              gagaCGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGaa                                                                         <  1:4885940/150‑1 (MQ=255)
                                                                                  cGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGt                                                                     <  1:2251320/150‑1 (MQ=255)
                                                                                  cGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGt                                                                     <  1:3817784/150‑1 (MQ=255)
                                                                                    tttGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAg                                                                   <  1:2199039/150‑1 (MQ=255)
                                                                                        agaCCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATATTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACg                                                               >  2:515025/1‑150 (MQ=255)
                                                                                        agaCCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATATCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACg                                                               >  1:3988830/1‑150 (MQ=255)
                                                                                          aCCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGc                                                             >  2:2190956/1‑150 (MQ=255)
                                                                                            cccAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCtt                                                           <  2:4206611/150‑1 (MQ=255)
                                                                                             ccAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTa                                                          >  1:1854334/1‑150 (MQ=255)
                                                                                             ccAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTa                                                          >  1:3804126/1‑150 (MQ=255)
                                                                                              cAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTAc                                                         <  2:2652636/150‑1 (MQ=255)
                                                                                              cAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTAc                                                         <  1:5115110/150‑1 (MQ=255)
                                                                                              cAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTAc                                                         <  1:4306794/150‑1 (MQ=255)
                                                                                               aGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGACCGGCTTACg                                                        >  1:3475737/1‑150 (MQ=255)
                                                                                                  ttttAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGt                                                     >  2:4965686/1‑150 (MQ=255)
                                                                                                  ttttAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGt                                                     >  1:90746/1‑150 (MQ=255)
                                                                                                    ttAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTaa                                                   <  1:3264959/150‑1 (MQ=255)
                                                                                                    ttAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTaa                                                   >  2:4199974/1‑150 (MQ=255)
                                                                                                      aaaCATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGt                                                 >  2:2413174/1‑150 (MQ=255)
                                                                                                       aaCATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTc                                                >  1:529827/1‑150 (MQ=255)
                                                                                                         cATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCaa                                              >  1:3388016/1‑150 (MQ=255)
                                                                                                              gggCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAac                                         >  2:116090/1‑150 (MQ=255)
                                                                                                              gggCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAac                                         <  2:3810885/150‑1 (MQ=255)
                                                                                                               ggCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAaca                                        >  1:283052/1‑150 (MQ=255)
                                                                                                               ggCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATATCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAaca                                        <  1:3076475/150‑1 (MQ=255)
                                                                                                               ggCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATATCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAaca                                        <  1:213878/150‑1 (MQ=255)
                                                                                                                  aTGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAacacac                                     >  1:5051259/1‑150 (MQ=255)
                                                                                                                     gagaAAGTACTGATTGATTGCGACAAGCTCATGTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCCCCTTTCATATTATGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTcc                                  <  2:3914082/150‑1 (MQ=255)
                                                                                                                      agaAAGTACTGATTGATTGCGACAAGCTCATATCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCg                                 >  2:330979/1‑150 (MQ=255)
                                                                                                                        aaaGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGtt                               <  2:2119142/150‑1 (MQ=255)
                                                                                                                         aaGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGttt                              >  2:2962555/1‑150 (MQ=255)
                                                                                                                           gTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTAt                            <  2:3595054/150‑1 (MQ=255)
                                                                                                                           gTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTAt                            <  1:5477655/150‑1 (MQ=255)
                                                                                                                           gTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCATTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTAt                            <  2:2761038/150‑1 (MQ=255)
                                                                                                                              cTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCg                         >  1:969718/1‑150 (MQ=255)
                                                                                                                                   tgattgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTGTATGCGGACTc                    >  1:651048/1‑150 (MQ=255)
                                                                                                                                    gattgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAATCACACTCCGTTTATGCGGACTcc                   >  2:819651/1‑150 (MQ=255)
                                                                                                                                    gattgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTcc                   >  2:3667056/1‑150 (MQ=255)
                                                                                                                                    gattgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTcc                   >  1:1480259/1‑150 (MQ=255)
                                                                                                                                      ttgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCtt                 <  1:5588629/150‑1 (MQ=255)
                                                                                                                                      ttgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCtt                 >  2:4775618/1‑150 (MQ=255)
                                                                                                                                       tgCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTa                >  1:4844510/1‑150 (MQ=255)
                                                                                                                                         cGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTaaa              >  2:5236033/1‑150 (MQ=255)
                                                                                                                                           aCAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACt            <  1:1015743/150‑1 (MQ=255)
                                                                                                                                           aCAAGCTCATGTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACt            <  1:4924563/150‑1 (MQ=255)
                                                                                                                                            cAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTc           >  2:2817043/1‑150 (MQ=255)
                                                                                                                                             aaGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTcc          >  1:4993712/1‑150 (MQ=255)
                                                                                                                                                cTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGc       >  2:3758995/1‑150 (MQ=255)
                                                                                                                                                 tCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGAGGACTCCTTAAACTCCGGCa      >  2:4295287/1‑150 (MQ=255)
                                                                                                                                                     ttCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTGGGGGCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTg  >  1:440045/1‑150 (MQ=255)
                                                                                                                                                     ttCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTg  >  1:4819004/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTG  >  CP009273/3459887‑3460185

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: