Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,798,045 G→T 5.6% pseudogene (271/1476 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,798,0450GT 5.6% 1218.4 / 27.0 304pseudogene (271/1476 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base G (131/155);  new base T (12/6);  total (143/161)
Fisher's exact test for biased strand distribution p-value = 9.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

CATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTG  >  CP009273/1797896‑1798194
                                                                                                                                                     |                                                                                                                                                     
cATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAgg                                                                                                                                                       <  1:1736721/150‑1 (MQ=255)
   tGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtt                                                                                                                                                    >  1:4095245/1‑150 (MQ=255)
   tGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtt                                                                                                                                                    >  1:1335312/1‑150 (MQ=255)
   tGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtt                                                                                                                                                    >  2:1233057/1‑150 (MQ=255)
   tGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATACATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTTGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtt                                                                                                                                                    >  2:1328734/1‑150 (MQ=255)
     agaAATGGAAAAGGCAAAGTCCACCTGATGAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtttt                                                                                                                                                  <  2:3121796/150‑1 (MQ=255)
     agaAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtttt                                                                                                                                                  <  2:4832688/150‑1 (MQ=255)
     agaAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtttt                                                                                                                                                  <  2:2573606/150‑1 (MQ=255)
       aaaTGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTgg                                                                                                                                                <  2:1461003/150‑1 (MQ=255)
         aTGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCt                                                                                                                                              >  2:1174898/1‑150 (MQ=255)
         aTGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCt                                                                                                                                              >  1:2000209/1‑150 (MQ=255)
          tGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                                                                                                             >  2:2490587/1‑150 (MQ=255)
           ggAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTc                                                                                                                                            <  2:4687316/150‑1 (MQ=255)
           ggAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTc                                                                                                                                            <  2:2137869/150‑1 (MQ=255)
            gAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCg                                                                                                                                           <  1:2401045/150‑1 (MQ=255)
             aaaaGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGc                                                                                                                                          >  1:4361322/1‑150 (MQ=255)
                aGGCACAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCttt                                                                                                                                       <  2:3489790/150‑1 (MQ=255)
                  gCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAg                                                                                                                                     <  2:2054941/150‑1 (MQ=255)
                   cAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGa                                                                                                                                    <  1:1741896/150‑1 (MQ=255)
                      aGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGatat                                                                                                                                 >  1:3829321/1‑150 (MQ=255)
                      aGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGatat                                                                                                                                 >  1:4893486/1‑150 (MQ=255)
                              tGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCtttt                                                                                                                         <  1:2490587/150‑1 (MQ=255)
                               gATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTg                                                                                                                        >  2:4922841/1‑150 (MQ=255)
                               gATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTg                                                                                                                        <  2:2848492/150‑1 (MQ=255)
                               gATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTg                                                                                                                        <  1:3443589/150‑1 (MQ=255)
                                aTTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGt                                                                                                                       <  2:1053510/150‑1 (MQ=255)
                                  tAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTGCGCTTTAGATATTTTCTTTTGTCg                                                                                                                     >  2:1392702/1‑150 (MQ=255)
                                    aaCATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGaa                                                                                                                   >  1:1099417/1‑150 (MQ=255)
                                     aCATTAATCATCCGTTGAATCCATACAGCAGCAAGTCACATTAGTCCGGGACTCTGAGTTCAGATCATTACTTCGCCATCATCAAGCAGGGAGTCGCGTTATAGCGAGTAATGATTTTGGCTTCGATTTAGATATATTCATTTGTCTAat                                                                                                                  >  1:2536464/1‑150 (MQ=255)
                                      cATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAatg                                                                                                                 <  1:2587040/150‑1 (MQ=255)
                                            tcctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATTTTTTCTTTTGTCGAATGATGGCa                                                                                                           <  1:1709042/150‑1 (MQ=255)
                                            tcctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCa                                                                                                           >  1:1434894/1‑150 (MQ=255)
                                            tcctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCa                                                                                                           <  1:4818120/150‑1 (MQ=255)
                                            tcctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCa                                                                                                           >  2:3537910/1‑150 (MQ=255)
                                             cctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaa                                                                                                          >  2:19744/1‑150 (MQ=255)
                                             cctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaa                                                                                                          <  1:4456777/150‑1 (MQ=255)
                                              ctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaaa                                                                                                         <  2:4038769/150‑1 (MQ=255)
                                                 cGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACaa                                                                                                      <  2:320500/150‑1 (MQ=255)
                                                   ttGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATg                                                                                                    >  1:1670632/1‑150 (MQ=255)
                                                   ttGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATg                                                                                                    <  2:384418/150‑1 (MQ=255)
                                                    tGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgt                                                                                                   <  2:2382456/150‑1 (MQ=255)
                                                      aaTCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtca                                                                                                 <  2:2998217/150‑1 (MQ=255)
                                                      aaTCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtca                                                                                                 >  1:71296/1‑150 (MQ=255)
                                                      aaTCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtca                                                                                                 >  2:114704/1‑150 (MQ=255)
                                                       aTCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtcag                                                                                                <  1:549045/150‑1 (MQ=255)
                                                             aCAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCatat                                                                                          >  1:628065/1‑150 (MQ=255)
                                                              cagcagCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATg                                                                                         >  2:586046/1‑150 (MQ=255)
                                                               agcagcATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGt                                                                                        >  2:4467074/1‑150 (MQ=255)
                                                                gcagcaTGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTc                                                                                       <  2:3343294/150‑1 (MQ=255)
                                                                  agcaTGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCcg                                                                                     <  1:1189425/150‑1 (MQ=255)
                                                                   gcaTGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCcgc                                                                                    >  1:3933012/1‑150 (MQ=255)
                                                                       gTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGAtt                                                                                >  2:1687715/1‑150 (MQ=255)
                                                                        tCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGAttt                                                                               <  1:4039506/150‑1 (MQ=255)
                                                                            cTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtata                                                                           >  2:4051098/1‑150 (MQ=255)
                                                                             ttAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatat                                                                          >  1:2197324/1‑150 (MQ=255)
                                                                               aGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatatat                                                                        <  1:2803983/150‑1 (MQ=255)
                                                                                 tCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAg                                                                      <  2:2664613/150‑1 (MQ=255)
                                                                                   cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                                                                    <  2:4775497/150‑1 (MQ=255)
                                                                                   cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                                                                    <  2:4933516/150‑1 (MQ=255)
                                                                                   cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                                                                    >  1:265576/1‑150 (MQ=255)
                                                                                    gtgACGCTGAGTTCAGAGCATTACTTCGACATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATTGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACa                                                                   <  1:4653622/148‑1 (MQ=255)
                                                                                    gggACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACa                                                                   >  1:2354155/1‑150 (MQ=255)
                                                                                       aCGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAAATATGTCAGTCATATGTCCGCGATTTTATATATAGACAATc                                                                <  2:1531894/150‑1 (MQ=255)
                                                                                        cgcgGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCa                                                               >  1:1890586/1‑150 (MQ=255)
                                                                                          cgGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATa                                                             >  2:370069/1‑150 (MQ=255)
                                                                                              gTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCTAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACttt                                                         <  1:183708/150‑1 (MQ=255)
                                                                                              gTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACttt                                                         <  1:5027064/150‑1 (MQ=255)
                                                                                               ttCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTg                                                        >  1:2053516/1‑150 (MQ=255)
                                                                                               ttCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTg                                                        >  1:140006/1‑150 (MQ=255)
                                                                                                 cAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTc                                                      <  2:3865432/150‑1 (MQ=255)
                                                                                                   gagCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAg                                                    <  2:685921/150‑1 (MQ=255)
                                                                                                   gagCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAg                                                    <  1:3754784/150‑1 (MQ=255)
                                                                                                        ttACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAgg                                               <  1:2951558/150‑1 (MQ=255)
                                                                                                          aCTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCt                                             <  2:1790980/150‑1 (MQ=255)
                                                                                                            ttCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTAc                                           >  2:390282/1‑150 (MQ=255)
                                                                                                             tCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACg                                          >  2:2600520/1‑150 (MQ=255)
                                                                                                               gCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGaa                                        <  1:1858403/150‑1 (MQ=255)
                                                                                                                ccATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAAt                                       >  2:41600/1‑150 (MQ=255)
                                                                                                                 cATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATa                                      >  1:4505675/1‑150 (MQ=255)
                                                                                                                  aTCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaa                                     <  2:2000209/150‑1 (MQ=255)
                                                                                                                   tcctccTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaaa                                    <  2:3324099/150‑1 (MQ=255)
                                                                                                                   tcctccTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaaa                                    <  2:1886229/150‑1 (MQ=255)
                                                                                                                   tcctccTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaaa                                    >  2:670283/1‑150 (MQ=255)
                                                                                                                         tGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGcagtc                              <  1:1328636/150‑1 (MQ=255)
                                                                                                                          gCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGcagtca                             <  1:3349454/150‑1 (MQ=255)
                                                                                                                             gggAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAgtg                          <  1:635169/150‑1 (MQ=255)
                                                                                                                                aGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTAt                       <  1:1160351/150‑1 (MQ=255)
                                                                                                                                aGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTAt                       <  2:2710000/150‑1 (MQ=255)
                                                                                                                                   cgcgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGa                    >  1:4211236/1‑150 (MQ=255)
                                                                                                                                   cgcgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGa                    <  2:881165/150‑1 (MQ=255)
                                                                                                                                    gcgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGAc                   <  1:437131/150‑1 (MQ=255)
                                                                                                                                     cgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATATACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACa                  >  2:3194182/1‑150 (MQ=255)
                                                                                                                                      gTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACaa                 >  1:1475522/1‑150 (MQ=255)
                                                                                                                                             cgAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACaattaaatt          >  1:3525404/1‑150 (MQ=255)
                                                                                                                                               aGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAt        >  1:3034484/1‑150 (MQ=255)
                                                                                                                                                gTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAtg       <  1:4562098/150‑1 (MQ=255)
                                                                                                                                                   aGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAtgttg    >  2:33968/1‑150 (MQ=255)
                                                                                                                                                   aGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAtgttg    >  2:5079064/1‑150 (MQ=255)
                                                                                                                                                    ggATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAtgttgt   >  2:2283805/1‑150 (MQ=255)
                                                                                                                                                     gATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTtgtg  >  1:4846754/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTG  >  CP009273/1797896‑1798194

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: