Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 743,797 G→T 100% R214R (CGC→CGA ybgO ← putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273743,7970GT97.8% 719.0 / 4.6 225R214R (CGC→CGAybgOputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (98/120);  total (101/124)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 226 total aligned reads displayed.

GGTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGAC  >  CP009273/743650‑743943
                                                                                                                                                   |                                                                                                                                                  
    cATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATg                                                                                                                                              >  2:1682156/1‑150 (MQ=255)
    cATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATg                                                                                                                                              <  1:1789018/150‑1 (MQ=255)
     aTTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGt                                                                                                                                             >  2:726451/1‑150 (MQ=255)
     aTTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGt                                                                                                                                             <  2:885357/150‑1 (MQ=255)
      ttGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtt                                                                                                                                            >  1:112481/1‑150 (MQ=255)
        gATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtttt                                                                                                                                          <  1:1457109/150‑1 (MQ=255)
        gATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtttt                                                                                                                                          <  2:2334185/150‑1 (MQ=255)
         aTCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTc                                                                                                                                         <  2:2013992/150‑1 (MQ=255)
         aTCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTc                                                                                                                                         >  2:1256225/1‑150 (MQ=255)
           cGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCa                                                                                                                                       >  2:1737500/1‑150 (MQ=255)
                 ggCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAg                                                                                                                                 <  2:738522/150‑1 (MQ=255)
                     ctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAat                                                                                                                             <  1:2150255/150‑1 (MQ=255)
                        ggTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATg                                                                                                                          <  1:1664889/150‑1 (MQ=255)
                         gTCGTTTTAATACTGAACGTTTTCGACATTTTGTGGCGACGTATATCCACCAAATTAAATTTACCAAAATCAATCAACTGTGTTTCCGGTGAAATGGAAAAAATTACTATGCAATCCAGCACACGAATCTTTTCCAGACCAGTAATATGa                                                                                                                         >  1:4109774/1‑150 (MQ=255)
                          tCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGat                                                                                                                        <  2:3004921/150‑1 (MQ=255)
                            gTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGatat                                                                                                                      >  1:3363129/1‑150 (MQ=255)
                                tAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCa                                                                                                                  >  2:631571/1‑150 (MQ=255)
                                 aaTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAg                                                                                                                 >  2:1494425/1‑150 (MQ=255)
                                 aaTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAg                                                                                                                 >  1:2658485/1‑150 (MQ=255)
                                 aaTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTGTCCAGACCAGTAATATGATAGTTCAg                                                                                                                 >  2:4734609/1‑150 (MQ=255)
                                   tACTGAACCTTAACTACAATGACTGGCGACGTATATCCAGCAAATAAAATTTACACAAATCAATCACCTGTGTTTCCGGTGAAATGCAAAAATTGACGCTGCAATCCAGCACACGAATGTTTTCCAGACCAGTAATATGATAGTTCAGGt                                                                                                               >  2:470773/1‑150 (MQ=255)
                                      tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtttt                                                                                                            <  2:904317/150‑1 (MQ=255)
                                      tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtttt                                                                                                            >  2:3320796/1‑150 (MQ=255)
                                      tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtttt                                                                                                            >  1:831624/1‑150 (MQ=255)
                                       gAACGTTTTCGACATTGTGTTGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCAGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                           <  1:3219951/150‑1 (MQ=255)
                                       gAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                           <  1:2594807/150‑1 (MQ=255)
                                          cGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGc                                                                                                        <  1:390107/150‑1 (MQ=255)
                                           gTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCa                                                                                                       >  1:2314526/1‑150 (MQ=255)
                                                cGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAgg                                                                                                  >  2:3653502/1‑150 (MQ=255)
                                                 taCATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAggg                                                                                                 >  1:1912338/2‑150 (MQ=255)
                                                  aCATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGt                                                                                                <  2:442209/150‑1 (MQ=255)
                                                  aCATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGt                                                                                                >  2:2587677/1‑150 (MQ=255)
                                                         gtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGtt                                                                                         <  1:220197/150‑1 (MQ=255)
                                                          tgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTa                                                                                        <  1:4142473/150‑1 (MQ=255)
                                                           ggCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTaa                                                                                       <  1:2574568/150‑1 (MQ=255)
                                                            gCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAAc                                                                                      <  2:2181808/150‑1 (MQ=255)
                                                            gCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAAc                                                                                      <  2:4889527/150‑1 (MQ=255)
                                                             cGACGTATATCCAGCAAATTAAATTTACCAAATTCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGAAGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACg                                                                                     <  1:2128551/150‑1 (MQ=255)
                                                             cGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACg                                                                                     <  1:4339518/150‑1 (MQ=255)
                                                               aCGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCt                                                                                   <  2:2078113/150‑1 (MQ=255)
                                                               aCGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCt                                                                                   <  1:2777976/150‑1 (MQ=255)
                                                                cGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                  >  2:3987735/1‑150 (MQ=255)
                                                                cGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                  >  2:4322075/1‑150 (MQ=255)
                                                                cGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                  <  1:5094995/150‑1 (MQ=255)
                                                                  tataTCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTAcc                                                                                <  1:744918/150‑1 (MQ=255)
                                                                  tataTCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTAcc                                                                                <  2:537249/150‑1 (MQ=255)
                                                                  tataTCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTAcc                                                                                >  1:2857396/1‑150 (MQ=255)
                                                                   atatCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGTTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCg                                                                               >  2:2945365/1‑150 (MQ=255)
                                                                    tatCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCgg                                                                              <  1:1342164/150‑1 (MQ=255)
                                                                    tatCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCgg                                                                              >  1:4734696/1‑150 (MQ=255)
                                                                     atCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGc                                                                             >  1:3626445/1‑150 (MQ=255)
                                                                            aaattaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCa                                                                      <  1:1935507/150‑1 (MQ=255)
                                                                                taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAAtt                                                                  <  1:1597620/150‑1 (MQ=255)
                                                                                taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAAtt                                                                  <  2:3243570/150‑1 (MQ=255)
                                                                                taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAAtt                                                                  <  1:2048236/150‑1 (MQ=255)
                                                                                        ccAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacgacg                                                          <  2:2854248/150‑1 (MQ=255)
                                                                                             atcaatcaCCTGTGTTTCAGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGAtata                                                     >  2:1222648/1‑150 (MQ=255)
                                                                                                aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTc                                                  <  2:2674960/150‑1 (MQ=255)
                                                                                                aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTc                                                  <  1:2357449/150‑1 (MQ=255)
                                                                                                 atcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCg                                                 >  1:3852868/1‑150 (MQ=255)
                                                                                                   caCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCt                                               <  2:662355/150‑1 (MQ=255)
                                                                                                      cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAg                                            <  1:1979921/150‑1 (MQ=255)
                                                                                                      cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAg                                            >  1:4968102/1‑150 (MQ=255)
                                                                                                        gtgtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCg                                          >  2:4313028/1‑150 (MQ=255)
                                                                                                         tgtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGt                                         >  2:2320970/1‑150 (MQ=255)
                                                                                                         tgtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGt                                         <  2:157056/150‑1 (MQ=255)
                                                                                                          gtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTa                                        >  1:3745779/1‑150 (MQ=255)
                                                                                                          gtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTa                                        <  2:2809655/150‑1 (MQ=255)
                                                                                                            ttCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACt                                      <  1:4391988/150‑1 (MQ=255)
                                                                                                              ccGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACttt                                    <  1:1836511/150‑1 (MQ=255)
                                                                                                              ccGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACttt                                    >  1:3976420/1‑150 (MQ=255)
                                                                                                                    aaaTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAAt                              >  1:3955481/1‑150 (MQ=255)
                                                                                                                     aaTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATc                             <  2:4109774/150‑1 (MQ=255)
                                                                                                                      aTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCa                            <  1:4749682/150‑1 (MQ=255)
                                                                                                                       tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCAt                           <  2:1101100/150‑1 (MQ=255)
                                                                                                                        ggAAAAATTGACGCTGCAATCCAGCACTCGATTGTTTTCCAGACCAGTATTATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGATACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATTATCATg                          <  1:2966172/150‑1 (MQ=255)
                                                                                                                        ggAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATg                          <  1:2419158/150‑1 (MQ=255)
                                                                                                                          aaaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcg                        <  2:2691719/150‑1 (MQ=255)
                                                                                                                          aaaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcg                        <  2:106087/150‑1 (MQ=255)
                                                                                                                           aaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcgg                       >  1:2698102/1‑150 (MQ=255)
                                                                                                                           aaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcgg                       >  2:309247/1‑150 (MQ=255)
                                                                                                                           aaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcgg                       <  2:3613450/150‑1 (MQ=255)
                                                                                                                             aaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcggcg                     <  1:726451/150‑1 (MQ=255)
                                                                                                                              aTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGcggcgg                    >  1:4896984/1‑150 (MQ=255)
                                                                                                                              aTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGcggcgg                    >  2:3210987/1‑150 (MQ=255)
                                                                                                                               ttGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCggcggc                   <  2:2570792/150‑1 (MQ=255)
                                                                                                                                tGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCa                  <  2:2060726/150‑1 (MQ=255)
                                                                                                                                 gACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                 >  2:3794745/1‑150 (MQ=255)
                                                                                                                                 gACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                 >  2:4618110/1‑150 (MQ=255)
                                                                                                                                   cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAttt               >  2:2414116/1‑150 (MQ=255)
                                                                                                                                    gCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtc              <  1:138259/150‑1 (MQ=255)
                                                                                                                                       gCAATCCAGGACTCGAATGTTTTCCAGTCCTGTAATATGATATTTCAGGTTTTGGGCAGTATGTTCTTCGTTAACTGTACAGGCACCGTCAAATTGCACGAGGATTTATTCGCTAAGCGTACTTTTATAATCATGCGGCGGCATTTCTCg           <  2:4279896/150‑1 (MQ=255)
                                                                                                                                       gCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCg           >  1:5002615/1‑150 (MQ=255)
                                                                                                                                         aaTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa         <  1:4593973/150‑1 (MQ=255)
                                                                                                                                           tCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTAAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGAGTACTATGATAATCATGCGGAGGCATTTCTCGTAtt       >  1:3108664/1‑150 (MQ=255)
                                                                                                                                            ccAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAttt      <  1:3142868/150‑1 (MQ=255)
                                                                                                                                             cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtttt     >  2:1981124/1‑150 (MQ=255)
                                                                                                                                             cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtttt     >  2:4066852/1‑150 (MQ=255)
                                                                                                                                              aGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTg    >  1:4666091/1‑150 (MQ=255)
                                                                                                                                                cACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGAc  >  2:3654136/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGAC  >  CP009273/743650‑743943

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: