Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 524,983 G→T 12.5% pseudogene (28/93 nt) ylbI → pseudogene, Rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,9830GT12.5% 744.0 / 57.4 232pseudogene (28/93 nt)ylbIpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (111/92);  new base T (19/10);  total (130/102)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 235 total aligned reads displayed.

ACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAGGTTAT  >  CP009273/524834‑525132
                                                                                                                                                     |                                                                                                                                                     
 cATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGa                                                                                                                                                      >  2:4779912/1‑150 (MQ=255)
    atTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACAt                                                                                                                                                   >  2:3577912/1‑150 (MQ=255)
     ttCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATg                                                                                                                                                  <  1:822678/150‑1 (MQ=255)
      tCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGc                                                                                                                                                 >  2:1607045/1‑150 (MQ=255)
         gCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGcccc                                                                                                                                              >  1:1783295/1‑150 (MQ=255)
          ccATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCa                                                                                                                                             >  1:3432676/1‑150 (MQ=255)
          ccATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCa                                                                                                                                             >  2:3819833/1‑150 (MQ=255)
           cATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAg                                                                                                                                            >  2:2310260/1‑150 (MQ=255)
           cATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTATTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAg                                                                                                                                            >  1:538901/1‑150 (MQ=255)
             tGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTa                                                                                                                                          >  2:1451604/1‑150 (MQ=255)
                  taattaaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTAcggcg                                                                                                                                     >  1:1749106/1‑150 (MQ=255)
                  taattaaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTAcggcg                                                                                                                                     >  1:667951/1‑150 (MQ=255)
                  taattaaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTAcggcg                                                                                                                                     >  1:4328166/1‑150 (MQ=255)
                    attaaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACggcggc                                                                                                                                   >  2:66426/1‑150 (MQ=255)
                    attaaAATGTTACGTGTATAATGTGAGTCCCTATCTGAAAATATATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTCGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACggcggc                                                                                                                                   <  2:5072661/150‑1 (MQ=255)
                      taaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGccc                                                                                                                                 <  2:4183977/150‑1 (MQ=255)
                      taaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGccc                                                                                                                                 >  1:373982/1‑150 (MQ=255)
                       aaaaTGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCg                                                                                                                                <  2:5096539/150‑1 (MQ=255)
                         aaTGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGAt                                                                                                                              <  1:4841706/150‑1 (MQ=255)
                         aaTGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGAt                                                                                                                              <  2:3710186/150‑1 (MQ=255)
                          aTGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATa                                                                                                                             <  2:399826/150‑1 (MQ=255)
                             ttACGTGTTTAATGTGAGTCCCGATCTGAAAATAAATAAGCCTTCCGGATTAAAAGAACTTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTc                                                                                                                          <  2:4222070/150‑1 (MQ=255)
                              tACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTcc                                                                                                                         >  1:1609282/1‑150 (MQ=255)
                              tACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTcc                                                                                                                         <  1:1611759/150‑1 (MQ=255)
                               aCGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCa                                                                                                                        <  2:1814669/150‑1 (MQ=255)
                                cGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAg                                                                                                                       >  2:4325474/1‑150 (MQ=255)
                                cGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAg                                                                                                                       <  2:210292/150‑1 (MQ=255)
                                 gtgtTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAGTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAgg                                                                                                                      <  2:5096170/150‑1 (MQ=255)
                                  tgtTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAggg                                                                                                                     >  2:1571964/1‑150 (MQ=255)
                                   gtTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGt                                                                                                                    >  1:5102681/1‑150 (MQ=255)
                                    tttAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGtt                                                                                                                   >  2:2678731/1‑150 (MQ=255)
                                      tAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCa                                                                                                                 <  1:363215/150‑1 (MQ=255)
                                      tAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCa                                                                                                                 <  1:363257/150‑1 (MQ=255)
                                       aaTGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAg                                                                                                                >  2:2823648/1‑150 (MQ=255)
                                         tgtgAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCg                                                                                                              >  2:3254062/1‑150 (MQ=255)
                                                cccTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTAc                                                                                                       >  1:1011814/1‑150 (MQ=255)
                                                   tATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGa                                                                                                    <  1:389263/150‑1 (MQ=255)
                                                    aTCTGAAAATAAATAATCCTTCCGGATTAAAATAAAGTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGAt                                                                                                   <  2:832993/150‑1 (MQ=255)
                                                         aaaaTAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGTGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcg                                                                                              <  2:4944044/150‑1 (MQ=255)
                                                          aaataaatAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgc                                                                                             <  1:3170306/150‑1 (MQ=255)
                                                             taaataaTCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGcc                                                                                          <  1:4916329/150‑1 (MQ=255)
                                                                 taatCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACt                                                                                      >  1:4254283/1‑150 (MQ=255)
                                                                   atCCTTCCGGATTAAAATAAATTCGTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGc                                                                                    <  1:2117786/150‑1 (MQ=255)
                                                                    tCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAGTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCa                                                                                   <  1:2441260/150‑1 (MQ=255)
                                                                      cTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATc                                                                                 <  1:233320/150‑1 (MQ=255)
                                                                       ttCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCg                                                                                >  1:1887906/1‑150 (MQ=255)
                                                                                  aataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATa                                                                     >  1:1326743/1‑150 (MQ=255)
                                                                                  aataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATa                                                                     >  1:2301755/1‑150 (MQ=255)
                                                                                     aaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACaa                                                                  <  1:1601622/150‑1 (MQ=255)
                                                                                      aatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAAc                                                                 >  2:1995343/1‑150 (MQ=255)
                                                                                       atTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACt                                                                >  2:4942071/1‑150 (MQ=255)
                                                                                       atTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACt                                                                >  1:3994091/1‑150 (MQ=255)
                                                                                       atTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACt                                                                >  2:1393717/1‑150 (MQ=255)
                                                                                       atTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACt                                                                <  1:305748/150‑1 (MQ=255)
                                                                                         tCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttt                                                              <  1:2088515/150‑1 (MQ=255)
                                                                                         tCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttt                                                              >  1:2950953/1‑150 (MQ=255)
                                                                                          cTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCTAGAACTCCATACAACtttt                                                             <  2:4893715/150‑1 (MQ=255)
                                                                                          cTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACtttt                                                             >  2:923283/1‑150 (MQ=255)
                                                                                           ttGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttttt                                                            >  1:4490353/1‑150 (MQ=255)
                                                                                            tGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTa                                                           >  1:4431493/1‑150 (MQ=255)
                                                                                            tGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTa                                                           >  2:3767257/1‑150 (MQ=255)
                                                                                            tGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTa                                                           >  1:4086992/1‑150 (MQ=255)
                                                                                              ccGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACg                                                         >  2:747968/1‑150 (MQ=255)
                                                                                               cGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACgg                                                        <  1:4298637/150‑1 (MQ=255)
                                                                                               cGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACgg                                                        >  1:2261138/1‑150 (MQ=255)
                                                                                                gggAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGa                                                       >  1:1071084/1‑150 (MQ=255)
                                                                                                 ggAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACggta                                                      <  1:1451604/150‑3 (MQ=255)
                                                                                                  gaaagaaaGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAAc                                                     <  2:3123703/150‑1 (MQ=255)
                                                                                                   aaagaaagAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACg                                                    >  2:202793/1‑150 (MQ=255)
                                                                                                     agaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTc                                                  <  2:3432966/150‑1 (MQ=255)
                                                                                                      gaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCa                                                 >  1:2195517/1‑150 (MQ=255)
                                                                                                      gaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCa                                                 <  2:1889259/150‑1 (MQ=255)
                                                                                                         agagGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCt                                              >  2:247324/1‑150 (MQ=255)
                                                                                                         agagGAAATAAACCATGAGCGGAAAACCAGCGGCACGCCAGGGTGACTTGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGAGGGGGTACTGATAGTCGCGCCGACTGCATCGAGATCTCCATACAACTTTTTACGGAACGTCACCt                                              <  1:646958/150‑1 (MQ=255)
                                                                                                             gAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATa                                          <  1:4732436/150‑1 (MQ=255)
                                                                                                              aaaTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATaa                                         <  2:2408631/150‑1 (MQ=255)
                                                                                                                aTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATaaaa                                       <  1:733241/150‑1 (MQ=255)
                                                                                                                 tAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAAt                                      >  1:2401748/1‑150 (MQ=255)
                                                                                                                 tAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAAt                                      <  1:980625/150‑1 (MQ=255)
                                                                                                                 tAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAAt                                      <  1:980648/150‑1 (MQ=255)
                                                                                                                  aaaCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGAGAAAAtt                                     >  1:4322982/1‑150 (MQ=255)
                                                                                                                   aaCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAAttt                                    >  2:4116934/1‑150 (MQ=255)
                                                                                                                    aCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTg                                   <  1:2808505/150‑1 (MQ=255)
                                                                                                                    aCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTg                                   <  1:2697747/150‑1 (MQ=255)
                                                                                                                     ccATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGt                                  >  2:703810/1‑150 (MQ=255)
                                                                                                                        ttAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCa                               >  1:1786975/1‑150 (MQ=255)
                                                                                                                         tAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAg                              >  2:1239764/1‑150 (MQ=255)
                                                                                                                         tAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAg                              <  2:3627594/150‑1 (MQ=255)
                                                                                                                            cGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTAt                           >  2:4154992/1‑150 (MQ=255)
                                                                                                                            cGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTAt                           >  2:2744421/1‑150 (MQ=255)
                                                                                                                              gAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCg                         >  2:201885/1‑150 (MQ=255)
                                                                                                                               aaaaCCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGt                        >  1:3433705/1‑150 (MQ=255)
                                                                                                                                aaaCCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTc                       <  1:4779912/150‑1 (MQ=255)
                                                                                                                                  aCCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCtt                     <  1:3511047/150‑1 (MQ=255)
                                                                                                                                     aGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTg                  <  2:3512635/150‑1 (MQ=255)
                                                                                                                                      gCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTga                 <  1:5016502/150‑1 (MQ=255)
                                                                                                                                       cGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTgag                >  1:2627300/1‑150 (MQ=255)
                                                                                                                                        ggCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTgaga               >  2:3278014/1‑150 (MQ=255)
                                                                                                                                            cGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtct           >  1:1070355/1‑150 (MQ=255)
                                                                                                                                            cGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtct           >  1:4205714/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAGGTTAT  >  CP009273/524834‑525132

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: