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,1160CA98.9% 530.1 / ‑2.8 171A523A (GCG→GCT) pbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (98/71);  total (100/71)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
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 171 total aligned reads displayed.

CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638972‑2639260
                                                                                                                                                |                                                                                                                                                
cAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGcc                                                                                                                                               <  1:67167/148‑1 (MQ=255)
cAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGcc                                                                                                                                               <  2:3443570/148‑1 (MQ=255)
cAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGcc                                                                                                                                               <  2:3443976/148‑1 (MQ=255)
 aGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGccc                                                                                                                                              <  1:10157/148‑1 (MQ=255)
 aGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGccc                                                                                                                                              <  1:10124/148‑1 (MQ=255)
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             >  2:2925070/1‑148 (MQ=255)
       gCACTGCGCGACAGTAAGATATTATTGAACTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCACTAGCCCAGGCGt                                                                                                                                        <  2:2596269/148‑1 (MQ=255)
        cACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTc                                                                                                                                       >  1:1203648/1‑148 (MQ=255)
        cACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTc                                                                                                                                       <  2:2320896/148‑1 (MQ=255)
        cACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTc                                                                                                                                       >  1:825513/1‑148 (MQ=255)
          cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     <  2:102416/148‑1 (MQ=255)
                gACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGc                                                                                                                               <  2:791335/148‑1 (MQ=255)
                gACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGc                                                                                                                               <  1:790710/148‑1 (MQ=255)
                 aCAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGcc                                                                                                                              >  1:1128171/1‑148 (MQ=255)
                   aGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGAGAGCCAt                                                                                                                            <  1:2265040/148‑1 (MQ=255)
                    gTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATa                                                                                                                           >  1:1593971/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          >  2:2520181/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          <  1:1871968/148‑1 (MQ=255)
                      aaGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGc                                                                                                                         >  2:1281365/1‑148 (MQ=255)
                        gATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTg                                                                                                                       >  1:922807/1‑148 (MQ=255)
                          tattatTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGt                                                                                                                     >  2:448215/1‑148 (MQ=255)
                           attattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTg                                                                                                                    >  1:1253828/1‑148 (MQ=255)
                           attattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTg                                                                                                                    >  1:1253961/1‑148 (MQ=255)
                            ttattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGc                                                                                                                   >  1:1551133/1‑148 (MQ=255)
                            ttattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGc                                                                                                                   <  1:651869/148‑1 (MQ=255)
                              attGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGAt                                                                                                                                  >  2:395554/1‑131 (MQ=255)
                              attGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGAt                                                                                                                                  <  1:395554/131‑1 (MQ=255)
                              attGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGccc                                                                                                                 >  1:147140/1‑148 (MQ=255)
                               ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                <  2:963499/148‑1 (MQ=255)
                                tGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGt                                                                                                               <  1:2306141/148‑1 (MQ=255)
                                tGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGt                                                                                                               <  2:2449322/148‑1 (MQ=255)
                                tGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGt                                                                                                               <  2:1619454/148‑1 (MQ=255)
                                tGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGt                                                                                                               >  1:1670348/1‑148 (MQ=255)
                                tGACCAGATTCAAAAATGGTACGGCACTGGCAAAGCCAAACTGACCCACAACGGGCGTGCCGACCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGt                                                                                                               >  2:702684/1‑148 (MQ=255)
                                 gACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGtt                                                                                                              >  2:2761300/1‑148 (MQ=255)
                                 gACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGtt                                                                                                              >  2:2096472/1‑148 (MQ=255)
                                   ccTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGtttt                                                                                                            <  1:1788153/148‑1 (MQ=255)
                                     tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          >  1:3694740/1‑148 (MQ=255)
                                     tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          >  2:2809630/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  1:2636746/1‑148 (MQ=255)
                                           aacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCa                                                                                                    >  1:2994105/1‑148 (MQ=255)
                                            acaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAg                                                                                                   >  2:1332220/1‑148 (MQ=255)
                                             caaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGt                                                                                                  >  1:3469676/1‑148 (MQ=255)
                                              aaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTg                                                                                                 <  2:3390818/148‑1 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                >  2:2795169/1‑148 (MQ=255)
                                                tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:2213263/1‑148 (MQ=255)
                                                 ggTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgggg                                                                                              >  1:3346894/1‑148 (MQ=255)
                                                      ggCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcg                                                                                         >  2:1654636/1‑148 (MQ=255)
                                                      ggCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcg                                                                                         >  2:3523897/1‑148 (MQ=255)
                                                      ggCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcg                                                                                         >  2:3549360/1‑148 (MQ=255)
                                                           tGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCgg                                                                                    <  1:3533331/148‑1 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 >  2:2792074/1‑148 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 <  2:3783938/148‑1 (MQ=255)
                                                                  gCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCa                                                                             <  1:2703218/148‑1 (MQ=255)
                                                                           aCCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGc                                                                    >  1:57816/1‑148 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                >  1:3178739/1‑148 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                <  1:2878070/148‑1 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                <  1:802286/148‑1 (MQ=255)
                                                                                caacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGGACTATCCGGCAAggg                                                               <  2:3237283/148‑1 (MQ=255)
                                                                                     ggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg                                                          >  2:1063111/1‑148 (MQ=255)
                                                                                      gCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcgc                                                         <  1:3014505/148‑1 (MQ=255)
                                                                                       cGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCt                                                        >  1:1539400/1‑148 (MQ=255)
                                                                                         tGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtc                                                      >  1:16686/1‑148 (MQ=255)
                                                                                          gCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtca                                                     >  2:1057962/1‑148 (MQ=255)
                                                                                              tCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAg                                                 >  1:2563276/1‑148 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  1:998048/1‑148 (MQ=255)
                                                                                                 ggTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCa                                              <  2:117027/148‑1 (MQ=255)
                                                                                                 ggTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGCGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCa                                              <  2:394419/148‑1 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                           >  1:3120768/1‑148 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                           <  2:3011948/148‑1 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGACCAGGCGTCACGATAGCAATAGCTTGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                           >  1:3120770/1‑148 (MQ=255)
                                                                                                     tGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAAt                                          <  1:640018/148‑1 (MQ=255)
                                                                                                       ccAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATcc                                        >  1:1278517/1‑148 (MQ=255)
                                                                                                            ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   >  1:2770209/1‑148 (MQ=255)
                                                                                                               gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                >  1:2304983/1‑148 (MQ=255)
                                                                                                                aTCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGa                               >  1:3777311/1‑148 (MQ=255)
                                                                                                                  cccAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATc                             <  2:1804523/148‑1 (MQ=255)
                                                                                                                  cccAATGACAGACCCAGCCTTAACCCCCAGAGCCCAGGCGTCACGAGAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGTCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATc                             <  2:3695863/148‑1 (MQ=255)
                                                                                                                     aaTGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAc                          >  2:1485588/1‑148 (MQ=255)
                                                                                                                      aTGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcg                         >  1:776077/1‑148 (MQ=255)
                                                                                                                         aCATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgccg                      >  2:2327528/1‑148 (MQ=255)
                                                                                                                         aCATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgccg                      >  1:1975195/1‑148 (MQ=255)
                                                                                                                            tAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccc                   >  1:1852291/1‑148 (MQ=255)
                                                                                                                            tAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccc                   >  2:665344/1‑148 (MQ=255)
                                                                                                                             agcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                  >  1:102337/1‑148 (MQ=255)
                                                                                                                              gcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccc                 >  1:2126433/1‑148 (MQ=255)
                                                                                                                              gcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccc                 >  1:2126589/1‑148 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  2:2184080/148‑1 (MQ=255)
                                                                                                                                   cgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacga            >  2:193558/1‑148 (MQ=255)
                                                                                                                                   cgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacga            >  2:193581/1‑148 (MQ=255)
                                                                                                                                   cgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacga            >  1:1371852/1‑148 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCTTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATTATCCACGCCGCCCCCgacgac           <  1:702684/148‑1 (MQ=255)
                                                                                                                                     ttAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACa          >  1:2857454/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3059334/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         <  2:2734411/148‑1 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:1152402/1‑148 (MQ=255)
                                                                                                                                         ccccAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTaa      <  1:1749047/148‑1 (MQ=255)
                                                                                                                                         ccccAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTaa      <  1:536155/148‑1 (MQ=255)
                                                                                                                                          cccAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTaaa     <  2:1582562/148‑1 (MQ=255)
                                                                                                                                             aaTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGc  >  1:3730390/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: ‑