Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,218,022 C→A 31.7% E110* (GAA→TAA)  pbpG ← D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,218,0220CA31.7% 527.3 / 235.0 286E110* (GAA→TAA) pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (93/102);  new base A (40/51);  total (133/153)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

GTGTTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCAG  >  CP009273/2217874‑2218169
                                                                                                                                                    |                                                                                                                                                   
gtgtTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTc                                                                                                                                                     <  2:3342782/149‑1 (MQ=255)
      ttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTAtt                                                                                                                                               >  2:871709/1‑149 (MQ=255)
      ttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTAtt                                                                                                                                               <  2:3145473/149‑1 (MQ=255)
       tCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTc                                                                                                                                              <  2:4076594/149‑1 (MQ=255)
       tCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATATCGCTATTc                                                                                                                                              >  2:1259695/1‑149 (MQ=255)
        cATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCa                                                                                                                                             <  1:504424/149‑1 (MQ=255)
          ttCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGt                                                                                                                                           >  2:7351/1‑149 (MQ=255)
           tCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTc                                                                                                                                          >  1:2186120/1‑149 (MQ=255)
            ccGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCg                                                                                                                                         <  2:667763/149‑1 (MQ=255)
              gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTa                                                                                                                                       <  1:4618684/149‑1 (MQ=255)
                gCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTAcg                                                                                                                                     >  2:3811429/1‑149 (MQ=255)
                gCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTAcg                                                                                                                                     <  1:1333755/149‑1 (MQ=255)
                gCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTAcg                                                                                                                                     >  2:2829117/1‑149 (MQ=255)
                gCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTAcg                                                                                                                                     >  2:1994840/1‑149 (MQ=255)
                  gATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTAcgcg                                                                                                                                   >  1:3920946/1‑149 (MQ=255)
                  gATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTAcgcg                                                                                                                                   >  1:4367197/1‑149 (MQ=255)
                       cgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAAt                                                                                                                              >  1:2335522/1‑149 (MQ=255)
                       cgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAAt                                                                                                                              >  1:4731424/1‑149 (MQ=255)
                       cgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAAt                                                                                                                              <  2:4485904/149‑1 (MQ=255)
                          tttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                           <  1:2903543/149‑1 (MQ=255)
                              gcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGAcccc                                                                                                                       >  2:2279274/1‑149 (MQ=255)
                               cattcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                      <  1:1789204/149‑1 (MQ=255)
                               cattcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                      <  1:1788520/149‑1 (MQ=255)
                                  tcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTCCGCGCGAATAGACCCCTTTc                                                                                                                   >  2:4605552/1‑149 (MQ=255)
                                   cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCa                                                                                                                  <  2:2230735/149‑1 (MQ=255)
                                   cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCa                                                                                                                  <  2:4355975/149‑1 (MQ=255)
                                   cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCa                                                                                                                  >  2:571583/1‑149 (MQ=255)
                                   cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCa                                                                                                                  >  1:3605400/1‑149 (MQ=255)
                                     ttGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                                <  1:736774/149‑1 (MQ=255)
                                     ttGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                                <  1:833302/149‑1 (MQ=255)
                                       gCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtctc                                                                                                              <  1:2980674/149‑1 (MQ=255)
                                          ttAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCgg                                                                                                           >  1:3717695/1‑149 (MQ=255)
                                           taataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGc                                                                                                          >  1:2489292/1‑149 (MQ=255)
                                             ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGt                                                                                                        >  2:2324472/1‑149 (MQ=255)
                                             ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGt                                                                                                        >  2:49528/1‑149 (MQ=255)
                                              taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTATCATCTCCGGCGTc                                                                                                       >  1:3810795/1‑149 (MQ=255)
                                              taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                       >  1:4568913/1‑149 (MQ=255)
                                              taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                       <  2:713760/149‑1 (MQ=255)
                                                 aGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTgg                                                                                                    >  2:3489617/1‑149 (MQ=255)
                                                  ggCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGc                                                                                                   >  1:941470/1‑149 (MQ=255)
                                                   gCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCt                                                                                                  <  2:3442289/149‑1 (MQ=255)
                                                   gCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCt                                                                                                  <  1:1476618/149‑1 (MQ=255)
                                                     tttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                                >  1:3059080/1‑149 (MQ=255)
                                                     tttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                                >  2:192922/1‑149 (MQ=255)
                                                         tAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATc                                                                                            >  1:2369504/1‑149 (MQ=255)
                                                          aaCCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATcc                                                                                           <  1:1287719/149‑1 (MQ=255)
                                                               ccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACttt                                                                                      <  2:119404/149‑1 (MQ=255)
                                                               ccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACttt                                                                                      <  1:3454514/149‑1 (MQ=255)
                                                                cGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACtttt                                                                                     <  1:4264239/149‑1 (MQ=255)
                                                                 gtgATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTTCGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                    <  1:1110369/147‑1 (MQ=255)
                                                                 gggATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                    >  2:2126515/1‑149 (MQ=255)
                                                                     tAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttt                                                                                <  1:3419044/149‑1 (MQ=255)
                                                                      aaTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtttt                                                                               >  1:2161158/1‑149 (MQ=255)
                                                                      aaTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTGTCATCTCCGGCGTCTGGCTGATCTCCACTTTTAGtttt                                                                               >  2:168028/1‑149 (MQ=255)
                                                                       aTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                              >  2:3941594/1‑149 (MQ=255)
                                                                        tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTc                                                                             <  1:2407481/149‑1 (MQ=255)
                                                                        tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTc                                                                             <  2:2383577/149‑1 (MQ=255)
                                                                            gcgcAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATcc                                                                         >  2:2533954/1‑149 (MQ=255)
                                                                              gcAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAg                                                                       >  2:1982373/1‑149 (MQ=255)
                                                                               cAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGc                                                                      <  2:4157267/149‑1 (MQ=255)
                                                                                aaGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCg                                                                     <  2:4480970/149‑1 (MQ=255)
                                                                                aaGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCg                                                                     <  1:3583882/149‑1 (MQ=255)
                                                                                   gCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCa                                                                  <  1:2558711/149‑1 (MQ=255)
                                                                                     ttGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGt                                                                <  1:4315148/149‑1 (MQ=255)
                                                                                      tGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGAGTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTc                                                               <  2:2441543/149‑1 (MQ=255)
                                                                                         cgccgcGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTATTTTTTCATCCAGCGGCAGTCGTg                                                            <  2:1314642/149‑1 (MQ=255)
                                                                                             gcgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATc                                                        >  1:579253/1‑149 (MQ=255)
                                                                                                 ggTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGc                                                    >  2:436482/1‑149 (MQ=255)
                                                                                                  gTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGca                                                   <  1:4125194/149‑1 (MQ=255)
                                                                                                     ttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACaa                                                <  1:2666860/149‑1 (MQ=255)
                                                                                                       cTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAAcc                                              >  1:3291194/1‑149 (MQ=255)
                                                                                                        tGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCa                                             <  1:2956725/149‑1 (MQ=255)
                                                                                                         gAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                            <  2:3690329/149‑1 (MQ=255)
                                                                                                         gAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                            >  1:1528549/1‑149 (MQ=255)
                                                                                                         gAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                            >  2:4409347/1‑149 (MQ=255)
                                                                                                         gAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                            >  2:259655/1‑149 (MQ=255)
                                                                                                            gACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgc                                         >  1:894512/1‑149 (MQ=255)
                                                                                                            gACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCGGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgc                                         >  2:2174265/1‑149 (MQ=255)
                                                                                                                 cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCa                                    <  2:1477373/149‑1 (MQ=255)
                                                                                                                 cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTGTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCa                                    >  1:4026024/1‑149 (MQ=255)
                                                                                                                  aGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCAt                                   >  1:1578192/1‑149 (MQ=255)
                                                                                                                   gcgcCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCAtt                                  <  2:1383383/149‑1 (MQ=255)
                                                                                                                       cAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAAtt                              <  2:3798118/149‑1 (MQ=255)
                                                                                                                       cAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAAGAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAAtt                              <  2:3030935/149‑1 (MQ=255)
                                                                                                                         gCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTg                            >  1:2565144/1‑149 (MQ=255)
                                                                                                                         gCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTg                            >  2:3088902/1‑149 (MQ=255)
                                                                                                                            aCAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTg                         >  2:3869850/1‑149 (MQ=255)
                                                                                                                              aGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGAt                       <  1:122932/149‑1 (MQ=255)
                                                                                                                              aGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGAt                       >  1:3156157/1‑149 (MQ=255)
                                                                                                                                cATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGACTAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAg                     <  2:2578906/149‑1 (MQ=255)
                                                                                                                                  tatCTGTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGAc                   <  2:1587043/149‑1 (MQ=255)
                                                                                                                                        tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATc             <  1:4109415/149‑1 (MQ=255)
                                                                                                                                         aCGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg            <  2:4731424/149‑1 (MQ=255)
                                                                                                                                         aCGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg            <  1:3564292/149‑1 (MQ=255)
                                                                                                                                              tGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGc       <  2:240912/149‑1 (MQ=255)
                                                                                                                                              tGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGc       >  2:802109/1‑149 (MQ=255)
                                                                                                                                                aTTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAc     >  2:2364324/1‑149 (MQ=255)
                                                                                                                                                aTTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAc     >  2:3610613/1‑149 (MQ=255)
                                                                                                                                                   tCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCAg  >  2:58082/1‑149 (MQ=255)
                                                                                                                                                   tCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCAg  >  2:3976819/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGTTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCAG  >  CP009273/2217874‑2218169

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: