Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,518,052 C→A 12.2% P608Q (CCG→CAG)  mrcA → fused penicillin‑binding protein 1a: murein transglycosylase/murein transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,518,0520CA12.2% 753.3 / 52.6 228P608Q (CCG→CAG) mrcAfused penicillin‑binding protein 1a: murein transglycosylase/murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (98/102);  new base A (12/16);  total (110/118)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 228 total aligned reads displayed.

CGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGAGCAG  >  CP009273/3517908‑3518197
                                                                                                                                                |                                                                                                                                                 
cgctggcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaa                                                                                                                                                >  2:3503876/1‑148 (MQ=255)
cgctggcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaa                                                                                                                                                <  2:1738224/148‑1 (MQ=255)
   tggcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGt                                                                                                                                             <  2:459157/148‑1 (MQ=255)
      cgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGc                                                                                                                                          <  2:800785/148‑1 (MQ=255)
          gggTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGc                                                                                                                                      >  1:3562086/1‑148 (MQ=255)
           ggTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGcc                                                                                                                                     >  2:2386955/1‑148 (MQ=255)
               cAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGa                                                                                                                                 <  2:2311130/148‑1 (MQ=255)
                 gCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAAt                                                                                                                               >  1:879540/1‑148 (MQ=255)
                   gTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGc                                                                                                                             >  1:719707/1‑148 (MQ=255)
                     ccTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGa                                                                                                                           >  2:2023804/1‑148 (MQ=255)
                      cTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGat                                                                                                                          <  1:1497307/148‑1 (MQ=255)
                      cTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGat                                                                                                                          <  2:458250/148‑1 (MQ=255)
                       ttCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGata                                                                                                                         <  2:195740/148‑1 (MQ=255)
                       ttCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGata                                                                                                                         <  2:2125346/148‑1 (MQ=255)
                         cACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATAtt                                                                                                                       <  2:3817632/148‑1 (MQ=255)
                          aCCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTc                                                                                                                      >  1:3686168/1‑148 (MQ=255)
                            cccAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCg                                                                                                                            >  2:3467465/1‑140 (MQ=255)
                            cccAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCg                                                                                                                            <  1:3467465/140‑1 (MQ=255)
                             ccAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCgg                                                                                                                   >  2:1364542/1‑148 (MQ=255)
                               aaTGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTg                                                                                                                 >  2:1617076/1‑148 (MQ=255)
                               aaTGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTg                                                                                                                 <  1:2233275/148‑1 (MQ=255)
                                aTGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGa                                                                                                                >  1:1822194/1‑148 (MQ=255)
                                    aGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTa                                                                                                            <  2:849133/148‑1 (MQ=255)
                                    aGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTa                                                                                                            <  2:1417713/148‑1 (MQ=255)
                                      gTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTCCg                                                                                                          >  2:2134120/1‑148 (MQ=255)
                                        ggCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGt                                                                                                        <  2:3390203/148‑1 (MQ=255)
                                         gcgcgcGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTg                                                                                                       <  2:1202993/148‑1 (MQ=255)
                                         gcgcgcGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTg                                                                                                       >  2:691383/1‑148 (MQ=255)
                                              cgGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAACCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACg                                                                                                  <  2:591437/148‑1 (MQ=255)
                                               ggCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGc                                                                                                 <  1:500287/148‑1 (MQ=255)
                                                   aCGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGaa                                                                                             <  1:1253370/148‑1 (MQ=255)
                                                    cgcgGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGaaa                                                                                            <  1:3355460/148‑1 (MQ=255)
                                                        gTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGa                                                                                        <  2:3808786/148‑1 (MQ=255)
                                                         tCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGaa                                                                                       <  2:3686168/148‑1 (MQ=255)
                                                         tCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGaa                                                                                       <  1:953322/148‑1 (MQ=255)
                                                            tGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGt                                                                                    >  2:2823011/1‑148 (MQ=255)
                                                             ggCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTg                                                                                   >  2:2225048/1‑148 (MQ=255)
                                                                   cggcggCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGaa                                                                             <  2:1701335/148‑1 (MQ=255)
                                                                    ggcggcTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGaaa                                                                            <  2:2803830/148‑1 (MQ=255)
                                                                    ggcggcTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGaaa                                                                            <  2:3450392/148‑1 (MQ=255)
                                                                    ggcggcTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGaaa                                                                            <  2:2876424/148‑1 (MQ=255)
                                                                        gcTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATaa                                                                        <  1:3379758/148‑1 (MQ=255)
                                                                           tCCTGGTGGACCCGTGGTTTATCATCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGAGAAACCGAAAGTAGCCTCCCCGGAATGCGATAATCCGGTGATTAACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGa                                                                     >  1:2525978/1‑148 (MQ=255)
                                                                           tCCTGGGGGACCCGTGGTTGATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGa                                                                     <  2:2346512/148‑1 (MQ=255)
                                                                             cTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATg                                                                   <  1:58993/148‑1 (MQ=255)
                                                                             cTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATg                                                                   <  1:2996807/148‑1 (MQ=255)
                                                                               ggtggACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGtt                                                                 >  2:3244783/1‑148 (MQ=255)
                                                                                   gACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAg                                                             >  2:3638961/1‑148 (MQ=255)
                                                                                      ccGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATg                                                          <  1:3503876/148‑1 (MQ=255)
                                                                                       cGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGt                                                         >  1:2786776/1‑148 (MQ=255)
                                                                                       cGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGAGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGt                                                         <  1:1569075/148‑1 (MQ=255)
                                                                                        gTGGTTTATCAGCAAAATTGAAAACAATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTc                                                        >  2:3060941/1‑148 (MQ=255)
                                                                                         tGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCg                                                       <  1:3001301/148‑1 (MQ=255)
                                                                                         tGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCg                                                       <  1:3181071/148‑1 (MQ=255)
                                                                                           gTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCt                                                     >  2:1171466/1‑148 (MQ=255)
                                                                                           gTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCt                                                     >  2:1889934/1‑148 (MQ=255)
                                                                                             ttATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTAt                                                   >  1:1614634/1‑148 (MQ=255)
                                                                                                 cAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTcc                                               <  1:2034471/148‑1 (MQ=255)
                                                                                                 cAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTcc                                               >  2:2547945/1‑148 (MQ=255)
                                                                                                    cAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCcgc                                            >  1:1947375/1‑148 (MQ=255)
                                                                                                      aaaTTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGa                                          >  2:646850/1‑148 (MQ=255)
                                                                                                        aTTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGagc                                        >  1:2116881/1‑148 (MQ=255)
                                                                                                        aTTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGagc                                        <  1:11087/148‑1 (MQ=255)
                                                                                                        aTTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGagc                                        <  1:1033861/148‑1 (MQ=255)
                                                                                                         ttGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGagca                                       >  2:3124425/1‑148 (MQ=255)
                                                                                                         ttGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGagca                                       <  2:3197461/148‑1 (MQ=255)
                                                                                                             aaaCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGcagcag                                   >  1:210441/1‑148 (MQ=255)
                                                                                                              aaCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGa                                  <  2:617762/148‑1 (MQ=255)
                                                                                                              aaCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGa                                  <  1:1812200/148‑1 (MQ=255)
                                                                                                              aaCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGa                                  <  1:1529327/148‑1 (MQ=255)
                                                                                                               aCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGaa                                 >  2:1038002/1‑148 (MQ=255)
                                                                                                               aCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGaa                                 >  2:2472931/1‑148 (MQ=255)
                                                                                                               aCGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGaa                                 >  1:1204188/1‑148 (MQ=255)
                                                                                                                 gATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATg                               <  2:1406333/148‑1 (MQ=255)
                                                                                                                 gATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATg                               <  1:1877818/148‑1 (MQ=255)
                                                                                                                  aTCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGt                              <  2:1294503/148‑1 (MQ=255)
                                                                                                                  aTCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGt                              <  2:3790859/148‑1 (MQ=255)
                                                                                                                    cAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGttt                            <  1:1089184/148‑1 (MQ=255)
                                                                                                                    cAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGttt                            >  1:2539359/1‑148 (MQ=255)
                                                                                                                    cAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGttt                            <  1:629355/148‑1 (MQ=255)
                                                                                                                      gggCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCt                          >  1:117171/1‑148 (MQ=255)
                                                                                                                      gggCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCt                          >  1:932972/1‑148 (MQ=255)
                                                                                                                       ggcggcGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTg                         >  2:1288334/1‑148 (MQ=255)
                                                                                                                         cggcgTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTa                       >  2:1167287/1‑148 (MQ=255)
                                                                                                                           gcgTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCCGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGAGAGCATCAGAATGTTTCTGTAcc                     >  1:3529959/1‑148 (MQ=255)
                                                                                                                            cgTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCa                    >  1:317763/1‑148 (MQ=255)
                                                                                                                            cgTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCa                    >  2:505737/1‑148 (MQ=255)
                                                                                                                            cgTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCa                    >  1:1200592/1‑148 (MQ=255)
                                                                                                                             gTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCaa                   <  1:2839678/148‑1 (MQ=255)
                                                                                                                             gTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCaa                   <  1:1857836/148‑1 (MQ=255)
                                                                                                                             gTGATTTTCGAAGCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCaa                   <  2:467165/148‑1 (MQ=255)
                                                                                                                              tGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAAt                  >  1:1914474/1‑148 (MQ=255)
                                                                                                                               gATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATg                 >  1:444513/1‑148 (MQ=255)
                                                                                                                                  tttCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCg              >  1:2034920/1‑148 (MQ=255)
                                                                                                                                  tttCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGAGGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCg              >  2:3781716/1‑148 (MQ=255)
                                                                                                                                   ttCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGc             <  2:1598690/148‑1 (MQ=255)
                                                                                                                                         gCGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTgg       >  2:1804716/1‑148 (MQ=255)
                                                                                                                                          cGAAACAGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGAGGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGa      >  2:3162672/1‑148 (MQ=255)
                                                                                                                                             aaCCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGAGCa   >  1:1722100/1‑148 (MQ=255)
                                                                                                                                              aCCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGAGCAg  >  2:2469494/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGAGCAG  >  CP009273/3517908‑3518197

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: