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,6310CA99.7% 1046.3 / ‑4.5 317R348L (CGC→CTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (160/156);  total (160/157)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-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 319 total aligned reads displayed.

AGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACC  >  CP009273/662485‑662776
                                                                                                                                                  |                                                                                                                                                 
  gCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTg                                                                                                                                              <  1:3080303/150‑1 (MQ=255)
    cACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCGGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGAt                                                                                                                                            >  1:3264279/1‑150 (MQ=255)
       cTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAc                                                                                                                                         <  1:29026/150‑1 (MQ=255)
         gATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAccc                                                                                                                                       >  2:3736125/1‑150 (MQ=255)
         gATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAccc                                                                                                                                       >  1:2708036/1‑150 (MQ=255)
         gATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAccc                                                                                                                                       <  2:2685784/150‑1 (MQ=255)
          aTAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAcccc                                                                                                                                      >  1:2563226/1‑150 (MQ=255)
          aTAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAcccc                                                                                                                                      >  2:4686930/1‑150 (MQ=255)
           tAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCg                                                                                                                                     <  2:327255/150‑1 (MQ=255)
            agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCgg                                                                                                                                    >  2:655510/1‑150 (MQ=255)
            agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCgg                                                                                                                                    >  2:137047/1‑150 (MQ=255)
             gaagaaGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgc                                                                                                                                   <  2:1454111/150‑1 (MQ=255)
               agaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgcgc                                                                                                                                 >  1:724659/1‑150 (MQ=255)
                  agGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCa                                                                                                                              >  2:4032542/1‑150 (MQ=255)
                      aTCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGc                                                                                                                          >  2:4541381/1‑150 (MQ=255)
                      aTCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGc                                                                                                                          >  1:1881599/1‑150 (MQ=255)
                        ccGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCg                                                                                                                        >  2:2877207/1‑150 (MQ=255)
                        ccGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCg                                                                                                                        <  2:199804/150‑1 (MQ=255)
                         cGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGa                                                                                                                       >  2:2198809/1‑150 (MQ=255)
                          gCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGaa                                                                                                                      >  2:5004100/1‑150 (MQ=255)
                           cAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGaaa                                                                                                                     >  1:1490697/1‑150 (MQ=255)
                            aGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAAc                                                                                                                    >  2:2243359/1‑150 (MQ=255)
                            aGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAAc                                                                                                                    >  2:578392/1‑150 (MQ=255)
                             gATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAAcc                                                                                                                   <  1:759864/150‑1 (MQ=255)
                                    ccAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt                                                                                                            >  1:3923818/1‑150 (MQ=255)
                                    ccAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt                                                                                                            <  2:490898/150‑1 (MQ=255)
                                     cAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATa                                                                                                           <  2:1260249/150‑1 (MQ=255)
                                         gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt                                                                                                       >  1:2034117/1‑150 (MQ=255)
                                         gATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt                                                                                                       >  2:4173577/1‑150 (MQ=255)
                                          aTCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGATTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGtt                                                                                                      <  2:4185875/150‑1 (MQ=255)
                                           tCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGttt                                                                                                     <  1:3007065/150‑1 (MQ=255)
                                            cTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTa                                                                                                    <  1:2691095/150‑1 (MQ=255)
                                            cTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTa                                                                                                    >  2:3296700/1‑150 (MQ=255)
                                              tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAAc                                                                                                  >  1:941844/1‑150 (MQ=255)
                                               gtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACt                                                                                                 <  1:2727787/150‑1 (MQ=255)
                                                tgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTg                                                                                                >  1:3809161/1‑150 (MQ=255)
                                                   cATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAg                                                                                             <  2:4295437/150‑1 (MQ=255)
                                                     ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAc                                                                                           <  2:1507149/150‑1 (MQ=255)
                                                        aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcgcg                                                                                        >  2:3773076/1‑150 (MQ=255)
                                                        aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcgcg                                                                                        <  1:1553396/150‑1 (MQ=255)
                                                           cGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGgga                                                                                     >  2:4817984/1‑150 (MQ=255)
                                                            gCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggag                                                                                    <  2:1435702/150‑1 (MQ=255)
                                                              ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagga                                                                                  >  2:3798088/1‑150 (MQ=255)
                                                              ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggagga                                                                                  >  2:3507111/1‑150 (MQ=255)
                                                               cGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGAt                                                                                 <  1:1962764/150‑1 (MQ=255)
                                                                 tGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaa                                                                               <  1:4942690/150‑1 (MQ=255)
                                                                 tGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaa                                                                               >  2:2601869/1‑150 (MQ=255)
                                                                      ccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAcccc                                                                          >  1:921144/1‑150 (MQ=255)
                                                                        aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                                        <  1:603978/150‑1 (MQ=255)
                                                                         ttttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtg                                                                       >  1:4091917/1‑150 (MQ=255)
                                                                                 aGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGgcgc                                                               >  1:4870294/1‑150 (MQ=255)
                                                                                       gagaaCGTTTTTCCGAACCTGGCCGTTGCCACCCGCCTGGGTCACACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCACTGCCGCCCCCGCCACATCGGGTTTAACTGGAGACGCGGGAGGATAAACCCCCTGGGTGGCGCGGTTCa                                                         <  2:1053724/147‑1 (MQ=255)
                                                                                       gATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCa                                                         >  1:3597774/1‑150 (MQ=255)
                                                                                        aTAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAGACCCCCTGTGTGGCGCGGTTCAc                                                        >  2:3823774/1‑150 (MQ=255)
                                                                                        aTAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAc                                                        <  1:4423346/150‑1 (MQ=255)
                                                                                        aTAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAc                                                        >  2:2660670/1‑150 (MQ=255)
                                                                                          aaCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCa                                                      <  1:2130904/150‑1 (MQ=255)
                                                                                          aaCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCa                                                      <  1:4781942/150‑1 (MQ=255)
                                                                                          aaCGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCa                                                      >  1:1067404/1‑150 (MQ=255)
                                                                                            cGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGc                                                    >  2:1576791/1‑150 (MQ=255)
                                                                                            cGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGc                                                    >  2:3383525/1‑150 (MQ=255)
                                                                                               ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                                 <  1:3689885/150‑1 (MQ=255)
                                                                                               ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCACCACATAGGGTTTAACTGTAGACGCGTGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgt                                                 >  1:1355525/1‑150 (MQ=255)
                                                                                                 ttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                               >  1:1763324/1‑150 (MQ=255)
                                                                                                  tCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTa                                              <  2:3146081/150‑1 (MQ=255)
                                                                                                   ccGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAt                                             >  1:3284487/1‑150 (MQ=255)
                                                                                                    cGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAtt                                            <  2:1876663/150‑1 (MQ=255)
                                                                                                    cGAACCTGGCAGTTGCCACCACCCTGGGTCAAACAGCGTCGTCTTGAGCGTGATCACCCCGGCGCTCAATTCAGAAACCGCCACATCGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAtt                                            <  1:3477684/150‑1 (MQ=255)
                                                                                                     gAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                           <  1:167734/150‑1 (MQ=255)
                                                                                                     gAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                           <  2:2299615/150‑1 (MQ=255)
                                                                                                     gAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                           >  2:2532549/1‑150 (MQ=255)
                                                                                                       aCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                                         <  1:2243359/150‑1 (MQ=255)
                                                                                                       aCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                                         >  2:3679180/1‑150 (MQ=255)
                                                                                                        ccTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGa                                        >  1:3131164/1‑150 (MQ=255)
                                                                                                         cTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGAt                                       >  2:4073613/1‑150 (MQ=255)
                                                                                                          tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      <  2:129341/150‑1 (MQ=255)
                                                                                                          tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      >  2:5225876/1‑150 (MQ=255)
                                                                                                          tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      >  1:4359458/1‑150 (MQ=255)
                                                                                                          tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                      <  2:3907903/150‑1 (MQ=255)
                                                                                                               gTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTca                                 <  2:81103/150‑1 (MQ=255)
                                                                                                               gTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTca                                 >  2:2133521/1‑150 (MQ=255)
                                                                                                               gTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTca                                 >  1:563365/1‑150 (MQ=255)
                                                                                                                ttGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaa                                <  1:5284160/150‑1 (MQ=255)
                                                                                                                    caccaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAg                            >  1:5170900/1‑150 (MQ=255)
                                                                                                                     accaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAgg                           <  2:2769448/150‑1 (MQ=255)
                                                                                                                        aGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCgg                        >  1:4705327/1‑150 (MQ=255)
                                                                                                                          ccTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaa                      <  1:4173577/150‑1 (MQ=255)
                                                                                                                            tGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaata                    <  1:4240425/150‑1 (MQ=255)
                                                                                                                            tGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaata                    >  2:5164291/1‑150 (MQ=255)
                                                                                                                             gggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGactac                   >  1:4641972/1‑149 (MQ=255)
                                                                                                                              ggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat                  <  2:3937913/150‑1 (MQ=255)
                                                                                                                                 cAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCtt               <  1:2493785/150‑1 (MQ=255)
                                                                                                                                 cAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCtt               >  1:190335/1‑150 (MQ=255)
                                                                                                                                  aaaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCttt              <  1:2877207/150‑1 (MQ=255)
                                                                                                                                  aaaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCttt              <  1:5091505/150‑1 (MQ=255)
                                                                                                                                  aaaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCttt              >  2:5190723/1‑150 (MQ=255)
                                                                                                                                           cgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt     <  1:2354873/150‑1 (MQ=255)
                                                                                                                                            gtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATa    <  2:1826240/150‑1 (MQ=255)
                                                                                                                                             tATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAc   <  1:3192933/150‑1 (MQ=255)
                                                                                                                                              aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc  >  1:2360548/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACC  >  CP009273/662485‑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: