Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 241,630 G→T 7.9% Y221* (TAC→TAA)  yafK ← L,D‑transpeptidase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273241,6300GT 7.9% 1062.9 / 46.4 280Y221* (TAC→TAA) yafKL,D‑transpeptidase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (136/121);  new base T (11/12);  total (147/133)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

ACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATAC  >  CP009273/241483‑241778
                                                                                                                                                   |                                                                                                                                                    
 cGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTa                                                                                                                                                    <  1:3490399/149‑1 (MQ=255)
    tGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACg                                                                                                                                                 >  2:1975279/1‑149 (MQ=255)
     gATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGa                                                                                                                                                >  1:3671517/1‑149 (MQ=255)
       ttGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGAcc                                                                                                                                              >  1:1264130/1‑149 (MQ=255)
         gCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAt                                                                                                                                            >  1:4783390/1‑149 (MQ=255)
         gCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAt                                                                                                                                            <  2:2334171/149‑1 (MQ=255)
         gCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAt                                                                                                                                            >  1:3558188/1‑149 (MQ=255)
          cAACACAGCCGCTGACGGGCAATGAACCCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAtt                                                                                                                                           <  1:2316426/149‑1 (MQ=255)
          cAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAtt                                                                                                                                           <  2:4891364/149‑1 (MQ=255)
          cAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAtt                                                                                                                                           <  2:2289285/149‑1 (MQ=255)
          cAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCAtt                                                                                                                                           >  2:2989716/1‑149 (MQ=255)
           aaCACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTg                                                                                                                                          <  2:425820/149‑1 (MQ=255)
           aaCACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTg                                                                                                                                          >  2:3127545/1‑149 (MQ=255)
                  ccGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAg                                                                                                                                   >  2:3361166/1‑149 (MQ=255)
                    gCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGaa                                                                                                                                 >  1:20428/1‑149 (MQ=255)
                    gCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGaa                                                                                                                                 >  1:21191/1‑149 (MQ=255)
                     cTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGaaa                                                                                                                                <  1:1965340/149‑1 (MQ=255)
                        aCGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACgg                                                                                                                             >  2:501294/1‑149 (MQ=255)
                        aCGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACgg                                                                                                                             <  1:3036740/149‑1 (MQ=255)
                        aCGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACgg                                                                                                                             <  2:761971/149‑1 (MQ=255)
                         cGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGt                                                                                                                            <  1:1639513/149‑1 (MQ=255)
                          gggCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGtt                                                                                                                           >  1:4575595/1‑149 (MQ=255)
                             cAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTtggt                                                                                                                        <  2:4610449/149‑1 (MQ=255)
                               aTGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTtggtgg                                                                                                                      >  1:901131/1‑149 (MQ=255)
                               aTGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTtggtgg                                                                                                                      >  1:1646148/1‑149 (MQ=255)
                                  aaaCCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCtt                                                                                                                   <  1:557479/149‑1 (MQ=255)
                                  aaaCCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCtt                                                                                                                   >  2:708239/1‑149 (MQ=255)
                                   aaCCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTa                                                                                                                  <  2:2629675/149‑1 (MQ=255)
                                   aaCCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTa                                                                                                                  >  2:4346056/1‑149 (MQ=255)
                                    aCCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTGTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTAc                                                                                                                 <  2:4758578/149‑1 (MQ=255)
                                      cTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACgg                                                                                                               >  1:2427449/1‑149 (MQ=255)
                                      cTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACgg                                                                                                               <  1:2063193/149‑1 (MQ=255)
                                        ggCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGATGGTGGCTTACGGGt                                                                                                             >  2:2795458/1‑149 (MQ=255)
                                         gCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGtt                                                                                                            <  2:544699/149‑1 (MQ=255)
                                           aaaaGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTg                                                                                                          >  1:1163514/1‑149 (MQ=255)
                                            aaaGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGc                                                                                                         <  2:1772643/149‑1 (MQ=255)
                                             aaGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCt                                                                                                        <  1:3116703/149‑1 (MQ=255)
                                                 ttATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCaaa                                                                                                    >  2:494318/1‑149 (MQ=255)
                                                   aTGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAagt                                                                                                  <  1:4193078/149‑1 (MQ=255)
                                                    tGCCAGGCGGATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAagta                                                                                                 >  1:2678844/1‑149 (MQ=255)
                                                    tGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAagta                                                                                                 >  1:3174294/1‑149 (MQ=255)
                                                     gCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAagtag                                                                                                <  2:3633674/149‑1 (MQ=255)
                                                     gCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAagtag                                                                                                >  1:340637/1‑149 (MQ=255)
                                                       cAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTc                                                                                              >  1:2973870/1‑149 (MQ=255)
                                                        aGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCg                                                                                             <  2:1779514/149‑1 (MQ=255)
                                                         ggCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGt                                                                                            >  1:3611170/1‑149 (MQ=255)
                                                          gCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTa                                                                                           <  1:3022639/149‑1 (MQ=255)
                                                           cGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTaa                                                                                          <  1:3233102/149‑1 (MQ=255)
                                                              aTGGCGCTTATTTTGCCGCGGGGAGCGTGTAGTTTGATTCCAGTTGTGGCTGCACAACATCGTGGCTTAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAAccc                                                                                       >  1:154251/1‑149 (MQ=255)
                                                                  cgcTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCt                                                                                   <  2:2708976/149‑1 (MQ=255)
                                                                      tATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAg                                                                               >  1:44452/1‑149 (MQ=255)
                                                                       aTTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGt                                                                              <  1:404142/149‑1 (MQ=255)
                                                                        ttttGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGtt                                                                             >  1:100429/1‑149 (MQ=255)
                                                                           tGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTCGGTGGCTTACGGGTTTGCTC‑AAGTAGTCGTAACCCGGCTTCAGTTGCTc                                                                         <  1:807865/149‑1 (MQ=255)
                                                                              cTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTccc                                                                       <  1:1230943/149‑1 (MQ=255)
                                                                                 ggggAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGa                                                                    >  2:440206/1‑149 (MQ=255)
                                                                                  gggAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGaa                                                                   >  2:1135388/1‑149 (MQ=255)
                                                                                    gAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGt                                                                 <  2:588108/149‑1 (MQ=255)
                                                                                     aGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTc                                                                <  2:1066052/149‑1 (MQ=255)
                                                                                       cGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCt                                                              <  1:3768160/149‑1 (MQ=255)
                                                                                            aGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAg                                                         <  2:3985784/149‑1 (MQ=255)
                                                                                             gTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGt                                                        <  1:708239/149‑1 (MQ=255)
                                                                                             gTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGt                                                        >  2:5162244/1‑149 (MQ=255)
                                                                                             gTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGt                                                        <  2:619703/149‑1 (MQ=255)
                                                                                             gTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGt                                                        >  1:3465447/1‑149 (MQ=255)
                                                                                              tttGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGtt                                                       <  1:3834739/149‑1 (MQ=255)
                                                                                              tttGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGtt                                                       <  2:3645422/149‑1 (MQ=255)
                                                                                                   tGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAat                                                  <  1:143179/149‑1 (MQ=255)
                                                                                                    gCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAata                                                 >  1:2214821/1‑149 (MQ=255)
                                                                                                    gCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAata                                                 <  2:3892846/149‑1 (MQ=255)
                                                                                                    gCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAata                                                 <  1:3990485/149‑1 (MQ=255)
                                                                                                      cAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATAtt                                               >  2:5175315/1‑149 (MQ=255)
                                                                                                       aGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATAttt                                              <  2:1550499/149‑1 (MQ=255)
                                                                                                       aGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATAttt                                              >  1:4503347/1‑149 (MQ=255)
                                                                                                        gTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTa                                             <  1:2742102/149‑1 (MQ=255)
                                                                                                         ttGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTAt                                            >  1:4905192/1‑149 (MQ=255)
                                                                                                           gtgGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATgg                                          <  2:1487900/149‑1 (MQ=255)
                                                                                                                  cacaACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCa                                   <  1:492284/149‑1 (MQ=255)
                                                                                                                  cacaACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCa                                   <  1:226789/149‑1 (MQ=255)
                                                                                                                    caACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCata                                 >  1:3525340/1‑149 (MQ=255)
                                                                                                                     aaCTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCatat                                <  1:1464730/149‑1 (MQ=255)
                                                                                                                      aCTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATAtt                               >  1:2447787/1‑149 (MQ=255)
                                                                                                                        ttCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTgg                             >  2:2408486/1‑149 (MQ=255)
                                                                                                                         tCGTGGCTTAACGGCTTGCTGACCACGTCACGACCATTGACCACAGAAACGGTTGGGGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGc                            <  1:4106036/149‑1 (MQ=255)
                                                                                                                         tCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGc                            <  1:1312155/149‑1 (MQ=255)
                                                                                                                           gTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGt                          >  1:4785775/1‑149 (MQ=255)
                                                                                                                              gCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCgg                       >  2:790820/1‑149 (MQ=255)
                                                                                                                                ttAACGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTc                     <  2:4038622/149‑1 (MQ=255)
                                                                                                                                 ttaCGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTTGTGGCTTACGTGTTTTCTCACAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCa                    <  1:4380938/147‑1 (MQ=255)
                                                                                                                                   aCGGCTTGCTGACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATg                  >  2:1954324/1‑149 (MQ=255)
                                                                                                                                   aCGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATg                  >  2:4115807/1‑149 (MQ=255)
                                                                                                                                    cGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGc                 >  2:691112/1‑149 (MQ=255)
                                                                                                                                        ttGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCGTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGaa             >  2:3464467/1‑149 (MQ=255)
                                                                                                                                          gCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCAGATTGGCGTCGGTCATGCGGAACg           >  2:2529929/1‑149 (MQ=255)
                                                                                                                                            tGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACggg         <  2:3796314/149‑1 (MQ=255)
                                                                                                                                             gACCACTTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGt        >  2:1731835/1‑149 (MQ=255)
                                                                                                                                                cACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTtaa     <  2:3674548/149‑1 (MQ=255)
                                                                                                                                                cACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTaaa     <  2:4503111/149‑1 (MQ=255)
                                                                                                                                                 aCGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAAt    >  2:5131826/1‑149 (MQ=255)
                                                                                                                                                  cGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCTGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATa   <  2:4541500/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATAC  >  CP009273/241483‑241778

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: ‑