Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 662,631 C→A 100% R348L (CGC→CTC)  mrdA ← transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273662,6310CA95.9% 1113.0 / 30.4 319R348L (CGC→CTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base C (4/9);  new base A (138/168);  total (142/177)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 323 total aligned reads displayed.

TAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACCG  >  CP009273/662484‑662777
                                                                                                                                                   |                                                                                                                                                  
tAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTgagc                                                                                                                                                  <  2:849682/150‑1 (MQ=255)
  ggCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGt                                                                                                                                                >  2:3646750/1‑150 (MQ=255)
       ccTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCa                                                                                                                                           <  2:4040568/150‑1 (MQ=255)
       ccTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCa                                                                                                                                           <  2:4466981/150‑1 (MQ=255)
           aTAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAcccc                                                                                                                                       <  1:417457/150‑1 (MQ=255)
           aTAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAcccc                                                                                                                                       <  2:4724236/150‑1 (MQ=255)
             agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCgg                                                                                                                                     >  2:1021203/1‑150 (MQ=255)
               aagaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgcg                                                                                                                                   >  1:4670005/1‑150 (MQ=255)
               aagaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgcg                                                                                                                                   <  1:1549855/150‑1 (MQ=255)
                 gaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCt                                                                                                                                 <  1:1044668/150‑1 (MQ=255)
                   agGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACTATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGAGTCGTATAGAGAGTGATCACCCCGGCGCTCa                                                                                                                               >  1:2147655/1‑150 (MQ=255)
                   agGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCa                                                                                                                               >  1:101575/1‑150 (MQ=255)
                      tATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATg                                                                                                                            >  1:1786668/1‑150 (MQ=255)
                        tCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGcc                                                                                                                          <  1:4491930/150‑1 (MQ=255)
                            cAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGaaa                                                                                                                      <  2:6880/150‑1 (MQ=255)
                            cAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGaaa                                                                                                                      <  2:4111574/150‑1 (MQ=255)
                                ttcttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGc                                                                                                                  >  2:4302310/1‑150 (MQ=255)
                                    tcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCcaca                                                                                                              >  2:2607618/1‑150 (MQ=255)
                                    tcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCcaca                                                                                                              <  1:1712156/150‑1 (MQ=255)
                                      cAGCGAGCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTAGTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACCTa                                                                                                            >  2:4623433/1‑150 (MQ=255)
                                       aGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAg                                                                                                           >  2:1994913/1‑150 (MQ=255)
                                          gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt                                                                                                        >  1:271736/1‑150 (MQ=255)
                                            tCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGttt                                                                                                      <  2:1665949/150‑1 (MQ=255)
                                             cTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATTACCCCGGCGCTCAATGCCGAAACCGCCACATCGGGTTTa                                                                                                     >  2:3100870/1‑150 (MQ=255)
                                               tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAAc                                                                                                   >  1:1334851/1‑150 (MQ=255)
                                                     aTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGa                                                                                             <  1:3652859/150‑1 (MQ=255)
                                                        cAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcgc                                                                                          >  1:685796/1‑150 (MQ=255)
                                                             gCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggag                                                                                     <  2:2305944/150‑1 (MQ=255)
                                                             gCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggag                                                                                     <  1:2059883/150‑1 (MQ=255)
                                                               ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagga                                                                                   <  2:2172173/150‑1 (MQ=255)
                                                                         aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                                         >  1:3378332/1‑150 (MQ=255)
                                                                         aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                                         <  1:1417672/150‑1 (MQ=255)
                                                                          ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                                        <  1:3449421/150‑1 (MQ=255)
                                                                          ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                                        <  1:911647/150‑1 (MQ=255)
                                                                          ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                                        <  2:2848944/150‑1 (MQ=255)
                                                                          ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                                        <  2:1518631/150‑1 (MQ=255)
                                                                            ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                                      >  1:1866859/1‑150 (MQ=255)
                                                                            ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                                      >  2:1896203/1‑150 (MQ=255)
                                                                            ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                                      >  2:2502925/1‑150 (MQ=255)
                                                                             tttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtgt                                                                     >  1:4241279/1‑150 (MQ=255)
                                                                             tttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtgt                                                                     <  2:4496605/150‑1 (MQ=255)
                                                                              ttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtgtg                                                                    <  1:757001/150‑1 (MQ=255)
                                                                                   gTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGcgcg                                                               <  1:3071710/150‑1 (MQ=255)
                                                                                   gTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGcgcg                                                               <  2:2043782/150‑1 (MQ=255)
                                                                                   gTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGcgcg                                                               <  2:4375400/150‑1 (MQ=255)
                                                                                    tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCgg                                                              >  1:3562338/1‑150 (MQ=255)
                                                                                      aCGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGtt                                                            <  1:2963844/150‑1 (MQ=255)
                                                                                          tAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAcc                                                        <  1:2396227/150‑1 (MQ=255)
                                                                                            aCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAg                                                      >  2:2710813/1‑150 (MQ=255)
                                                                                            aCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAg                                                      <  2:389382/150‑1 (MQ=255)
                                                                                              gTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCg                                                    <  1:4555384/150‑1 (MQ=255)
                                                                                                ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                                  <  1:95506/150‑1 (MQ=255)
                                                                                                ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                                  <  1:2468647/150‑1 (MQ=255)
                                                                                                 tttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtg                                                 <  2:2997017/150‑1 (MQ=255)
                                                                                                 tttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtg                                                 <  1:369839/150‑1 (MQ=255)
                                                                                                  ttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                                <  2:121040/150‑1 (MQ=255)
                                                                                                  ttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                                >  1:3239586/1‑150 (MQ=255)
                                                                                                  ttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGGTCACCAGCGgtgt                                                >  1:1713464/1‑150 (MQ=255)
                                                                                                    ccGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAt                                              <  2:3291862/150‑1 (MQ=255)
                                                                                                      gAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                            <  2:3482536/150‑1 (MQ=255)
                                                                                                        aCCTGGCAGGTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                                          <  1:1182891/150‑1 (MQ=255)
                                                                                                         ccTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                                         >  2:4219974/1‑150 (MQ=255)
                                                                                                         ccTGGCAGGTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                                         <  2:1263867/150‑1 (MQ=255)
                                                                                                          cTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                                        <  2:2777481/150‑1 (MQ=255)
                                                                                                          cTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                                        <  2:91390/150‑1 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                       <  2:3815207/150‑1 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                       <  1:854531/150‑1 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                       >  1:1386605/1‑150 (MQ=255)
                                                                                                             gCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGt                                     <  2:3323133/150‑1 (MQ=255)
                                                                                                              cAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGtt                                    <  1:3920150/150‑1 (MQ=255)
                                                                                                              cAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGtt                                    <  2:1786668/150‑1 (MQ=255)
                                                                                                               aGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                                   <  2:2653155/150‑1 (MQ=255)
                                                                                                               aGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                                   >  1:267769/1‑150 (MQ=255)
                                                                                                                    ccaccaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAGACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaacaa                              >  2:1687229/1‑150 (MQ=255)
                                                                                                                     caccaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAg                             >  1:1115544/1‑150 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                          >  1:1123075/1‑150 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                          >  1:4331511/1‑150 (MQ=255)
                                                                                                                           ccTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaa                       <  2:621987/150‑1 (MQ=255)
                                                                                                                           ccTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaa                       <  1:1493291/150‑1 (MQ=255)
                                                                                                                            cTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaat                      <  1:1636742/150‑1 (MQ=255)
                                                                                                                               ggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat                   <  1:164622/150‑1 (MQ=255)
                                                                                                                                gTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATc                  <  2:520904/150‑1 (MQ=255)
                                                                                                                                  cAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATACTCtt                >  1:4152444/1‑150 (MQ=255)
                                                                                                                                  cAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCtt                <  1:4018160/150‑1 (MQ=255)
                                                                                                                                  cAAACAGCGTCGTATTGAGAGTGATCACCCCGGCGCTCAATTCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGTCGCGGTTCACCAGAGGTGTATTCGGATCGTTCAACAAGGCGTAATAATCtt                >  2:1144718/1‑150 (MQ=255)
                                                                                                                                   aaaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCttt               <  2:2131836/150‑1 (MQ=255)
                                                                                                                                    aaCAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTg              <  2:3761404/150‑1 (MQ=255)
                                                                                                                                    aaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTg              >  1:486191/1‑150 (MQ=255)
                                                                                                                                      cAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCt            <  2:2223636/150‑1 (MQ=255)
                                                                                                                                      cAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCt            >  2:42396/1‑150 (MQ=255)
                                                                                                                                      cAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCt            <  2:459993/150‑1 (MQ=255)
                                                                                                                                       aGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTg           >  1:2232560/1‑150 (MQ=255)
                                                                                                                                         cgtcgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGga         <  2:2147209/150‑1 (MQ=255)
                                                                                                                                         cgtcgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGga         <  2:4298995/150‑1 (MQ=255)
                                                                                                                                             gtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATa     <  2:1442815/150‑1 (MQ=255)
                                                                                                                                              tATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAc    >  2:1008201/1‑150 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGTCGCTCAATTCCGAAACCGCCTAATAGTGTTTAACTGTAGACGCGGGAGTATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCTGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc   <  1:2021188/150‑1 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc   >  1:443062/1‑150 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc   >  1:634718/1‑150 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc   >  1:257052/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACCG  >  CP009273/662484‑662777

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: