Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,452,752 T→A 12.9% V352E (GTG→GAG)  gspE → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,452,7520TA12.9% 1018.7 / 120.6 326V352E (GTG→GAG) gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (140/144);  new base A (19/23);  total (159/167)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 331 total aligned reads displayed.

GGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATC  >  CP009273/3452604‑3452901
                                                                                                                                                    |                                                                                                                                                     
 gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGa                                                                                                                                                     >  2:1423565/1‑150 (MQ=255)
   cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGAcc                                                                                                                                                   >  1:3714728/1‑150 (MQ=255)
   cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGAcc                                                                                                                                                   >  2:4339753/1‑150 (MQ=255)
     gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCg                                                                                                                                                 <  1:1016289/150‑1 (MQ=255)
     gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCg                                                                                                                                                 <  1:735029/150‑1 (MQ=255)
        cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCt                                                                                                                                              >  1:1438435/1‑150 (MQ=255)
         gcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTc                                                                                                                                             >  2:1979357/1‑150 (MQ=255)
         gcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTc                                                                                                                                             >  2:4291125/1‑150 (MQ=255)
               cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGAc                                                                                                                                       <  1:534430/150‑1 (MQ=255)
               cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGAc                                                                                                                                       <  1:1659132/150‑1 (MQ=255)
                ttCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACa                                                                                                                                      <  2:1169582/150‑1 (MQ=255)
                 tCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACAt                                                                                                                                     <  2:3032787/150‑1 (MQ=255)
                  cGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATg                                                                                                                                    >  2:1104080/1‑150 (MQ=255)
                  cGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATg                                                                                                                                    <  1:3630531/150‑1 (MQ=255)
                    ccAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATggg                                                                                                                                  <  1:59469/150‑1 (MQ=255)
                    ccAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATggg                                                                                                                                  <  2:363713/150‑1 (MQ=255)
                    ccAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATggg                                                                                                                                  >  2:3483601/1‑150 (MQ=255)
                      aGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCg                                                                                                                                >  2:1022658/1‑150 (MQ=255)
                       ggACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGt                                                                                                                               <  1:3900884/150‑1 (MQ=255)
                        gACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTc                                                                                                                              <  2:3865643/150‑1 (MQ=255)
                          cccGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGa                                                                                                                            <  2:1657415/150‑1 (MQ=255)
                            cGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAAt                                                                                                                          <  1:1176583/150‑1 (MQ=255)
                             ggATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTCACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATc                                                                                                                         >  1:3336594/1‑150 (MQ=255)
                             ggATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATc                                                                                                                         <  2:651064/150‑1 (MQ=255)
                             ggATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATc                                                                                                                         >  1:502243/1‑150 (MQ=255)
                             ggATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATc                                                                                                                         <  2:3078892/150‑1 (MQ=255)
                                 gtcgtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTc                                                                                                                     >  1:3160036/1‑150 (MQ=255)
                                   cgtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCt                                                                                                                   <  1:2517961/150‑1 (MQ=255)
                                    gtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTg                                                                                                                  >  1:1845739/1‑150 (MQ=255)
                                    gtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTg                                                                                                                  >  2:2424687/1‑150 (MQ=255)
                                      cATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCt                                                                                                                <  1:4430667/150‑1 (MQ=255)
                                      cATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCt                                                                                                                >  1:609477/1‑150 (MQ=255)
                                        tggtggGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCttt                                                                                                              <  2:1707176/150‑1 (MQ=255)
                                        tggtggGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCttt                                                                                                              >  1:2952314/1‑150 (MQ=255)
                                                gAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCCTCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTcc                                                                                                      >  2:1185220/1‑150 (MQ=255)
                                                  aaTTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCt                                                                                                    <  1:791729/150‑1 (MQ=255)
                                                  aaTTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCt                                                                                                    >  1:1378379/1‑150 (MQ=255)
                                                   aTTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTg                                                                                                   <  2:3560490/150‑1 (MQ=255)
                                                   aTTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTg                                                                                                   <  2:553162/150‑1 (MQ=255)
                                                   aTTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTg                                                                                                   >  1:1144503/1‑150 (MQ=255)
                                                         gATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAggg                                                                                             >  1:2814717/1‑150 (MQ=255)
                                                           tACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGAt                                                                                           <  2:740893/150‑1 (MQ=255)
                                                           tACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGAt                                                                                           <  1:3615978/150‑1 (MQ=255)
                                                            aCAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGAtt                                                                                          <  1:2164412/150‑1 (MQ=255)
                                                              aGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTAt                                                                                        >  1:2626121/1‑150 (MQ=255)
                                                                aaaCCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcg                                                                                      <  1:3040838/150‑1 (MQ=255)
                                                                aaaCCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcg                                                                                      >  2:4024605/1‑150 (MQ=255)
                                                                 gaCCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgc                                                                                     >  1:2362039/2‑150 (MQ=255)
                                                                 aaCCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgc                                                                                     >  1:3369993/1‑150 (MQ=255)
                                                                   ccGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCgcgc                                                                                   >  1:131271/1‑150 (MQ=255)
                                                                      cgcAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAAc                                                                                <  1:3679256/150‑1 (MQ=255)
                                                                      cgcAGATTGCGGTTCAGGCCTCGCTCACCGGACATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAAc                                                                                <  1:3063262/150‑1 (MQ=255)
                                                                       gcAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACg                                                                               <  1:954317/150‑1 (MQ=255)
                                                                          gATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCt                                                                            <  1:4495315/150‑1 (MQ=255)
                                                                          gATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCt                                                                            <  1:2579536/150‑1 (MQ=255)
                                                                              gCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGtt                                                                        >  1:3837189/1‑150 (MQ=255)
                                                                                ggTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcg                                                                      <  1:789375/150‑1 (MQ=255)
                                                                                ggTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcg                                                                      >  2:2806995/1‑150 (MQ=255)
                                                                                ggTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcg                                                                      <  2:333617/150‑1 (MQ=255)
                                                                                  ttCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcgtc                                                                    <  2:1810981/150‑1 (MQ=255)
                                                                                   tCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcgtcg                                                                   <  2:2067228/150‑1 (MQ=255)
                                                                                     aGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGcc                                                                 >  1:1133381/1‑150 (MQ=255)
                                                                                      ggCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCt                                                                >  1:4254959/1‑150 (MQ=255)
                                                                                       gCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCtg                                                               >  2:3172673/1‑150 (MQ=255)
                                                                                        ccTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCtgt                                                              <  2:321441/150‑1 (MQ=255)
                                                                                           cGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTc                                                           <  2:1584842/150‑1 (MQ=255)
                                                                                                  cGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATg                                                    >  2:1489227/1‑150 (MQ=255)
                                                                                                  cGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATg                                                    >  2:3588846/1‑150 (MQ=255)
                                                                                                    gCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGcc                                                  <  1:4291125/150‑1 (MQ=255)
                                                                                                      cATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGa                                                >  1:1786334/1‑150 (MQ=255)
                                                                                                      cATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGa                                                >  1:1816307/1‑150 (MQ=255)
                                                                                                       aTCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGAc                                               <  2:617764/150‑1 (MQ=255)
                                                                                                        tCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACa                                              <  1:2886848/150‑1 (MQ=255)
                                                                                                         cTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACaa                                             >  2:270032/1‑150 (MQ=255)
                                                                                                            gTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTc                                          <  1:1739319/150‑1 (MQ=255)
                                                                                                            gTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGAGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTc                                          >  1:4470799/1‑150 (MQ=255)
                                                                                                              aCTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCAc                                        >  2:1906868/1‑150 (MQ=255)
                                                                                                              aCTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCAc                                        >  2:3838268/1‑150 (MQ=255)
                                                                                                               ctctCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACg                                       <  2:3274102/150‑1 (MQ=255)
                                                                                                                  tcAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGccc                                    >  1:3932824/1‑150 (MQ=255)
                                                                                                                  tcAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGccc                                    >  1:2528001/1‑150 (MQ=255)
                                                                                                                    aaCACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGt                                  >  1:2441899/1‑150 (MQ=255)
                                                                                                                     acacTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACCATTCACGCCCGTa                                 >  2:912942/1‑150 (MQ=255)
                                                                                                                         tCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTAGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAc                             >  2:707969/1‑150 (MQ=255)
                                                                                                                          ccACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAcc                            <  1:2763932/150‑1 (MQ=255)
                                                                                                                          ccACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAcc                            <  1:2424687/150‑1 (MQ=255)
                                                                                                                               cTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCcaac                       <  1:2012011/150‑1 (MQ=255)
                                                                                                                                   cAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAgc                   <  1:4111936/150‑1 (MQ=255)
                                                                                                                                   cAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAgc                   >  2:2695690/1‑150 (MQ=255)
                                                                                                                                    aGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACACTTCACGCCCGTATCACCCCAACAAgcg                  <  2:1273486/150‑1 (MQ=255)
                                                                                                                                     gTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAgcgc                 >  2:232107/1‑150 (MQ=255)
                                                                                                                                        cATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGa              <  2:2899947/150‑1 (MQ=255)
                                                                                                                                        cATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGa              <  2:1844581/150‑1 (MQ=255)
                                                                                                                                         aTCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGAt             <  2:2569733/150‑1 (MQ=255)
                                                                                                                                              gcgcAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTa        <  2:3145626/150‑1 (MQ=255)
                                                                                                                                               cgcAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTaa       >  2:4090543/1‑150 (MQ=255)
                                                                                                                                                 cAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGTTTTCGTCTTCCCTGGCAGGGATTATCGCGCAAGGTCTGGTTCGTCGCCTGTGTCCGAAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAat     <  1:1162817/150‑1 (MQ=255)
                                                                                                                                                 cAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAat     >  1:3489216/1‑150 (MQ=255)
                                                                                                                                                    tGACCCGGCTCCGCGACATGGGCTTCTAATCATTCGTGCTTTCGTCTTCCCTGGATGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATatc  <  1:3892989/150‑1 (MQ=255)
                                                                                                                                                    tGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATatc  <  1:1611233/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATC  >  CP009273/3452604‑3452901

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: