Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 181,035 C→A 5.4% P174P (CCG→CCT yaeI ← phosphodiesterase with model substrate bis‑pNPP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273181,0350CA 5.4% 963.4 / 23.0 250P174P (CCG→CCTyaeIphosphodiesterase with model substrate bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (112/124);  new base A (8/6);  total (120/130)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

AAAGGTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAAC  >  CP009273/180886‑181184
                                                                                                                                                     |                                                                                                                                                     
aaaGGTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGc                                                                                                                                                       >  1:3478102/1‑150 (MQ=255)
 aaGGTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCg                                                                                                                                                      >  1:2842191/1‑150 (MQ=255)
   ggTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGt                                                                                                                                                    >  2:2292610/1‑150 (MQ=255)
   ggTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGt                                                                                                                                                    >  1:1847813/1‑150 (MQ=255)
    gTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGtt                                                                                                                                                   >  2:4448389/1‑150 (MQ=255)
     ttCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGttt                                                                                                                                                  >  2:412374/1‑150 (MQ=255)
     ttCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGttt                                                                                                                                                  <  2:3592155/150‑1 (MQ=255)
      tCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTg                                                                                                                                                 >  2:2088848/1‑150 (MQ=255)
          caaccaacGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCAtt                                                                                                                                             <  2:507189/150‑1 (MQ=255)
            accaacGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGt                                                                                                                                           >  1:4081326/1‑150 (MQ=255)
            accaacGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGt                                                                                                                                           <  2:46796/150‑1 (MQ=255)
             ccaacGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTc                                                                                                                                          >  2:1536346/1‑150 (MQ=255)
             ccaacGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGAGGTTTGCATTGTc                                                                                                                                          <  1:1735892/150‑1 (MQ=255)
              caacGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTcc                                                                                                                                         <  2:4392623/150‑1 (MQ=255)
                   gTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGccc                                                                                                                                    >  1:3045270/1‑150 (MQ=255)
                    tACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCa                                                                                                                                   <  2:403923/150‑1 (MQ=255)
                     aCGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCAt                                                                                                                                  <  1:1031283/150‑1 (MQ=255)
                       gcgcAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATaa                                                                                                                                <  2:2484023/150‑1 (MQ=255)
                        cgcAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAg                                                                                                                               >  1:399971/1‑150 (MQ=255)
                        cgcAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAg                                                                                                                               >  1:3531997/1‑150 (MQ=255)
                        cgcAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAg                                                                                                                               >  2:1336223/1‑150 (MQ=255)
                           aGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCa                                                                                                                            >  1:604480/1‑150 (MQ=255)
                                cccGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGt                                                                                                                       >  2:559533/1‑150 (MQ=255)
                                 ccgccgTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTg                                                                                                                      <  2:3465249/150‑1 (MQ=255)
                                   gccgTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGcc                                                                                                                    >  1:337459/1‑150 (MQ=255)
                                   gccgTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGcc                                                                                                                    <  2:3706081/150‑1 (MQ=255)
                                   gccgTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGcc                                                                                                                    <  2:1102434/150‑1 (MQ=255)
                                    ccgTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCt                                                                                                                   <  1:1822930/150‑1 (MQ=255)
                                      gTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTAc                                                                                                                 >  2:1206220/1‑150 (MQ=255)
                                       tGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTAcc                                                                                                                <  2:2261921/150‑1 (MQ=255)
                                        gAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCa                                                                                                               <  2:2517056/150‑1 (MQ=255)
                                        gAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCa                                                                                                               >  1:3968204/1‑150 (MQ=255)
                                           tATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGtt                                                                                                            >  1:3437191/1‑150 (MQ=255)
                                           tATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGtt                                                                                                            <  1:2416806/150‑1 (MQ=255)
                                           tATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGtt                                                                                                            >  1:1420460/1‑150 (MQ=255)
                                           tATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCCGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGtt                                                                                                            <  1:2457681/150‑1 (MQ=255)
                                             tGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGAGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCg                                                                                                          >  2:2467268/1‑150 (MQ=255)
                                              ggCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGa                                                                                                         >  1:1718407/1‑150 (MQ=255)
                                                ccGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAAt                                                                                                       <  1:305308/150‑1 (MQ=255)
                                                 cGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAAtt                                                                                                      >  2:2796775/1‑150 (MQ=255)
                                                  gCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTg                                                                                                     >  2:155755/1‑150 (MQ=255)
                                                  gCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTg                                                                                                     >  2:542024/1‑150 (MQ=255)
                                                   cacaGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGc                                                                                                    >  2:4185126/1‑150 (MQ=255)
                                                    acaGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGcc                                                                                                   <  1:1089981/150‑1 (MQ=255)
                                                     caGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCt                                                                                                  <  1:4504430/150‑1 (MQ=255)
                                                     caGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCt                                                                                                  <  1:1152094/150‑1 (MQ=255)
                                                     caGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCt                                                                                                  <  1:1159669/150‑1 (MQ=255)
                                                     caGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCt                                                                                                  <  2:2748819/150‑1 (MQ=255)
                                                      aGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTg                                                                                                 >  1:1812544/1‑150 (MQ=255)
                                                         aTCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTc                                                                                              <  2:3258497/150‑1 (MQ=255)
                                                           cAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCgg                                                                                            >  2:2664898/1‑150 (MQ=255)
                                                                   aGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGa                                                                                    <  1:252115/150‑1 (MQ=255)
                                                                       ttCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCAc                                                                                >  1:2226294/1‑150 (MQ=255)
                                                                        tcatcaCGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCAcc                                                                               <  1:4334232/150‑1 (MQ=255)
                                                                          atcaCGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGt                                                                             >  1:2887444/1‑150 (MQ=255)
                                                                              atcATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGcc                                                                         >  1:938885/3‑150 (MQ=255)
                                                                               gCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCt                                                                        >  2:1218844/1‑150 (MQ=255)
                                                                               gCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCt                                                                        >  2:2569817/1‑150 (MQ=255)
                                                                               gCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGAGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCt                                                                        <  1:4282512/150‑1 (MQ=255)
                                                                                cATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTg                                                                       >  2:2622461/1‑150 (MQ=255)
                                                                                cATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTg                                                                       >  1:3119815/1‑150 (MQ=255)
                                                                                cATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTg                                                                       >  1:854242/1‑150 (MQ=255)
                                                                                 aTGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTgg                                                                      >  2:3046268/1‑150 (MQ=255)
                                                                                 aTGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTgg                                                                      <  1:2279973/150‑1 (MQ=255)
                                                                                        cttTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAaaca                                                               >  1:4046945/1‑150 (MQ=255)
                                                                                        cttTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAaaca                                                               >  2:3018401/1‑150 (MQ=255)
                                                                                        cttTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAaaca                                                               >  1:3304454/1‑150 (MQ=255)
                                                                                             cTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACAcc                                                          <  2:1297393/150‑1 (MQ=255)
                                                                                              tGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCg                                                         <  1:2251585/150‑1 (MQ=255)
                                                                                               gTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGt                                                        >  2:889003/1‑150 (MQ=255)
                                                                                                 cGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCTAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGa                                                      <  1:3607768/150‑1 (MQ=255)
                                                                                                  gggATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGAt                                                     <  1:79857/150‑1 (MQ=255)
                                                                                                  gggATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGAt                                                     <  2:3789778/150‑1 (MQ=255)
                                                                                                     aTTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGcc                                                  <  2:843802/150‑1 (MQ=255)
                                                                                                     aTTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGcc                                                  <  2:3956883/150‑1 (MQ=255)
                                                                                                            gcCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGAt                                           <  1:264496/150‑1 (MQ=255)
                                                                                                            gcCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGAGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGAt                                           <  1:4162130/150‑1 (MQ=255)
                                                                                                              cAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGAttt                                         <  2:3373587/150‑1 (MQ=255)
                                                                                                              cAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGAttt                                         >  2:1760743/1‑150 (MQ=255)
                                                                                                              cAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGAGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGAttt                                         >  1:164583/1‑150 (MQ=255)
                                                                                                                  aCCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAAc                                     >  1:3877242/1‑150 (MQ=255)
                                                                                                                      gTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGtct                                 <  2:1672702/150‑1 (MQ=255)
                                                                                                                        cTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCg                               >  2:101989/1‑150 (MQ=255)
                                                                                                                        cTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCg                               <  2:1920719/150‑1 (MQ=255)
                                                                                                                          tGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGcc                             <  1:1326683/150‑1 (MQ=255)
                                                                                                                           ggCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCa                            >  2:3346971/1‑150 (MQ=255)
                                                                                                                             cAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAAt                          >  1:4427294/1‑150 (MQ=255)
                                                                                                                             cAGATTGGCTTCGCTGGCGGGAGGAGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAAt                          >  1:1887163/1‑150 (MQ=255)
                                                                                                                              aGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAAtt                         <  2:953421/150‑1 (MQ=255)
                                                                                                                               gATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTa                        <  2:2982736/150‑1 (MQ=255)
                                                                                                                                     cTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGAt                  >  1:4226194/1‑150 (MQ=255)
                                                                                                                                      ttCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGAtt                 <  1:559533/150‑1 (MQ=255)
                                                                                                                                          cTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGAtttttt             <  2:3437191/150‑1 (MQ=255)
                                                                                                                                           tGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGAttttttt            >  2:2824582/1‑150 (MQ=255)
                                                                                                                                           tGGCGGGAGGAGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGAttttttt            >  1:3791151/1‑150 (MQ=255)
                                                                                                                                              cGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTg         <  1:906958/150‑1 (MQ=255)
                                                                                                                                              cGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTg         >  2:1148438/1‑150 (MQ=255)
                                                                                                                                                ggAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCAGTa       >  2:2900313/1‑150 (MQ=255)
                                                                                                                                                  aGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTAcc     <  2:337512/150‑1 (MQ=255)
                                                                                                                                                     cGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAAc  <  1:1536346/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGGTTCACCAACCAACGGTACGCGCAGTTGCCCGCCGTGAGTATGGCCGCACAGCATCAGATCCCAGGGTTCATCACGCATGACTTCTTTGCTGTCGGGATTATGCGCCAGCACCAGTCTTGGCAGATTGGCTTCGCTGGCGGGAGGCGGTTTGCATTGTCCGGCCCATAAGTCACCAGTGCCTACCAGTTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAAC  >  CP009273/180886‑181184

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: