Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 395,951 G→T 40.7% V143V (GTC→GTA ddlA ← D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273395,9510GT40.7% 261.0 / 257.2 221V143V (GTC→GTAddlAD‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base G (64/67);  new base T (53/37);  total (117/104)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

ACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACGGCTG  >  CP009273/395802‑396100
                                                                                                                                                     |                                                                                                                                                     
  agaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGAc                                                                                                                                                     <  1:167456/150‑1 (MQ=255)
   gaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACa                                                                                                                                                    >  1:755243/1‑150 (MQ=255)
   gaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACa                                                                                                                                                    <  1:657481/150‑1 (MQ=255)
       agCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCtt                                                                                                                                                >  2:2056866/1‑150 (MQ=255)
        gCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCttt                                                                                                                                               >  2:4085409/1‑150 (MQ=255)
        gCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCttt                                                                                                                                               >  1:3779067/1‑150 (MQ=255)
        gCCCTGATTAGCCGGGTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCttt                                                                                                                                               >  2:2496644/1‑150 (MQ=255)
         cccTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTg                                                                                                                                              >  1:3051098/1‑150 (MQ=255)
          ccTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGt                                                                                                                                             <  2:320685/150‑1 (MQ=255)
          ccTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGt                                                                                                                                             >  1:576217/1‑150 (MQ=255)
           cTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAATACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTc                                                                                                                                            >  2:1938052/1‑150 (MQ=255)
               ttAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATa                                                                                                                                        <  1:892630/150‑1 (MQ=255)
                 aGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACa                                                                                                                                      >  2:942792/1‑150 (MQ=255)
                  gCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAg                                                                                                                                     <  1:2450231/150‑1 (MQ=255)
                    cGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGgc                                                                                                                                   <  2:2246324/150‑1 (MQ=255)
                         ttACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGctgctg                                                                                                                              >  1:488859/1‑150 (MQ=255)
                         ttACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGctgctg                                                                                                                              <  1:2916753/150‑1 (MQ=255)
                           aCAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGaa                                                                                                                            >  1:1090480/1‑150 (MQ=255)
                               aCAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCa                                                                                                                        >  1:2876308/1‑150 (MQ=255)
                                cAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAg                                                                                                                       <  1:2209466/150‑1 (MQ=255)
                                   tGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAAc                                                                                                                    >  1:1857480/1‑150 (MQ=255)
                                       aaCCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCa                                                                                                                >  1:1823681/1‑150 (MQ=255)
                                        aCCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAg                                                                                                               <  1:3437577/150‑1 (MQ=255)
                                               tttAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTTCATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTAc                                                                                                        <  1:2586884/150‑1 (MQ=255)
                                               tttAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGCTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTAc                                                                                                        <  1:3803839/150‑1 (MQ=255)
                                                 tAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACaa                                                                                                      <  1:2684509/150‑1 (MQ=255)
                                                    aCTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACg                                                                                                   <  2:370226/150‑1 (MQ=255)
                                                           ttCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAAtt                                                                                            <  1:2361447/150‑1 (MQ=255)
                                                            tCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTg                                                                                           <  2:2844073/150‑1 (MQ=255)
                                                                 aaaaCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGAc                                                                                      >  2:736089/1‑150 (MQ=255)
                                                                 aaaaCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGAc                                                                                      <  2:2690966/150‑1 (MQ=255)
                                                                  aaaCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGAcc                                                                                     >  2:1781761/1‑150 (MQ=255)
                                                                   aaCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGAccc                                                                                    >  2:682721/1‑150 (MQ=255)
                                                                   aaCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGAccc                                                                                    >  2:802511/1‑150 (MQ=255)
                                                                      tAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCgca                                                                                 >  2:1715509/1‑150 (MQ=255)
                                                                      tAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCgca                                                                                 <  1:2005287/150‑1 (MQ=255)
                                                                          gttgtGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCAt                                                                             <  2:1230737/150‑1 (MQ=255)
                                                                          gttgtGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCAt                                                                             <  1:1133178/150‑1 (MQ=255)
                                                                          gttgtGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCAt                                                                             <  2:1047040/150‑1 (MQ=255)
                                                                             gtgACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTcc                                                                          <  1:3160668/150‑1 (MQ=255)
                                                                              tgACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTccc                                                                         >  2:3417906/1‑150 (MQ=255)
                                                                                 cGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGc                                                                      >  1:1160993/1‑150 (MQ=255)
                                                                                     tGGCGCGCGTCGGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAgg                                                                  >  1:1788375/1‑150 (MQ=255)
                                                                                     tGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAgg                                                                  <  2:620112/150‑1 (MQ=255)
                                                                                      ggCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGa                                                                 <  2:2333954/150‑1 (MQ=255)
                                                                                      ggCGAGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCCGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGa                                                                 <  2:1824572/150‑1 (MQ=255)
                                                                                          cgcgTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCa                                                             >  2:2131603/1‑150 (MQ=255)
                                                                                          cgcgTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCa                                                             <  2:622779/150‑1 (MQ=255)
                                                                                              tCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAAGTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCtt                                                         <  2:841019/150‑1 (MQ=255)
                                                                                               cAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTc                                                        <  1:1433891/150‑1 (MQ=255)
                                                                                               cAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTc                                                        <  2:3978129/150‑1 (MQ=255)
                                                                                                aGGGTAATAACTGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCg                                                       <  1:706431/150‑1 (MQ=255)
                                                                                                  ggTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGcc                                                     <  1:3913791/150‑1 (MQ=255)
                                                                                                   gTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGccc                                                    <  2:1768176/150‑1 (MQ=255)
                                                                                                    taataaATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCa                                                   >  2:3586878/1‑150 (MQ=255)
                                                                                                    taataaATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCa                                                   <  2:1477149/150‑1 (MQ=255)
                                                                                                     aataaATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAg                                                  >  2:3461892/1‑150 (MQ=255)
                                                                                                     aataaATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAg                                                  >  2:1011134/1‑150 (MQ=255)
                                                                                                          aTGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTAc                                             <  1:1666295/150‑1 (MQ=255)
                                                                                                          aTGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTAc                                             >  1:1554811/1‑150 (MQ=255)
                                                                                                           tGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTAcc                                            <  1:3357363/150‑1 (MQ=255)
                                                                                                           tGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTAcc                                            <  1:3357251/150‑1 (MQ=255)
                                                                                                             gcgcAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGt                                          <  1:3086762/150‑1 (MQ=255)
                                                                                                                 aaTGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGAc                                      >  1:3287941/1‑150 (MQ=255)
                                                                                                                 aaTGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGAc                                      >  2:153716/1‑150 (MQ=255)
                                                                                                                 aaTGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGAc                                      >  1:3807299/1‑150 (MQ=255)
                                                                                                                   tGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACaa                                    <  2:506759/150‑1 (MQ=255)
                                                                                                                    gTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAAt                                   <  1:611267/150‑1 (MQ=255)
                                                                                                                    gTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAAt                                   >  1:3094592/1‑150 (MQ=255)
                                                                                                                     ttCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATc                                  <  2:662805/150‑1 (MQ=255)
                                                                                                                     ttCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATc                                  <  2:610692/150‑1 (MQ=255)
                                                                                                                      tCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCg                                 <  2:2208534/150‑1 (MQ=255)
                                                                                                                       cAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCgg                                <  2:4151650/150‑1 (MQ=255)
                                                                                                                          cccGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGaa                             <  1:138494/150‑1 (MQ=255)
                                                                                                                          cccGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGaa                             <  1:2153147/150‑1 (MQ=255)
                                                                                                                            cGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGaaaa                           >  1:1491739/1‑150 (MQ=255)
                                                                                                                            cGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCAAACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGaaaa                           >  2:3474574/1‑150 (MQ=255)
                                                                                                                             ggCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAAt                          <  1:103638/150‑1 (MQ=255)
                                                                                                                                 tCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACa                      >  2:611314/1‑150 (MQ=255)
                                                                                                                                     gTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCa                  <  2:2632780/150‑1 (MQ=255)
                                                                                                                                      tAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCAc                 >  1:4075384/1‑150 (MQ=255)
                                                                                                                                       aGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCAcc                <  1:767987/150‑1 (MQ=255)
                                                                                                                                       aGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCAcc                >  2:3106902/1‑150 (MQ=255)
                                                                                                                                         cAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGt              >  2:3580946/1‑150 (MQ=255)
                                                                                                                                         cAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGt              >  1:3665300/1‑150 (MQ=255)
                                                                                                                                          aGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGACAATGACATCCACCGTc             >  2:1413901/1‑150 (MQ=255)
                                                                                                                                            aCGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCgg           <  2:3323986/150‑1 (MQ=255)
                                                                                                                                            aCGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCgg           <  1:811825/150‑1 (MQ=255)
                                                                                                                                             cGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGc          >  2:2423568/1‑150 (MQ=255)
                                                                                                                                              gTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCa         >  2:3353504/1‑150 (MQ=255)
                                                                                                                                               tttGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCaa        <  1:2835614/150‑1 (MQ=255)
                                                                                                                                                ttGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAAc       >  2:3029850/1‑150 (MQ=255)
                                                                                                                                                ttGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAAc       >  1:2374133/1‑150 (MQ=255)
                                                                                                                                                 tGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACg      >  1:69155/1‑150 (MQ=255)
                                                                                                                                                  ggTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACgg     >  2:3812488/1‑150 (MQ=255)
                                                                                                                                                  ggTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAACCGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACgg     <  2:2964356/150‑1 (MQ=255)
                                                                                                                                                  ggTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACgg     >  2:389402/1‑150 (MQ=255)
                                                                                                                                                  ggTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACgg     <  2:3515570/150‑1 (MQ=255)
                                                                                                                                                  ggTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACgg     <  2:875517/150‑1 (MQ=255)
                                                                                                                                                     gACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACGGCTg  <  2:1857133/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCAAGGAACCATCTTCGCCCAGCGTACCGTGGACAATCGGGAAAATGACATCCACCGTCGGCAACGGCTG  >  CP009273/395802‑396100

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: