Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,639,116 C→A 100% A523A (GCG→GCT pbpC ← penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,639,1160CA99.3% 559.2 / ‑3.9 186A523A (GCG→GCTpbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (102/81);  total (103/83)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 186 total aligned reads displayed.

CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638972‑2639260
                                                                                                                                                |                                                                                                                                                
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             >  1:133911/1‑148 (MQ=255)
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             >  1:134006/1‑148 (MQ=255)
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             <  2:3069665/148‑1 (MQ=255)
   aTTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAg                                                                                                                                            >  2:27283/1‑148 (MQ=255)
         aCTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCa                                                                                                                                      <  1:2443979/148‑1 (MQ=255)
          cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     >  2:1233922/1‑148 (MQ=255)
               cgACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAg                                                                                                                                >  2:3448660/1‑148 (MQ=255)
               cgACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAg                                                                                                                                <  1:1771281/148‑1 (MQ=255)
                gACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGc                                                                                                                               <  2:1567612/148‑1 (MQ=255)
                  cAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCa                                                                                                                             >  1:3154309/1‑148 (MQ=255)
                   aGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCAt                                                                                                                            <  2:3434214/148‑1 (MQ=255)
                    gTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATa                                                                                                                           >  2:241720/1‑148 (MQ=255)
                    gTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATa                                                                                                                           >  1:719005/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          >  1:2266563/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          >  1:594757/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          <  2:1634032/148‑1 (MQ=255)
                      aaGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGc                                                                                                                         >  1:2877241/1‑148 (MQ=255)
                      aaGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGc                                                                                                                         >  1:435979/1‑148 (MQ=255)
                       aGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCt                                                                                                                        <  2:2574279/148‑1 (MQ=255)
                              attGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGccc                                                                                                                 >  1:3210624/1‑148 (MQ=255)
                                    cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                           >  1:1268911/1‑148 (MQ=255)
                                     tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          <  1:1990148/148‑1 (MQ=255)
                                       aTTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAt                                                                                                        >  1:1733534/1‑148 (MQ=255)
                                        ttCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATg                                                                                                       >  1:3366604/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  1:3500805/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  1:944231/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  2:2647305/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  2:2747501/1‑148 (MQ=255)
                                           aacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCa                                                                                                    <  1:243889/148‑1 (MQ=255)
                                            acaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAg                                                                                                   >  2:1492599/1‑148 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                >  1:455609/1‑148 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                <  2:349823/148‑1 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                <  2:350668/148‑1 (MQ=255)
                                                tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:2877103/1‑148 (MQ=255)
                                                tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:827907/1‑148 (MQ=255)
                                                  gTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGc                                                                                             >  2:887954/1‑148 (MQ=255)
                                                     cGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAc                                                                                          >  2:1534187/1‑148 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       >  1:107597/1‑148 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       <  2:3015576/148‑1 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       >  1:1883035/1‑148 (MQ=255)
                                                          cTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcgcg                                                                                     >  2:118136/1‑148 (MQ=255)
                                                             gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                                                                  <  1:525842/148‑1 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 <  2:199920/148‑1 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 >  2:1267901/1‑148 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 <  2:594757/148‑1 (MQ=255)
                                                                   ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                            >  1:1057940/1‑148 (MQ=255)
                                                                    cAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACt                                                                           <  1:871074/148‑1 (MQ=255)
                                                                      aaCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAt                                                                         <  2:3238683/148‑1 (MQ=255)
                                                                      aaCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAt                                                                         <  2:197898/148‑1 (MQ=255)
                                                                       aCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATc                                                                        <  2:1604109/148‑1 (MQ=255)
                                                                            ccAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCa                                                                   <  1:1582911/148‑1 (MQ=255)
                                                                              aacaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAg                                                                 <  2:2083389/148‑1 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                >  1:3180192/1‑148 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                >  1:2241794/1‑148 (MQ=255)
                                                                                caacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAggg                                                               <  1:2610463/148‑1 (MQ=255)
                                                                                   cGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg                                                            <  2:2226015/148‑1 (MQ=255)
                                                                                   cGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg                                                            <  2:2226084/148‑1 (MQ=255)
                                                                                   cGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg                                                            >  1:1502282/1‑148 (MQ=255)
                                                                                    gggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgc                                                           <  2:2444719/148‑1 (MQ=255)
                                                                                     ggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg                                                          >  2:3140561/1‑148 (MQ=255)
                                                                                     ggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg                                                          <  2:2951790/148‑1 (MQ=255)
                                                                                      gCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcgc                                                         >  1:194520/1‑148 (MQ=255)
                                                                                          gCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtca                                                     >  1:1758033/1‑148 (MQ=255)
                                                                                              tCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAg                                                 >  1:1969699/1‑148 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  2:942818/1‑148 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  1:803622/1‑148 (MQ=255)
                                                                                                cGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGcc                                               <  2:1268911/148‑1 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                           <  1:117951/148‑1 (MQ=255)
                                                                                                      gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGGCTCTGGGGCGTGGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         <  1:3448218/148‑1 (MQ=255)
                                                                                                      gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         >  2:756664/1‑148 (MQ=255)
                                                                                                      gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         <  1:3490171/148‑1 (MQ=255)
                                                                                                        cAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCACTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCg                                       <  2:2528599/148‑1 (MQ=255)
                                                                                                         aGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCgg                                      <  1:2472897/148‑1 (MQ=255)
                                                                                                          gTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGc                                     >  1:1669540/1‑148 (MQ=255)
                                                                                                          gTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGc                                     >  1:209162/1‑148 (MQ=255)
                                                                                                            ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   <  2:1909918/148‑1 (MQ=255)
                                                                                                            ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   >  1:2742956/1‑148 (MQ=255)
                                                                                                            ccAGATCCCAATGACAGAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   <  1:457837/148‑1 (MQ=255)
                                                                                                                aTCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGa                               >  1:573180/1‑148 (MQ=255)
                                                                                                                 tCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGAt                              <  1:3030995/148‑1 (MQ=255)
                                                                                                                   ccAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATcc                            <  2:3118059/148‑1 (MQ=255)
                                                                                                                    cAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCa                           >  1:233501/1‑148 (MQ=255)
                                                                                                                      aTGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcg                         <  1:27283/148‑1 (MQ=255)
                                                                                                                      aTGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcg                         <  1:442074/148‑1 (MQ=255)
                                                                                                                        gACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgcc                       >  1:49602/1‑148 (MQ=255)
                                                                                                                         aCATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgccg                      >  1:2908180/1‑148 (MQ=255)
                                                                                                                             agcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                  >  2:2516036/1‑148 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  1:3238553/148‑1 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  1:2142740/148‑1 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  1:3238557/148‑1 (MQ=255)
                                                                                                                                  gcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacg             <  1:999005/148‑1 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCgacgac           <  2:134006/148‑1 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCgacgac           <  2:435979/148‑1 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCgacgac           <  2:133911/148‑1 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3177138/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3483741/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3037442/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:1747324/1‑148 (MQ=255)
                                                                                                                                            cAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg   <  2:719005/148‑1 (MQ=255)
                                                                                                                                            cAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg   >  1:2837098/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638972‑2639260

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: