Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 395,611 C→A 21.6% D257Y (GAC→TAC)  ddlA ← D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273395,6110CA21.6% 728.1 / 162.7 276D257Y (GAC→TAC) ddlAD‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base C (96/120);  new base A (21/39);  total (117/159)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGAC  >  CP009273/395462‑395758
                                                                                                                                                     |                                                                                                                                                   
   tAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTa                                                                                                                                                  <  1:475308/150‑1 (MQ=255)
   tAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTa                                                                                                                                                  >  2:2884823/1‑150 (MQ=255)
       aaCACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCa                                                                                                                                              >  2:1841720/1‑150 (MQ=255)
        acacGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCAt                                                                                                                                             <  1:2805269/150‑1 (MQ=255)
        acacGTCTACACGCGCCATGCCCGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCAt                                                                                                                                             <  1:1966482/150‑1 (MQ=255)
           cGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGa                                                                                                                                          >  2:4761110/1‑150 (MQ=255)
              cTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAAt                                                                                                                                       >  2:4798570/1‑150 (MQ=255)
                acacGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCg                                                                                                                                     <  1:4782016/150‑1 (MQ=255)
                  acGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCt                                                                                                                                   <  2:464459/150‑1 (MQ=255)
                  acGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCATAGGCATAGAAATCGCt                                                                                                                                   >  1:1879128/1‑150 (MQ=255)
                   cgcgCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTg                                                                                                                                  <  2:189636/150‑1 (MQ=255)
                       ccATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGa                                                                                                                              >  1:1442315/1‑150 (MQ=255)
                       ccATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGa                                                                                                                              >  2:2553109/1‑150 (MQ=255)
                               gcgcATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGAtc                                                                                                                      <  2:3044753/150‑1 (MQ=255)
                                cgcATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGAtct                                                                                                                     >  1:4597823/1‑150 (MQ=255)
                                cgcATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGAtct                                                                                                                     <  1:3926367/150‑1 (MQ=255)
                                  cATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGAtctct                                                                                                                   <  2:1189414/150‑1 (MQ=255)
                                    tCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTcc                                                                                                                 >  1:4498945/1‑150 (MQ=255)
                                    tCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTcc                                                                                                                 >  1:4497815/1‑150 (MQ=255)
                                    tCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTcc                                                                                                                 <  2:4458728/150‑1 (MQ=255)
                                      ccAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCcac                                                                                                               >  2:3566075/1‑150 (MQ=255)
                                      ccAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCcac                                                                                                               >  2:2926949/1‑150 (MQ=255)
                                      ccAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCcac                                                                                                               <  2:176377/150‑1 (MQ=255)
                                        aaCGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAg                                                                                                             >  2:4139970/1‑150 (MQ=255)
                                        aaCGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAg                                                                                                             >  1:4058848/1‑150 (MQ=255)
                                        aaCGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAg                                                                                                             >  1:1601139/1‑150 (MQ=255)
                                          cGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGt                                                                                                           <  1:1650740/150‑1 (MQ=255)
                                              tGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTg                                                                                                       <  1:488431/150‑1 (MQ=255)
                                              tGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTg                                                                                                       <  2:3262256/150‑1 (MQ=255)
                                                aTAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGc                                                                                                     >  2:3869274/1‑150 (MQ=255)
                                                 tAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCt                                                                                                    >  1:3065446/1‑150 (MQ=255)
                                                  aaGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCtt                                                                                                   <  1:4586692/150‑1 (MQ=255)
                                                  aaGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCtt                                                                                                   <  1:2434502/150‑1 (MQ=255)
                                                  aaGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCtt                                                                                                   >  2:3389503/1‑150 (MQ=255)
                                                     ccTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCg                                                                                                >  1:3662480/1‑150 (MQ=255)
                                                     ccTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCg                                                                                                <  2:828931/150‑1 (MQ=255)
                                                      cTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCgg                                                                                               >  1:3855654/1‑150 (MQ=255)
                                                      cTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCgg                                                                                               >  1:2470594/1‑150 (MQ=255)
                                                        gAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGAt                                                                                             <  1:3714235/150‑1 (MQ=255)
                                                        gAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGAt                                                                                             >  2:506457/1‑150 (MQ=255)
                                                         aaCGGCACTCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGAtt                                                                                            <  2:409820/150‑1 (MQ=255)
                                                          aCGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTg                                                                                           <  1:102628/150‑1 (MQ=255)
                                                            ggCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTc                                                                                         >  2:1618300/1‑150 (MQ=255)
                                                            ggCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTc                                                                                         <  1:1393602/150‑1 (MQ=255)
                                                            ggCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTc                                                                                         >  1:763099/1‑150 (MQ=255)
                                                              cAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGt                                                                                       <  2:284365/150‑1 (MQ=255)
                                                              cAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGt                                                                                       <  2:1368989/150‑1 (MQ=255)
                                                               aaTCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCTGAACTAACACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATTGAAATCGCTGGTGTGTACGCTCTCTCAAAAGGTGCTGGCTTGCGGATTGTCGtt                                                                                      <  1:2627643/150‑1 (MQ=255)
                                                                 tCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTAGTTGc                                                                                    <  2:1314687/150‑1 (MQ=255)
                                                                    cccGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCa                                                                                 >  1:1796747/1‑150 (MQ=255)
                                                                       ggATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGaa                                                                              >  2:50851/1‑150 (MQ=255)
                                                                          tCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTg                                                                           <  1:3836246/150‑1 (MQ=255)
                                                                            ttATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTgcg                                                                         <  1:2886966/150‑1 (MQ=255)
                                                                               tCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCAt                                                                      <  2:4575728/150‑1 (MQ=255)
                                                                                cGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCAtt                                                                     <  2:3097549/150‑1 (MQ=255)
                                                                                 gTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTc                                                                    <  1:2767541/150‑1 (MQ=255)
                                                                                     aTTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGaac                                                                >  1:916059/1‑148 (MQ=255)
                                                                                     aTTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGAtc                                                                <  2:2244786/150‑1 (MQ=255)
                                                                                      tttCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGAtct                                                               <  1:652127/150‑1 (MQ=255)
                                                                                         cTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCAc                                                            <  2:3617649/150‑1 (MQ=255)
                                                                                           ggCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGa                                                          >  2:4034547/1‑150 (MQ=255)
                                                                                            gcgcAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGAc                                                         <  1:556717/150‑1 (MQ=255)
                                                                                              gcAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCt                                                       <  1:1269385/150‑1 (MQ=255)
                                                                                              gcAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCt                                                       >  1:3854509/1‑150 (MQ=255)
                                                                                               cAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCtt                                                      <  1:204112/150‑1 (MQ=255)
                                                                                                 aTGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTg                                                    <  1:2491162/150‑1 (MQ=255)
                                                                                                  tGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGa                                                   >  2:2819964/1‑150 (MQ=255)
                                                                                                    gCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATc                                                 <  1:3293579/150‑1 (MQ=255)
                                                                                                    gCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATc                                                 <  2:709078/150‑1 (MQ=255)
                                                                                                     cTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATcc                                                <  1:2119627/150‑1 (MQ=255)
                                                                                                     cTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATcc                                                <  1:2117956/150‑1 (MQ=255)
                                                                                                     cTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATcc                                                <  1:1090619/150‑1 (MQ=255)
                                                                                                     cTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATcc                                                <  1:4139970/150‑1 (MQ=255)
                                                                                                      tGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATccc                                               <  1:3972112/150‑1 (MQ=255)
                                                                                                         cGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTg                                            <  1:4385032/150‑1 (MQ=255)
                                                                                                         cGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTg                                            <  1:4047138/150‑1 (MQ=255)
                                                                                                          ggAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGc                                           >  2:1489384/1‑150 (MQ=255)
                                                                                                          ggAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGc                                           >  2:2018901/1‑150 (MQ=255)
                                                                                                            aaCTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTc                                         <  2:1223204/150‑1 (MQ=255)
                                                                                                              cTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCaa                                       >  1:951566/1‑150 (MQ=255)
                                                                                                                aCCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACg                                     <  2:2413114/150‑1 (MQ=255)
                                                                                                                 ccACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGa                                    <  1:495745/150‑1 (MQ=255)
                                                                                                                  cACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGAt                                   >  1:2559742/1‑150 (MQ=255)
                                                                                                                   aCTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATc                                  <  1:3530588/150‑1 (MQ=255)
                                                                                                                    cTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCa                                 <  2:3228951/150‑1 (MQ=255)
                                                                                                                    cTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCa                                 >  2:3726806/1‑150 (MQ=255)
                                                                                                                     tttCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCAc                                >  2:2133063/1‑150 (MQ=255)
                                                                                                                          cgcCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTAt                           <  1:3206954/150‑1 (MQ=255)
                                                                                                                                cTTCGTCAATGTACTTGGTGTCGTAGGCGTAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcga                     >  1:1529802/1‑150 (MQ=255)
                                                                                                                                   cGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaact                  <  2:3442543/150‑1 (MQ=255)
                                                                                                                                      cAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcga               >  2:260009/1‑150 (MQ=255)
                                                                                                                                      cAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcga               >  1:4125311/1‑150 (MQ=255)
                                                                                                                                      cAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcga               >  2:2675662/1‑150 (MQ=255)
                                                                                                                                        aTGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcgaac             >  2:151360/1‑150 (MQ=255)
                                                                                                                                        aTGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCACAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcgaac             >  2:1019127/1‑150 (MQ=255)
                                                                                                                                          gTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGc           >  2:63120/1‑150 (MQ=255)
                                                                                                                                              ttGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGa       >  1:583623/1‑150 (MQ=255)
                                                                                                                                               tGGTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGAt      >  2:4671230/1‑150 (MQ=255)
                                                                                                                                                ggTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATc     >  1:1943948/1‑150 (MQ=255)
                                                                                                                                                ggTGTAGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATc     <  2:1775911/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGAC  >  CP009273/395462‑395758

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: