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,6310CA96.1% 1104.2 / 29.4 332R348L (CGC→CTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base A (149/170);  total (155/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 337 total aligned reads displayed.

TAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACC  >  CP009273/662484‑662776
                                                                                                                                                   |                                                                                                                                                 
 aGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGagcg                                                                                                                                                >  2:183769/1‑150 (MQ=255)
   gCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTg                                                                                                                                              >  2:1566022/1‑150 (MQ=255)
   gCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTg                                                                                                                                              >  1:1514186/1‑150 (MQ=255)
   gCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTg                                                                                                                                              >  2:1978184/1‑150 (MQ=255)
   gCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTg                                                                                                                                              <  1:2006575/150‑1 (MQ=255)
   gCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTACCAGTCACGATAACGTTTTTCCGAACCTGGAAGTTTCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTg                                                                                                                                              >  2:4531900/1‑150 (MQ=255)
    ccACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGa                                                                                                                                             >  2:5326/1‑150 (MQ=255)
    ccACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCAATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACAAGCCTGGGTCAAACATCGTCGTATTTAGCGTGa                                                                                                                                             >  1:3352139/1‑150 (MQ=255)
     cACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGAt                                                                                                                                            >  1:4031813/1‑150 (MQ=255)
       ccTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCa                                                                                                                                          >  1:1118032/1‑150 (MQ=255)
       ccTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCa                                                                                                                                          >  1:7846/1‑150 (MQ=255)
        cTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAc                                                                                                                                         >  1:4591302/1‑150 (MQ=255)
            tAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCg                                                                                                                                     <  1:4189044/150‑1 (MQ=255)
            tAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCg                                                                                                                                     <  1:294088/150‑1 (MQ=255)
             agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCgg                                                                                                                                    >  2:4448588/1‑150 (MQ=255)
             agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCgg                                                                                                                                    >  2:110930/1‑150 (MQ=255)
               aagaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgcg                                                                                                                                  >  1:1214532/1‑150 (MQ=255)
                 gaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCt                                                                                                                                <  1:1566489/150‑1 (MQ=255)
                  aagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTc                                                                                                                               >  1:4419308/1‑150 (MQ=255)
                  aagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTc                                                                                                                               <  2:2808534/150‑1 (MQ=255)
                     gTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAAt                                                                                                                            >  1:3165543/1‑150 (MQ=255)
                        tCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGcc                                                                                                                         <  2:3889776/150‑1 (MQ=255)
                             aGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAAc                                                                                                                    <  1:2619488/150‑1 (MQ=255)
                               aTTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCg                                                                                                                  >  1:1466427/1‑150 (MQ=255)
                               aTTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCg                                                                                                                  <  1:2680385/150‑1 (MQ=255)
                                  cttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCca                                                                                                               <  2:4561558/150‑1 (MQ=255)
                                  cttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCca                                                                                                               >  1:345811/1‑150 (MQ=255)
                                   ttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCcac                                                                                                              >  1:4178759/1‑150 (MQ=255)
                                    tcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCcaca                                                                                                             >  1:4580086/1‑150 (MQ=255)
                                       aGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAg                                                                                                          <  1:314660/150‑1 (MQ=255)
                                         cGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAggg                                                                                                        <  2:614199/150‑1 (MQ=255)
                                          gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt                                                                                                       >  1:708228/1‑150 (MQ=255)
                                           aTCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                                      <  2:2019664/150‑1 (MQ=255)
                                           aTCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                                      <  1:4489751/150‑1 (MQ=255)
                                              ttGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTaa                                                                                                   <  1:2164038/150‑1 (MQ=255)
                                                  gACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGt                                                                                               <  1:3580336/150‑1 (MQ=255)
                                                  gACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGt                                                                                               <  1:679367/150‑1 (MQ=255)
                                                    cATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAg                                                                                             <  1:3998060/150‑1 (MQ=255)
                                                      ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                                           >  1:88978/1‑150 (MQ=255)
                                                      ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                                           <  1:5107606/150‑1 (MQ=255)
                                                         aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCACTGCCGAAACCGCCACAGAGGGTTTAACTGTACAcgcg                                                                                        >  1:4289463/1‑150 (MQ=255)
                                                         aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGGTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcgcg                                                                                        <  1:1165904/150‑1 (MQ=255)
                                                              cccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagg                                                                                   <  2:3212921/150‑1 (MQ=255)
                                                                 gTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATa                                                                                >  2:713591/1‑150 (MQ=255)
                                                                   ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaaa                                                                              >  2:2911562/1‑150 (MQ=255)
                                                                   ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaaa                                                                              >  1:2459060/1‑150 (MQ=255)
                                                                       ccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAcccc                                                                          >  1:4052170/1‑150 (MQ=255)
                                                                        cATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAccccc                                                                         <  2:3075354/150‑1 (MQ=255)
                                                                           tttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgt                                                                      <  2:1614806/150‑1 (MQ=255)
                                                                            ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                                     <  2:871511/150‑1 (MQ=255)
                                                                            ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                                     <  2:2379372/150‑1 (MQ=255)
                                                                                ccAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGgc                                                                 <  2:1214532/150‑1 (MQ=255)
                                                                                    tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCgg                                                             >  1:4103427/1‑150 (MQ=255)
                                                                                    tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCgg                                                             <  1:3087608/150‑1 (MQ=255)
                                                                                    tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAGAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCgg                                                             >  2:704948/1‑150 (MQ=255)
                                                                                     cACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGt                                                            <  1:3328224/150‑1 (MQ=255)
                                                                                      aCGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGtt                                                           >  2:2760044/1‑150 (MQ=255)
                                                                                      aCGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGtt                                                           <  1:209066/150‑1 (MQ=255)
                                                                                        gATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCa                                                         <  1:2764252/150‑1 (MQ=255)
                                                                                          tAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAcc                                                       >  1:4547572/1‑150 (MQ=255)
                                                                                           aaCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCa                                                      >  2:67886/1‑150 (MQ=255)
                                                                                           aaCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCa                                                      <  2:1251467/150‑1 (MQ=255)
                                                                                             cGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGc                                                    >  2:4695319/1‑150 (MQ=255)
                                                                                              gTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCg                                                   >  2:2497204/1‑150 (MQ=255)
                                                                                                ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                                 >  1:3614978/1‑150 (MQ=255)
                                                                                                 tttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtg                                                <  1:1299028/150‑1 (MQ=255)
                                                                                                      gAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                           <  1:2383275/150‑1 (MQ=255)
                                                                                                      gAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                           <  1:1824709/150‑1 (MQ=255)
                                                                                                       aaCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCg                                          <  1:1719236/150‑1 (MQ=255)
                                                                                                       aaCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCg                                          >  2:1818253/1‑150 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      >  2:2185093/1‑150 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      >  2:4978310/1‑150 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      >  2:4348010/1‑150 (MQ=255)
                                                                                                             gCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGt                                    >  2:4628223/1‑150 (MQ=255)
                                                                                                               aGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                                  <  1:917794/150‑1 (MQ=255)
                                                                                                                 ttGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaa                                <  2:4802011/150‑1 (MQ=255)
                                                                                                                     caccaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAg                            >  1:3801596/1‑150 (MQ=255)
                                                                                                                       ccaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGc                          <  1:4701887/150‑1 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                         >  2:34301/1‑150 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                         <  2:4902246/150‑1 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                         >  2:3762342/1‑150 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                         >  2:1881197/1‑150 (MQ=255)
                                                                                                                         aGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCgg                        <  1:1984953/150‑1 (MQ=255)
                                                                                                                           ccTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaa                      <  1:5085405/150‑1 (MQ=255)
                                                                                                                              gggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGactaa                   >  1:3883305/1‑150 (MQ=255)
                                                                                                                              gggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataa                   <  1:4780789/150‑1 (MQ=255)
                                                                                                                                gTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATc                 >  1:563840/1‑150 (MQ=255)
                                                                                                                                 tCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGATACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCt                <  2:4097971/150‑1 (MQ=255)
                                                                                                                                   aaaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCttt              <  2:3352139/150‑1 (MQ=255)
                                                                                                                                     aCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGc            >  2:1716814/1‑150 (MQ=255)
                                                                                                                                      cAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCt           <  1:2030914/150‑1 (MQ=255)
                                                                                                                                      cAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCt           <  2:2611893/150‑1 (MQ=255)
                                                                                                                                       aGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTg          <  2:2587734/150‑1 (MQ=255)
                                                                                                                                       aGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTg          <  2:1451660/150‑1 (MQ=255)
                                                                                                                                          gtcgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGgag       >  2:795045/1‑150 (MQ=255)
                                                                                                                                              tATTGAGCGTGATCACCCCGGCGCTCAATGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAc   <  2:3566687/150‑1 (MQ=255)
                                                                                                                                              tATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAc   <  2:2876011/150‑1 (MQ=255)
                                                                                                                                              tATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAc   <  1:4507151/150‑1 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc  >  1:450211/1‑150 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc  <  2:3113829/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACC  >  CP009273/662484‑662776

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: