Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 662,631 C→A 91.2% 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,6310CA91.2% 864.9 / 69.2 283R348L (CGC→CTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base C (10/15);  new base A (142/116);  total (152/131)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 285 total aligned reads displayed.

TAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATACCG  >  CP009273/662484‑662777
                                                                                                                                                   |                                                                                                                                                  
tAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTgagc                                                                                                                                                  >  2:268259/1‑150 (MQ=255)
    ccACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGa                                                                                                                                              >  1:3178999/1‑150 (MQ=255)
     cACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGAt                                                                                                                                             <  1:1051731/150‑1 (MQ=255)
     cACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGAt                                                                                                                                             <  1:2589143/150‑1 (MQ=255)
     cACCGGATAGAAGAAGGTATCCGCAGATTCTTCCCGCGATCTTGTGACCTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTGTCCGAACCTGGCAGTTGCCACCCGCCTGGGTCAAACAGCGTCGTATTGAGCGTGAt                                                                                                                                             <  2:2420997/150‑1 (MQ=255)
       ccTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCa                                                                                                                                           >  1:2945314/1‑150 (MQ=255)
       acTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCa                                                                                                                                           >  1:551104/2‑150 (MQ=255)
        cTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAc                                                                                                                                          >  2:2593295/1‑150 (MQ=255)
         tGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCAcc                                                                                                                                         >  2:1939854/1‑150 (MQ=255)
            tAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCg                                                                                                                                      <  1:3209140/150‑1 (MQ=255)
             agaagaAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCgg                                                                                                                                     >  1:4399653/1‑150 (MQ=255)
              gaagaaGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAGCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgc                                                                                                                                    >  1:16794/1‑150 (MQ=255)
               aagaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgcg                                                                                                                                   >  1:689094/1‑150 (MQ=255)
                agaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGgcgc                                                                                                                                  >  1:903433/1‑150 (MQ=255)
                 gaagGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCt                                                                                                                                 >  1:3415264/1‑150 (MQ=255)
                  aagGTATCCACAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTc                                                                                                                                >  1:3113477/1‑150 (MQ=255)
                        tCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGcc                                                                                                                          >  1:473795/1‑150 (MQ=255)
                             aGATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAAc                                                                                                                     <  2:63887/150‑1 (MQ=255)
                              gATTCTTCCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAAcc                                                                                                                    <  1:1005676/150‑1 (MQ=255)
                                  cttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCca                                                                                                                >  1:3335624/1‑150 (MQ=255)
                                  cttcCAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCca                                                                                                                >  1:2958700/1‑150 (MQ=255)
                                     ccAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt                                                                                                             <  1:3268254/150‑1 (MQ=255)
                                      cAGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATa                                                                                                            >  2:44178/1‑150 (MQ=255)
                                       aGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAACCCGCCACATAg                                                                                                           >  2:2663092/1‑150 (MQ=255)
                                       aGCGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAg                                                                                                           >  2:2890909/1‑150 (MQ=255)
                                         cGATCTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAggg                                                                                                         <  2:1998958/150‑1 (MQ=255)
                                             cTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTa                                                                                                     <  1:1469010/150‑1 (MQ=255)
                                             cTTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTa                                                                                                     <  1:3034956/150‑1 (MQ=255)
                                              ttGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTaa                                                                                                    >  1:17216/1‑150 (MQ=255)
                                              ttGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTaa                                                                                                    >  1:2420092/1‑150 (MQ=255)
                                                   aCATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTa                                                                                               >  2:2174777/1‑150 (MQ=255)
                                                   aCATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTa                                                                                               >  2:3298466/1‑150 (MQ=255)
                                                   aCATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTa                                                                                               >  2:4447505/1‑150 (MQ=255)
                                                    cATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAg                                                                                              <  2:1103564/150‑1 (MQ=255)
                                                         aGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGAcgcg                                                                                         >  2:2720398/1‑150 (MQ=255)
                                                           aCGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCggg                                                                                       <  2:4518421/150‑1 (MQ=255)
                                                            cGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGgga                                                                                      <  1:2512292/150‑1 (MQ=255)
                                                             gCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGggag                                                                                     >  2:1020173/1‑150 (MQ=255)
                                                                  tGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaa                                                                                >  1:4396335/1‑150 (MQ=255)
                                                                   ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATaaa                                                                               >  1:410516/1‑150 (MQ=255)
                                                                    gCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAc                                                                              <  1:2603423/150‑1 (MQ=255)
                                                                     ccccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAcc                                                                             >  1:1163859/1‑150 (MQ=255)
                                                                     ccccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAcc                                                                             >  1:50239/1‑150 (MQ=255)
                                                                      cccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAccc                                                                            <  1:3190616/150‑1 (MQ=255)
                                                                        cATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAccccc                                                                          <  2:1028597/150‑1 (MQ=255)
                                                                         aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                                         >  1:908508/1‑150 (MQ=255)
                                                                         aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                                         >  2:2547145/1‑150 (MQ=255)
                                                                         aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCt                                                                         >  2:4473550/1‑150 (MQ=255)
                                                                            ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCtgtg                                                                      <  2:4023170/150‑1 (MQ=255)
                                                                               tCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTgg                                                                   >  1:1295205/1‑150 (MQ=255)
                                                                                      aCGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGtt                                                            >  1:1929864/1‑150 (MQ=255)
                                                                                       cGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                                           >  2:2513080/1‑150 (MQ=255)
                                                                                       cGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                                           <  2:214684/150‑1 (MQ=255)
                                                                                       cGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTc                                                           >  1:2992795/1‑150 (MQ=255)
                                                                                        gATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCa                                                          >  2:3681961/1‑150 (MQ=255)
                                                                                        gATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCa                                                          >  2:1361958/1‑150 (MQ=255)
                                                                                          tAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCAcc                                                        <  2:3890907/150‑1 (MQ=255)
                                                                                              gTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCg                                                    >  2:1089699/1‑150 (MQ=255)
                                                                                                   tCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTa                                               >  1:3938012/1‑150 (MQ=255)
                                                                                                     cGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTAtt                                             >  2:1457201/1‑150 (MQ=255)
                                                                                                       aaCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCg                                           >  1:1143446/1‑150 (MQ=255)
                                                                                                        aCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                                          >  1:3284203/1‑150 (MQ=255)
                                                                                                        aCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                                          <  1:3719220/150‑1 (MQ=255)
                                                                                                        aCCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCgg                                          >  1:822731/1‑150 (MQ=255)
                                                                                                           tGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATc                                       <  1:1388852/150‑1 (MQ=255)
                                                                                                               aGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTc                                   >  2:56127/1‑150 (MQ=255)
                                                                                                                   gCCACCAGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaaca                               <  2:4243581/150‑1 (MQ=255)
                                                                                                                    ccaccaGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTcaacaa                              >  1:22346/1‑150 (MQ=255)
                                                                                                                       ccaGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGc                           <  2:3094552/150‑1 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                          <  1:2074799/150‑1 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCg                          <  2:3546852/150‑1 (MQ=255)
                                                                                                                        caGCCTGGGTCAAACAGCGTCGTAATGAGAGTGATCACCCCGGCGCTCAAAGCCGAAAACCCCACAAAGTGTTTAACTGTAGACGCCGGAGGATAAACCCCCTGTGTGGCGAGGTTCACCAGCGGTGTATACGGATCGTTCAAAAAGGCg                          >  1:2802794/1‑150 (MQ=255)
                                                                                                                         aGCCTGGGTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCgg                         >  1:3086880/1‑150 (MQ=255)
                                                                                                                              gggTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataa                    <  1:3106293/150‑1 (MQ=255)
                                                                                                                               ggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat                   >  1:440645/1‑150 (MQ=255)
                                                                                                                               ggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat                   <  1:228326/150‑1 (MQ=255)
                                                                                                                               ggTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGaataat                   <  2:3747776/150‑1 (MQ=255)
                                                                                                                                gTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATc                  <  1:3586726/150‑1 (MQ=255)
                                                                                                                                gTCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGGATTCGGATCGTTCAACAAGGCGGACTAATc                  >  2:1148878/1‑150 (MQ=255)
                                                                                                                                 tCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTCGACGCGGGAGGATCAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCt                 >  2:3022644/1‑150 (MQ=255)
                                                                                                                                 tCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGACTAATCt                 >  2:3023024/1‑150 (MQ=255)
                                                                                                                                 tCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCt                 <  1:1602934/150‑1 (MQ=255)
                                                                                                                                 tCAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCt                 <  1:1207677/150‑1 (MQ=255)
                                                                                                                                  cAAACAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCtt                <  1:2496306/150‑1 (MQ=255)
                                                                                                                                   aaaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCttt               <  2:2271825/150‑1 (MQ=255)
                                                                                                                                    aaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTg              >  2:919282/1‑150 (MQ=255)
                                                                                                                                    aaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTg              <  2:602479/150‑1 (MQ=255)
                                                                                                                                    aaCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTg              >  2:3802472/1‑150 (MQ=255)
                                                                                                                                     aCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGc             <  1:4259464/150‑1 (MQ=255)
                                                                                                                                     aCAGCGTCGTATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGc             <  2:2496117/150‑1 (MQ=255)
                                                                                                                                       aGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTg           >  1:183277/1‑150 (MQ=255)
                                                                                                                                            cgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt      <  2:1388670/150‑1 (MQ=255)
                                                                                                                                            cgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt      >  1:3734711/1‑150 (MQ=255)
                                                                                                                                            cgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt      >  1:3452327/1‑150 (MQ=255)
                                                                                                                                            cgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt      >  1:2749463/1‑150 (MQ=255)
                                                                                                                                            cgtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt      <  1:2111767/150‑1 (MQ=255)
                                                                                                                                             gtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATa     <  2:1552490/150‑1 (MQ=255)
                                                                                                                                             gtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATa     <  2:2126818/150‑1 (MQ=255)
                                                                                                                                             gtATTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATa     >  2:657709/1‑150 (MQ=255)
                                                                                                                                               aTTGAGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGATCGTTCAACAAGGCGGAATAATCTTTGCTGGAGATAcc   >  2:2225290/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: