Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,287,829 C→A 46.8% S185S (TCG→TCT ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,287,8290CA46.8% 167.9 / inf 288S185S (TCG→TCTccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (69/84);  new base A (59/76);  total (128/160)
Fisher's exact test for biased strand distribution p-value = 9.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

CCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAA  >  CP009273/2287680‑2287977
                                                                                                                                                     |                                                                                                                                                    
ccGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCAcc                                                                                                                                                      <  2:4215823/150‑1 (MQ=255)
  gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                                    <  2:1457355/150‑1 (MQ=255)
  gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                                    >  2:1222718/1‑150 (MQ=255)
  gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                                    >  2:1222387/1‑150 (MQ=255)
  gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                                    >  1:3198888/1‑150 (MQ=255)
    gTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAATATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAAAGGCAAAAGCAAACGCCACCgaga                                                                                                                                                  >  1:2786658/1‑150 (MQ=255)
      tCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGaaa                                                                                                                                                <  2:4328636/150‑1 (MQ=255)
         tAATCGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCTGGAAGATCCACGCTGCCAGCGTCCCCGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAccc                                                                                                                                             <  1:4728570/150‑1 (MQ=255)
         tAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAAccc                                                                                                                                             >  1:4137805/1‑150 (MQ=255)
           aTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCAc                                                                                                                                           <  1:1147170/150‑1 (MQ=255)
                 ccATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAccc                                                                                                                                     >  2:114468/1‑150 (MQ=255)
                 ccATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAccc                                                                                                                                     <  2:4486026/150‑1 (MQ=255)
                      cgtaaCCGAGCACGTTGCCGAGCTTCTGTAATATCCACGCGGCCAGCGTCCACGGACTCGAAAAACGGGAATACGTGCTGTCCAGACGCCCGATCAGAAAAGAAGAAATGGCAAAAGAAAACGCCACCGAGAAACCCACGTACCCCatat                                                                                                                                <  1:1656642/147‑1 (MQ=255)
                      cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGATTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCatat                                                                                                                                <  2:858552/150‑1 (MQ=255)
                       ggAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGTCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCAtata                                                                                                                               >  1:1004175/1‑150 (MQ=255)
                        gAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATaa                                                                                                                              >  2:53162/1‑150 (MQ=255)
                          aCCGAGCACGTTGCCGAGGGTCTGGAAGATCCAGGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCATTGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAg                                                                                                                            <  1:2205009/150‑1 (MQ=255)
                          aCCGAGCACGATGCGGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAg                                                                                                                            <  2:2606261/150‑1 (MQ=255)
                          aCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAg                                                                                                                            >  2:4085114/1‑150 (MQ=255)
                          aCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAg                                                                                                                            >  1:4351902/1‑150 (MQ=255)
                          aCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGAAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAg                                                                                                                            <  2:1568153/150‑1 (MQ=255)
                             gAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCag                                                                                                                         <  2:3472456/150‑1 (MQ=255)
                                cACGATGCCGAGCGTCAGTAATATCCACGCTGCCAGCGTCCACGGACTATTAAAACGCGAATAAGTGCTGTCCATACGCCCGCTCAGCAAAGAAGCAATGGAAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAgg                                                                                                                      <  2:3648401/150‑1 (MQ=255)
                                  cGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCg                                                                                                                    <  1:2749247/150‑1 (MQ=255)
                                    aTGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGATCAGAAAAGAAGAAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGa                                                                                                                  <  1:1865737/150‑1 (MQ=255)
                                     tGCCGAGCGTCAGGACGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGAt                                                                                                                 <  2:4734861/150‑1 (MQ=255)
                                     tGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGAt                                                                                                                 <  2:4028600/150‑1 (MQ=255)
                                      gCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATg                                                                                                                <  2:1789081/150‑1 (MQ=255)
                                        cGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGa                                                                                                              >  2:2181351/1‑150 (MQ=255)
                                         gAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGaa                                                                                                             >  1:1459219/1‑150 (MQ=255)
                                          aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAg                                                                                                            >  1:3456353/1‑150 (MQ=255)
                                            cGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                                                                          <  1:1018606/150‑1 (MQ=255)
                                             gTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                         <  1:1410941/150‑1 (MQ=255)
                                                    ataTCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCg                                                                                                  >  1:3737468/3‑150 (MQ=255)
                                                    ataTCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCg                                                                                                  >  1:3739271/3‑150 (MQ=255)
                                                        ccACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATc                                                                                              >  1:3984136/1‑150 (MQ=255)
                                                          aCGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCt                                                                                            >  1:1766672/1‑150 (MQ=255)
                                                             cTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCa                                                                                         <  1:1178769/150‑1 (MQ=255)
                                                                ccAGCGTCCACGGACGATTAAAACGCGCATAAGTGCTGTCCACACGCCCGCTCACCAAAGAAGCATTGGAAAAACCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATa                                                                                      <  1:1214285/150‑1 (MQ=255)
                                                                ccAGCGTCCACGGACGAGTAACCCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATa                                                                                      <  2:758617/150‑1 (MQ=255)
                                                                ccAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATa                                                                                      >  1:3601382/1‑150 (MQ=255)
                                                                  aGCGTCCACGGACGAGTAAACCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGc                                                                                    <  2:817030/150‑1 (MQ=255)
                                                                  aGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGc                                                                                    <  1:4374235/150‑1 (MQ=255)
                                                                    cGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCgg                                                                                  <  2:3066492/150‑1 (MQ=255)
                                                                         aCGGACGAGTAAACCGCGCATAAGTGCTGCCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTaa                                                                             <  2:3418692/150‑1 (MQ=255)
                                                                         aCGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTaa                                                                             <  2:2347152/150‑1 (MQ=255)
                                                                           ggACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGa                                                                           <  2:3346370/150‑1 (MQ=255)
                                                                              cGAGTAAATCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAATCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcg                                                                        <  1:3513654/150‑1 (MQ=255)
                                                                              cGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcg                                                                        <  2:1378020/150‑1 (MQ=255)
                                                                               gAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcgc                                                                       >  2:3983077/1‑150 (MQ=255)
                                                                                  tAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                                    >  2:1177616/1‑150 (MQ=255)
                                                                                    aaaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCt                                                                  <  2:331788/150‑1 (MQ=255)
                                                                                    aaaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCt                                                                  >  2:3879869/1‑150 (MQ=255)
                                                                                     aaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                                 >  1:3899852/1‑150 (MQ=255)
                                                                                       cgcgCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCa                                                               >  1:3249877/1‑150 (MQ=255)
                                                                                         cgcATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAAt                                                             >  1:156192/1‑150 (MQ=255)
                                                                                          tcATATGTGCTGTCCAGACGCCCGCTCATCATCTAATAATTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATc                                                            <  1:1326928/149‑1 (MQ=255)
                                                                                          gcATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATc                                                            <  1:1144648/150‑1 (MQ=255)
                                                                                             tAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCggg                                                         >  2:4204197/1‑150 (MQ=255)
                                                                                                   cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                   <  2:1710784/150‑1 (MQ=255)
                                                                                                   cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                   <  2:2774344/150‑1 (MQ=255)
                                                                                                   cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                   <  1:3121683/150‑1 (MQ=255)
                                                                                                   cTGTCCAAACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                   <  2:3967841/150‑1 (MQ=255)
                                                                                                       ccAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCa                                               >  2:2206777/1‑150 (MQ=255)
                                                                                                        cAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                              >  2:205350/1‑150 (MQ=255)
                                                                                                        cAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                              >  2:2270699/1‑150 (MQ=255)
                                                                                                        cAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                              <  1:4166145/150‑1 (MQ=255)
                                                                                                        cAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                              >  2:4443646/1‑150 (MQ=255)
                                                                                                         aGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAAc                                             >  1:639844/1‑150 (MQ=255)
                                                                                                            cGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTg                                          <  1:2081459/150‑1 (MQ=255)
                                                                                                             gCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCACTAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCACCGTGc                                         >  2:1020792/1‑150 (MQ=255)
                                                                                                             gCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCACTAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGc                                         >  1:3282893/1‑150 (MQ=255)
                                                                                                             gCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGc                                         >  1:2468636/1‑150 (MQ=255)
                                                                                                                cGCTCAGCAAAGAAGCAATGGAAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgag                                      <  1:1999113/150‑1 (MQ=255)
                                                                                                                 gCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgaga                                     <  2:3793455/150‑1 (MQ=255)
                                                                                                                  cTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                                    <  2:1127261/150‑1 (MQ=255)
                                                                                                                    cAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                                  >  2:1465504/1‑150 (MQ=255)
                                                                                                                     aGCAAAGAAGCAATCGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                                 <  2:4626932/150‑1 (MQ=255)
                                                                                                                      gCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACg                                <  2:3244105/150‑1 (MQ=255)
                                                                                                                        aaaGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACggg                              <  2:2811468/150‑1 (MQ=255)
                                                                                                                          agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGtt                            <  2:2543484/150‑1 (MQ=255)
                                                                                                                          agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGtt                            <  1:4193152/150‑1 (MQ=255)
                                                                                                                           gaagCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTa                           >  1:2653404/1‑150 (MQ=255)
                                                                                                                             agCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTCCCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTaga                         <  1:1496466/150‑1 (MQ=255)
                                                                                                                                aaTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                      <  1:1042383/150‑1 (MQ=255)
                                                                                                                                aaTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                      >  2:3750341/1‑150 (MQ=255)
                                                                                                                                 gtGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTa                     >  2:3137631/2‑150 (MQ=255)
                                                                                                                                  tGGCAAAAGCACACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaa                    <  2:2867740/150‑1 (MQ=255)
                                                                                                                                  tGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaa                    <  1:946380/150‑1 (MQ=255)
                                                                                                                                    gCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAag                  <  1:4513744/150‑1 (MQ=255)
                                                                                                                                    gCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAag                  >  2:662460/1‑150 (MQ=255)
                                                                                                                                    gCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAag                  >  2:3546582/1‑150 (MQ=255)
                                                                                                                                      aaaaGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAagag                >  1:482377/1‑150 (MQ=255)
                                                                                                                                         aGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAAt             <  2:4034709/150‑1 (MQ=255)
                                                                                                                                          gCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATg            <  1:4723752/150‑1 (MQ=255)
                                                                                                                                           cAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGa           >  1:2998212/1‑150 (MQ=255)
                                                                                                                                            aaaCGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGaa          >  1:3639813/1‑150 (MQ=255)
                                                                                                                                             aaCGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAAc         <  2:1948028/150‑1 (MQ=255)
                                                                                                                                              aCGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACa        <  1:1792730/150‑1 (MQ=255)
                                                                                                                                                    ccGAGAAACCCACGCACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCaaaa  <  1:1679272/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAA  >  CP009273/2287680‑2287977

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: