Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 523,781 G→T 13.4% pseudogene (374/708 nt) rhsH → pseudogene, Rhs familly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273523,7810GT13.4% 950.8 / 78.3 285pseudogene (374/708 nt)rhsHpseudogene, Rhs familly protein
Reads supporting (aligned to +/- strand):  ref base G (129/116);  new base T (18/22);  total (147/138)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATG  >  CP009273/523632‑523929
                                                                                                                                                     |                                                                                                                                                    
  gTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGAt                                                                                                                                                    <  1:3012409/150‑1 (MQ=255)
    cTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTa                                                                                                                                                  >  2:2417042/1‑150 (MQ=255)
       tATGCGTATCCGCTGAATCCGGTGAATTTTATTGATCCATTAGGGTTAAGTCCCTCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAAt                                                                                                                                               <  2:854647/150‑1 (MQ=255)
          gCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTa                                                                                                                                            <  2:2673480/150‑1 (MQ=255)
            gTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAgg                                                                                                                                          >  1:3589943/1‑150 (MQ=255)
               tCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTgg                                                                                                                                       <  1:231010/150‑1 (MQ=255)
                 cGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGa                                                                                                                                     <  1:2631835/150‑1 (MQ=255)
                  gCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGAc                                                                                                                                    >  2:2590661/1‑150 (MQ=255)
                   cTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACt                                                                                                                                   <  1:474328/150‑1 (MQ=255)
                    tGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTg                                                                                                                                  <  1:1917752/150‑1 (MQ=255)
                      aaTCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGAt                                                                                                                                <  2:1978096/150‑1 (MQ=255)
                       aTCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATc                                                                                                                               <  2:2503363/150‑1 (MQ=255)
                       aTCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATc                                                                                                                               <  2:1260813/150‑1 (MQ=255)
                         ccGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCaa                                                                                                                             >  2:7218/1‑150 (MQ=255)
                           ggTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAAtt                                                                                                                           >  1:826862/1‑150 (MQ=255)
                               aaTGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTc                                                                                                                       <  2:1147361/150‑1 (MQ=255)
                                  ggTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCAt                                                                                                                    >  2:2295104/1‑150 (MQ=255)
                                   gTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCAtt                                                                                                                   <  2:3331255/150‑1 (MQ=255)
                                    tATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTg                                                                                                                  >  1:3499504/1‑150 (MQ=255)
                                    tATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTg                                                                                                                  >  2:1309800/1‑150 (MQ=255)
                                         aTCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGACGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATgg                                                                                                             <  2:3475070/150‑1 (MQ=255)
                                           ccATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCa                                                                                                           >  2:3685218/1‑150 (MQ=255)
                                            cATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCAt                                                                                                          >  1:2032214/1‑150 (MQ=255)
                                            cATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCAt                                                                                                          <  1:3220018/150‑1 (MQ=255)
                                             aTTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCAtt                                                                                                         >  1:3879610/1‑150 (MQ=255)
                                               tAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCAtttt                                                                                                       >  2:3282471/1‑150 (MQ=255)
                                               tAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCAtttt                                                                                                       >  2:1600299/1‑150 (MQ=255)
                                               tAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCAtttt                                                                                                       <  1:2182040/150‑1 (MQ=255)
                                                 gggTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGt                                                                                                     <  2:3113750/150‑1 (MQ=255)
                                                  ggTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTc                                                                                                    <  1:2701365/150‑1 (MQ=255)
                                                   gTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCg                                                                                                   >  1:672510/1‑150 (MQ=255)
                                                   gTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCg                                                                                                   >  2:3478570/1‑150 (MQ=255)
                                                    ttAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGa                                                                                                  <  2:1367843/150‑1 (MQ=255)
                                                    ttAAGTCCCGCAGATGTAGCGATAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGa                                                                                                  <  2:762932/150‑1 (MQ=255)
                                                      aaGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAgt                                                                                                <  2:1895591/150‑1 (MQ=255)
                                                      aaGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAgt                                                                                                <  2:500080/150‑1 (MQ=255)
                                                          cccGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAgtgtat                                                                                            >  2:3407232/1‑148 (MQ=255)
                                                          cccGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCt                                                                                            <  2:3319055/150‑1 (MQ=255)
                                                            cGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaa                                                                                          >  2:1790263/1‑150 (MQ=255)
                                                              cAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaaat                                                                                        >  1:3969030/1‑150 (MQ=255)
                                                                gaagTATCGCTAATAAGAAGAAAATATCAACAAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaaatta                                                                                      >  1:267178/4‑150 (MQ=255)
                                                                gATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaaatta                                                                                      >  1:266847/1‑150 (MQ=255)
                                                                 aTGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaaattaa                                                                                     >  2:3986363/1‑150 (MQ=255)
                                                                  tGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaaattaaa                                                                                    >  2:2934087/1‑150 (MQ=255)
                                                                  tGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCtaaattaaa                                                                                    <  1:2746183/150‑1 (MQ=255)
                                                                    tAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATg                                                                                  <  1:2265627/150‑1 (MQ=255)
                                                                       cgcgAATAAGAAGAAAAGATCAACTAACCCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACg                                                                               <  2:2994423/150‑1 (MQ=255)
                                                                       cgcTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACg                                                                               <  2:3732651/150‑1 (MQ=255)
                                                                             taaGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGgtg                                                                         >  1:126163/1‑150 (MQ=255)
                                                                               agaagaAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGtt                                                                       >  2:1862225/1‑150 (MQ=255)
                                                                                gaagaaAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTa                                                                      <  2:825288/150‑1 (MQ=255)
                                                                                  agaaAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGc                                                                    >  2:915333/1‑150 (MQ=255)
                                                                                   gaaAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGcc                                                                   >  2:285185/1‑150 (MQ=255)
                                                                                   gaaAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGcc                                                                   <  2:1379442/150‑1 (MQ=255)
                                                                                    aaaaGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCg                                                                  >  2:3774401/1‑150 (MQ=255)
                                                                                         aTCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCgg                                                             >  1:354505/1‑150 (MQ=255)
                                                                                             aCTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGaa                                                         >  2:58415/1‑150 (MQ=255)
                                                                                             aCTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGaa                                                         >  2:2966025/1‑150 (MQ=255)
                                                                                             aCTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGaa                                                         >  1:942003/1‑150 (MQ=255)
                                                                                              cTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGaaa                                                        <  1:19255/150‑1 (MQ=255)
                                                                                               tAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAg                                                       >  2:2008624/1‑150 (MQ=255)
                                                                                                aaaCCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAgg                                                      <  1:1640281/150‑1 (MQ=255)
                                                                                                 aaCCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAggg                                                     >  2:318555/1‑150 (MQ=255)
                                                                                                  aCCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAATGGGc                                                    >  2:1702717/1‑150 (MQ=255)
                                                                                                     aTCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTgg                                                 <  2:1503686/150‑1 (MQ=255)
                                                                                                      tCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGGGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTggg                                                <  2:3588088/150‑1 (MQ=255)
                                                                                                        aaaGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGtt                                              >  2:3618761/1‑150 (MQ=255)
                                                                                                        aaaGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGtt                                              >  1:1066760/1‑150 (MQ=255)
                                                                                                          agagCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTAt                                            <  1:2417042/150‑1 (MQ=255)
                                                                                                           gagCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATg                                           <  1:1454749/150‑1 (MQ=255)
                                                                                                             gCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaa                                         <  1:332425/150‑1 (MQ=255)
                                                                                                                 gggATATATTATCTGATACTTATGAAGATATGAAGAGATTACATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAag                                     <  2:1106195/150‑1 (MQ=255)
                                                                                                                 gggATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAag                                     >  1:1141582/1‑150 (MQ=255)
                                                                                                                 gggATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAag                                     <  2:4161285/150‑1 (MQ=255)
                                                                                                                   gATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagag                                   >  2:2376804/1‑150 (MQ=255)
                                                                                                                     tatatTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagat                                 <  2:225829/150‑1 (MQ=255)
                                                                                                                         ttatCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattaga                             <  2:879546/150‑1 (MQ=255)
                                                                                                                          tatCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattagag                            <  1:2830830/150‑1 (MQ=255)
                                                                                                                           atCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattagaga                           >  2:1929614/1‑150 (MQ=255)
                                                                                                                             cTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattagagatt                         >  1:1013307/1‑150 (MQ=255)
                                                                                                                             cTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattagagatt                         <  2:1874794/150‑1 (MQ=255)
                                                                                                                              tGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAgagattagagatta                        <  2:1464985/150‑1 (MQ=255)
                                                                                                                               gATACTTATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTAc                       >  1:210143/1‑150 (MQ=255)
                                                                                                                               gATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGTGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTAc                       <  2:2297934/150‑1 (MQ=255)
                                                                                                                               gATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTAc                       <  2:2091300/150‑1 (MQ=255)
                                                                                                                               gATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTAc                       >  1:709810/1‑150 (MQ=255)
                                                                                                                                    ttATGAAGATATGAAGATATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCTAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGtt                  >  1:308011/1‑150 (MQ=255)
                                                                                                                                    ttATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGtt                  >  1:2510976/1‑150 (MQ=255)
                                                                                                                                    ttATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGtt                  >  1:741907/1‑150 (MQ=255)
                                                                                                                                       tgaagatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTaaa               >  1:101394/1‑150 (MQ=255)
                                                                                                                                        gaagatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAAt              <  2:2347345/150‑1 (MQ=255)
                                                                                                                                        gaagatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAAt              <  1:2907132/150‑1 (MQ=255)
                                                                                                                                            atatgaagaTATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGt          >  1:1770426/1‑150 (MQ=255)
                                                                                                                                             tatgaagaTATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGtt         >  1:3713722/1‑150 (MQ=255)
                                                                                                                                               tgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCg       >  1:3411018/1‑150 (MQ=255)
                                                                                                                                                 aagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTTTGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGt     <  1:3644459/150‑1 (MQ=255)
                                                                                                                                                 aagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGt     >  1:2932760/1‑150 (MQ=255)
                                                                                                                                                 aagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGt     >  2:1004686/1‑150 (MQ=255)
                                                                                                                                                  agatATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTa    <  2:3425734/150‑1 (MQ=255)
                                                                                                                                                    agaTTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATg  >  1:49363/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATATTATCTGATACTTATGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATG  >  CP009273/523632‑523929

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: