Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,523,931 T→G 5.7% V354G (GTC→GGC)  mrcA → penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,523,9310TG 5.7% 335.0 / 3.0 115V354G (GTC→GGC) mrcApenicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (63/43);  new base G (3/5);  total (67/48)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

AATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGC  >  NC_000913/3523783‑3524077
                                                                                                                                                    |                                                                                                                                                  
aataaCGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTc                                                                                                                                                   >  1:1322855/1‑150 (MQ=255)
   aaCGTGCTGGACTTCGACATGAGCCACGGCTATCGCGGACCGGCAAATGTGCTGTGGAAAGTTGGCGAGTCGGCGTTGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCAcc                                                                                                                                                <  2:1591279/150‑1 (MQ=255)
   aaCGTGCTGGACTACGACATGCGCCTCGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCAcc                                                                                                                                                <  2:739070/150‑1 (MQ=255)
    aCGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCCGTCAccc                                                                                                                                               >  1:436694/1‑149 (MQ=255)
     cGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAg                                                                                                                                              >  1:1624592/1‑150 (MQ=255)
          tGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCa                                                                                                                                         <  1:1531541/150‑1 (MQ=255)
           ggACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCaa                                                                                                                                        >  2:19963/1‑150 (MQ=255)
             aCTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATc                                                                                                                                      >  1:398998/1‑150 (MQ=255)
              cTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATcc                                                                                                                                     >  1:867551/1‑150 (MQ=255)
                aCGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTc                                                                                                                                   <  1:1282460/150‑1 (MQ=255)
                aCGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTc                                                                                                                                   <  1:915619/150‑1 (MQ=255)
                aCGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTc                                                                                                                                   <  2:608276/150‑1 (MQ=255)
                    cATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCa                                                                                                                               <  1:1089387/150‑1 (MQ=255)
                        cgcCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCg                                                                                                                           <  2:1322855/150‑1 (MQ=255)
                             cGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGGTGCATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGGCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGc                                                                                                                      <  1:438528/150‑1 (MQ=255)
                              ggCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCCGCAAGCGACGGCg                                                                                                                     >  1:122135/1‑150 (MQ=255)
                                 tATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATg                                                                                                                  <  2:597307/150‑1 (MQ=255)
                                  aTCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATga                                                                                                                 >  1:7453/1‑149 (MQ=255)
                                     gcgGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGTCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTgg                                                                                                              >  2:615575/1‑150 (MQ=255)
                                      gatctcGGGAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGc                                                                                                             <  1:585043/145‑1 (MQ=255)
                                       ggCCCGGCAAATGTGCTGTGGAAAGGGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCg                                                                                                            <  1:1470880/150‑1 (MQ=255)
                                         cccGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGa                                                                                                          >  1:944116/1‑150 (MQ=255)
                                         cccGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGa                                                                                                          >  2:1239593/1‑150 (MQ=255)
                                         cccGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGa                                                                                                          >  2:1016572/1‑150 (MQ=255)
                                            ggCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACgg                                                                                                       >  1:700016/1‑150 (MQ=255)
                                              cAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGt                                                                                                     <  2:485131/150‑1 (MQ=255)
                                               aaaTGTGCTGTGGAAAGGGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTc                                                                                                    <  1:557576/150‑1 (MQ=255)
                                                   gtgCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGAcc                                                                                                >  2:1529189/1‑150 (MQ=255)
                                                     gCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGt                                                                                              >  1:1255796/1‑150 (MQ=255)
                                                       tggggAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCg                                                                                            <  2:523029/147‑1 (MQ=255)
                                                        gtgGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGc                                                                                           >  1:759610/1‑150 (MQ=255)
                                                         tgGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACAGCGATGCTGGCGGACGGGTCGACCGTCGCa                                                                                          >  2:1243011/1‑150 (MQ=255)
                                                         gggAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGCAGGCGCTGCCAACCTATGGGCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCa                                                                                          <  1:903210/149‑1 (MQ=255)
                                                           gAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCAtt                                                                                        >  2:1135958/1‑150 (MQ=255)
                                                            aaaGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTg                                                                                       >  2:1484297/1‑150 (MQ=255)
                                                              aGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAg                                                                                     <  1:981678/150‑1 (MQ=255)
                                                               gTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGt                                                                                    <  1:9920/150‑1 (MQ=255)
                                                               gTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGAGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGCCGGGTCGACCGTCGCATTGAGt                                                                                    >  2:170187/1‑150 (MQ=255)
                                                                tGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTa                                                                                   <  1:537960/150‑1 (MQ=255)
                                                                    cGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGa                                                                               <  1:144141/150‑1 (MQ=255)
                                                                    cGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCAGTGAGTTTGGa                                                                               >  2:455449/1‑150 (MQ=255)
                                                                      aGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAg                                                                             >  1:1088662/1‑150 (MQ=255)
                                                                       gTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAgg                                                                            <  2:17983/150‑1 (MQ=255)
                                                                       gTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAgg                                                                            <  1:1016572/150‑1 (MQ=255)
                                                                          ggCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGt                                                                         <  2:867551/150‑1 (MQ=255)
                                                                          ggCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGt                                                                         >  1:1045924/1‑150 (MQ=255)
                                                                          ggCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTAAGGTCAGCTGCTTCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGAGACGCCGATGATGCCGGACGGGTCGACGGGCGCATTGAGGTTTGAAGGcct                                                                         >  2:613175/1‑148 (MQ=255)
                                                                           gCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACTGGTCGACCGTCGCATTGAGTATGGAAGGCGtt                                                                        >  2:312446/1‑150 (MQ=255)
                                                                              tGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGc                                                                     >  2:429049/1‑150 (MQ=255)
                                                                               gggATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCCGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCt                                                                    >  1:1542277/1‑150 (MQ=255)
                                                                               gggATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGTCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCt                                                                    >  1:710734/1‑150 (MQ=255)
                                                                               gggATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGCAGGCGTTCGCt                                                                    >  1:1535514/1‑150 (MQ=255)
                                                                                ggATAACTACACGTATACAGATACGCTGAAGGCGCTGCCAACCCATGGGCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCCGCTAGCGACGGCGATGCTGGCGGACGGGGCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTg                                                                   <  1:328764/150‑1 (MQ=255)
                                                                                ggATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTTTGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTg                                                                   <  1:1411687/150‑1 (MQ=255)
                                                                                 gATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTgg                                                                  >  1:485193/1‑150 (MQ=255)
                                                                                 gATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCACATTGAGTATGGAAGGCGTTCGCTgg                                                                  <  2:1484176/150‑1 (MQ=255)
                                                                                    aacaacAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGCTGCTGGCGGACGGGTCGACCGTCGCAGTGAGTATGGAAGGCGTTCGCTGGgcg                                                               >  2:1148565/1‑150 (MQ=255)
                                                                                      caacaaGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGCCGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGcgcg                                                             >  1:447349/1‑150 (MQ=255)
                                                                                      caacaaGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGCCGGGTCGCCCGTCGCATTGAGGATGGAAGGCGTTCGCTGGGcgcg                                                             >  2:850261/1‑150 (MQ=255)
                                                                                      caacaaGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGcgcg                                                             >  2:1547749/1‑150 (MQ=255)
                                                                                      caacaaGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAAACCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGcgcg                                                             <  2:825777/150‑1 (MQ=255)
                                                                                       aacaaGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATACTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGt                                                            >  2:762581/1‑150 (MQ=255)
                                                                                            gATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTa                                                       >  1:939956/1‑150 (MQ=255)
                                                                                            gATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTCCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGACCTTa                                                       >  2:822791/1‑150 (MQ=255)
                                                                                              ttACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTAcc                                                     >  1:1000640/1‑150 (MQ=255)
                                                                                                aCCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGt                                                   >  1:1133102/1‑150 (MQ=255)
                                                                                                  cGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCACGCGACGGCGATGCTGGCGGCCGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGtgc                                                 >  1:1604431/1‑148 (MQ=255)
                                                                                                  cGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGCCCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTc                                                 >  2:759732/1‑150 (MQ=255)
                                                                                                   gATCCGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCg                                                <  1:442530/150‑1 (MQ=255)
                                                                                                    aTACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCgg                                               <  2:700016/150‑1 (MQ=255)
                                                                                                    aTACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGGCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCgg                                               <  1:850261/150‑1 (MQ=255)
                                                                                                        gCTGAAGGCGCTGCCAACCTATGGTACGCTGCTGACAGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCAGACGGACGGGTCGAACGTCGCATTAAGAAAGGAAGGCGTTCGCTGGGCGCGTCCTTACCGATCGGATAc                                           >  1:1490648/1‑150 (MQ=255)
                                                                                                         cTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACt                                          >  1:205183/1‑150 (MQ=255)
                                                                                                          tGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCGGCCGCAGGCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTc                                         <  1:547931/150‑1 (MQ=255)
                                                                                                           gAAGGCGCGGCCAACCTATGGACCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGTGTCGAACGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGGCCTTACCGTTCGGATACTCa                                        <  2:1134038/150‑1 (MQ=255)
                                                                                                            aaGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAg                                       <  2:1255796/150‑1 (MQ=255)
                                                                                                             aGGCGCTGCCAACCACTGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCCGACGGGGCGACCGTCGCATTGAGTATGGAAGGCGCTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGc                                      <  2:436694/150‑1 (MQ=255)
                                                                                                              ggCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCa                                     <  1:1239593/150‑1 (MQ=255)
                                                                                                              ggCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGGCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCa                                     <  1:1077404/150‑1 (MQ=255)
                                                                                                                cgcTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAg                                   >  2:249193/1‑150 (MQ=255)
                                                                                                                  cTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGa                                 <  1:789501/150‑1 (MQ=255)
                                                                                                                  cTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGa                                 >  2:1395832/1‑150 (MQ=255)
                                                                                                                   tGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGAc                                >  2:1455755/1‑150 (MQ=255)
                                                                                                                    gCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGAcc                               >  2:1677073/1‑150 (MQ=255)
                                                                                                                     ccAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCg                              >  2:634908/1‑150 (MQ=255)
                                                                                                                     ccAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCg                              >  2:1559050/1‑150 (MQ=255)
                                                                                                                     ccAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAACAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTg                                                                   <  1:27077/113‑1 (MQ=255)
                                                                                                                      cAACCTATGGTCCGCTGCTGCCTGCCGCGGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGa                             >  1:1189668/1‑150 (MQ=255)
                                                                                                                           tATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCTGATACTCAGCAAGGACCGAcgccg                        <  1:762581/150‑1 (MQ=255)
                                                                                                                                ggcGCTGGTGACTGCGGCAGGCAGCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTaa                   >  1:1228703/3‑150 (MQ=255)
                                                                                                                                  cGCTGGTGACTGCGGCAGGCAGCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAg                 >  1:1484814/1‑150 (MQ=255)
                                                                                                                                  cGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCCGCAAGGACCGACGCCGCGTAAAg                 >  1:879936/1‑150 (MQ=255)
                                                                                                                                  cGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAg                 >  2:1325214/1‑150 (MQ=255)
                                                                                                                                    ctgctgCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGCCGACGCGTAAAGTg               >  1:140745/1‑150 (MQ=255)
                                                                                                                                      ggtgactgcgGCAGGCAGCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGAc             >  2:738403/6‑150 (MQ=255)
                                                                                                                                           ctgccGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATg        >  1:1530925/1‑150 (MQ=255)
                                                                                                                                           cggccGCAGGCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGGCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATg        <  2:179630/148‑1 (MQ=255)
                                                                                                                                              ccGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTc     <  2:122135/150‑1 (MQ=255)
                                                                                                                                              ccGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTc     >  1:800640/1‑150 (MQ=255)
                                                                                                                                                 cAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCAACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGc  <  1:615575/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGC  >  NC_000913/3523783‑3524077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: