Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,641,686 G→A 5.1% pseudogene (208/279 nt) ydfX ← pseudogene, Qin prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,641,6860GA 5.1% 1132.6 / 31.3 311pseudogene (208/279 nt)ydfXpseudogene, Qin prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (152/143);  new base A (6/10);  total (158/153)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 312 total aligned reads displayed.

TTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATT  >  CP009273/1641537‑1641834
                                                                                                                                                     |                                                                                                                                                    
 tAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGTAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGt                                                                                                                                                     >  1:1399025/1‑150 (MQ=255)
 tAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGt                                                                                                                                                     >  1:3122479/1‑150 (MQ=255)
  aaCATAGTATTTTGAGCGTAACTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGtt                                                                                                                                                    <  2:489529/150‑1 (MQ=255)
  aaCATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGtt                                                                                                                                                    >  1:5082572/1‑150 (MQ=255)
   aCATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTc                                                                                                                                                   >  1:1281312/1‑150 (MQ=255)
    cATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCa                                                                                                                                                  >  2:4172415/1‑150 (MQ=255)
     aTAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCaa                                                                                                                                                 >  1:325786/1‑150 (MQ=255)
      tAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCaaa                                                                                                                                                <  1:410249/150‑1 (MQ=255)
      tAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCaaa                                                                                                                                                >  2:1292674/1‑150 (MQ=255)
        gTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAata                                                                                                                                              <  2:3234847/150‑1 (MQ=255)
           ttttGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCt                                                                                                                                           <  2:4959183/150‑1 (MQ=255)
               gAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCAc                                                                                                                                       <  1:4558284/150‑1 (MQ=255)
                 gCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTc                                                                                                                                     <  2:3901482/150‑1 (MQ=255)
                 gCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTc                                                                                                                                     >  2:5340313/1‑150 (MQ=255)
                 gCGTAAATCAAGACGGGCATGAGCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTc                                                                                                                                     <  2:3077315/150‑1 (MQ=255)
                     aaaTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGAtt                                                                                                                                 >  2:3955657/1‑150 (MQ=255)
                       aTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTaa                                                                                                                               >  1:4868668/1‑150 (MQ=255)
                         cAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACg                                                                                                                             <  2:3446220/150‑1 (MQ=255)
                         cAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACg                                                                                                                             >  2:3554156/1‑150 (MQ=255)
                            tACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCa                                                                                                                          <  2:30028/150‑1 (MQ=255)
                            tACGGGCATTATCCACTCGTTCATGGACAATCTTTTCATCATCCAGTGGTAGCATACATAATTGCAGGCACTAGCGCCGTTAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCa                                                                                                                          <  1:519355/150‑1 (MQ=255)
                             aCGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCaa                                                                                                                         <  2:2069070/150‑1 (MQ=255)
                               gggCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAgc                                                                                                                       <  1:5556418/150‑1 (MQ=255)
                                ggCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAgcg                                                                                                                      >  1:685479/1‑150 (MQ=255)
                                  cATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCCAGcgct                                                                                                                    >  2:4358283/1‑150 (MQ=255)
                                   aTTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctc                                                                                                                   >  1:3962475/1‑150 (MQ=255)
                                     tATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgc                                                                                                                 >  2:1786622/1‑150 (MQ=255)
                                      aTCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgct                                                                                                                <  2:5321195/150‑1 (MQ=255)
                                      aTCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgct                                                                                                                >  1:5489320/1‑150 (MQ=255)
                                          aCTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATc                                                                                                            >  1:2484902/1‑150 (MQ=255)
                                             cGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAAt                                                                                                         <  1:1656398/150‑1 (MQ=255)
                                             cGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAAt                                                                                                         <  2:5475956/150‑1 (MQ=255)
                                               ttCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcg                                                                                                       <  1:3029970/150‑1 (MQ=255)
                                                tCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcgc                                                                                                      >  1:1054847/1‑150 (MQ=255)
                                                   tGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGcc                                                                                                   >  1:2229918/1‑150 (MQ=255)
                                                    ggACACTCTTATCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCa                                                                                                  <  2:2909386/150‑1 (MQ=255)
                                                       cacTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCt                                                                                               >  2:4144492/1‑150 (MQ=255)
                                                        acTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCtt                                                                                              >  2:1804747/1‑150 (MQ=255)
                                                           ctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCat                                                                                           >  2:1039964/1‑150 (MQ=255)
                                                             tttCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCatat                                                                                         <  1:3848387/150‑1 (MQ=255)
                                                              ttCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATAtg                                                                                        >  1:852835/1‑150 (MQ=255)
                                                               tcatcaTCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTACATAtgt                                                                                       >  1:5023242/1‑150 (MQ=255)
                                                                catcatCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATAtgtg                                                                                      >  1:994817/1‑150 (MQ=255)
                                                                    atcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGacac                                                                                  <  1:4194193/150‑1 (MQ=255)
                                                                     tcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACg                                                                                 >  1:1963819/1‑150 (MQ=255)
                                                                           ggTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAAc                                                                           >  1:3555661/1‑150 (MQ=255)
                                                                             tAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACAt                                                                         >  2:5674306/1‑150 (MQ=255)
                                                                              aGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATc                                                                        >  2:1412427/1‑150 (MQ=255)
                                                                               gCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCg                                                                       >  1:5576274/1‑150 (MQ=255)
                                                                                ccTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGa                                                                      >  2:3270641/1‑150 (MQ=255)
                                                                                 cTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGAt                                                                     <  1:2439393/150‑1 (MQ=255)
                                                                                  tCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATg                                                                    >  2:810973/1‑150 (MQ=255)
                                                                                     aTAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTc                                                                 <  2:3354785/150‑1 (MQ=255)
                                                                                     aTAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCAGATGTGACACGCTGAACATCGATGTTc                                                                 >  2:5571250/1‑150 (MQ=255)
                                                                                        aTTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCggg                                                              >  2:4851350/1‑150 (MQ=255)
                                                                                        aTTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCggg                                                              <  2:3122479/150‑1 (MQ=255)
                                                                                         ttGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGgc                                                             <  2:5616529/150‑1 (MQ=255)
                                                                                          tGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGgcg                                                            <  2:425570/150‑1 (MQ=255)
                                                                                          tGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGgcg                                                            >  1:2863605/1‑150 (MQ=255)
                                                                                               gCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCtt                                                       <  2:714238/150‑1 (MQ=255)
                                                                                                cACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTTTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCttt                                                      <  2:2277553/150‑1 (MQ=255)
                                                                                                cACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCttt                                                      >  2:2374263/1‑150 (MQ=255)
                                                                                                 aCTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTc                                                     >  2:2582796/1‑150 (MQ=255)
                                                                                                 aCTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTc                                                     <  2:1719887/150‑1 (MQ=255)
                                                                                                    aGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCaaa                                                  >  2:2447047/1‑150 (MQ=255)
                                                                                                    aGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCaaa                                                  <  2:29155/150‑1 (MQ=255)
                                                                                                        ccGTGAACTTTTCACGCTTATCACTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGccc                                              <  1:4760532/150‑1 (MQ=255)
                                                                                                         cGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCa                                             >  2:1134905/1‑150 (MQ=255)
                                                                                                          gTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAg                                            >  1:3341914/1‑150 (MQ=255)
                                                                                                          gTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAg                                            >  1:4916625/1‑150 (MQ=255)
                                                                                                           tGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGa                                           <  1:3942948/150‑1 (MQ=255)
                                                                                                             aaCTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAAt                                         <  1:4650431/150‑1 (MQ=255)
                                                                                                                   tCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCa                                   >  2:526262/1‑150 (MQ=255)
                                                                                                                    cACGCTTATCCCTGGTGTAGATAGCATTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAg                                  <  2:4512509/150‑1 (MQ=255)
                                                                                                                     aCGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGt                                 >  1:3741065/1‑150 (MQ=255)
                                                                                                                         ttATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTAGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGAtt                             <  1:1360979/150‑1 (MQ=255)
                                                                                                                         ttATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGAtt                             >  2:4333974/1‑150 (MQ=255)
                                                                                                                          tATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTa                            <  1:2790149/150‑1 (MQ=255)
                                                                                                                          tATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTa                            >  1:4467069/1‑150 (MQ=255)
                                                                                                                               cTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtc                       >  1:1822024/1‑150 (MQ=255)
                                                                                                                               cTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtc                       <  1:5381033/150‑1 (MQ=255)
                                                                                                                                tGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcg                      >  2:5036345/1‑150 (MQ=255)
                                                                                                                                tGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcg                      <  1:3972014/150‑1 (MQ=255)
                                                                                                                                 ggTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcgt                     <  2:516717/150‑1 (MQ=255)
                                                                                                                                 ggTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcgt                     >  2:4087135/1‑150 (MQ=255)
                                                                                                                                   tgtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTGCAGTGATTAATAtcgtcg                   >  2:3295631/1‑150 (MQ=255)
                                                                                                                                   tgtCGAAAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTAGCTGATCAATCGCGCCCCCTTCATATTTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGAGTAATAtcgtcg                   >  2:3804946/1‑150 (MQ=255)
                                                                                                                                    gtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                  <  2:3670622/150‑1 (MQ=255)
                                                                                                                                      cGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCTATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCTAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCt                <  2:685479/150‑1 (MQ=255)
                                                                                                                                       gATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTg               >  2:2168011/1‑150 (MQ=255)
                                                                                                                                         tAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCt             >  1:3621056/1‑150 (MQ=255)
                                                                                                                                          aGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTc            <  1:5607680/150‑1 (MQ=255)
                                                                                                                                             cTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCtg         >  1:4826362/1‑150 (MQ=255)
                                                                                                                                               tCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCGTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCtgtg       >  2:1835763/1‑150 (MQ=255)
                                                                                                                                               tCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCGTCATATGAGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAATGCCCAGAATGCTTCAGTGATTAATATCTTCGCCTGCTCCtgtg       >  2:4624424/1‑150 (MQ=255)
                                                                                                                                                 cAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTc     <  2:2790000/150‑1 (MQ=255)
                                                                                                                                                  aGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCa    <  2:882449/150‑1 (MQ=255)
                                                                                                                                                    cGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCAtt  >  2:1634315/1‑150 (MQ=255)
                                                                                                                                                    cGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCAtt  <  2:4453127/150‑1 (MQ=255)
                                                                                                                                                    cGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCAtt  >  1:2005666/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATT  >  CP009273/1641537‑1641834

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: