Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,524,184 C→T A438A (GCC→GCT 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,524,1840CT95.6% 1260.3 / 31.7 388A438A (GCC→GCTmrcApenicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (6/11);  new base T (164/205);  total (170/218)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 388 total aligned reads displayed.

TCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAAC  >  NC_000913/3524037‑3524331
                                                                                                                                                   |                                                                                                                                                   
 cAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGtt                                                                                                                                                  <  5:606713/150‑1 (MQ=255)
 cAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGtt                                                                                                                                                  >  2:779259/1‑150 (MQ=255)
 cAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGt                                                                                                                                                   <  4:218358/149‑1 (MQ=255)
 cAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGt                                                                                                                                                   <  3:684033/149‑1 (MQ=255)
 cAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTa                                                                                                                                                 <  3:156845/151‑1 (MQ=255)
   gCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATg                                                                                                                                               <  5:368756/151‑1 (MQ=255)
    cAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATgg                                                                                                                                              <  6:868124/151‑1 (MQ=255)
      aGGCCCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGACGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGCCGCGGGTGTCGATCCATCCGCAAAACGGTGCTGTTATGgcg                                                                                                                                            >  3:405119/1‑151 (MQ=255)
         aCCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTg                                                                                                                                         <  6:1630975/151‑1 (MQ=255)
            gACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTc                                                                                                                                      <  5:1403037/151‑1 (MQ=255)
             aCGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCg                                                                                                                                     >  8:716773/1‑151 (MQ=255)
             aCGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCg                                                                                                                                     <  4:28974/151‑1 (MQ=255)
               gccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGATCTCGGGGCTGGTGTCGATCAATCCGCTAAACGGTGCTGTTATGGCGCTGGTCGGt                                                                                                                                   <  6:941255/151‑1 (MQ=255)
                 cgcgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTgg                                                                                                                                 <  4:1635671/151‑1 (MQ=255)
                 cgcgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTgg                                                                                                                                 <  7:1165036/151‑1 (MQ=255)
                 cgcgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTgg                                                                                                                                 >  6:1049925/1‑151 (MQ=255)
                 cgcgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTgg                                                                                                                                 <  5:655631/151‑1 (MQ=255)
                     tAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCttt                                                                                                                             >  5:502936/1‑151 (MQ=255)
                      aaaGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTg                                                                                                                            <  7:498724/151‑1 (MQ=255)
                      aaaGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTg                                                                                                                            <  5:1142371/151‑1 (MQ=255)
                      aaaGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTCGTCTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGCTCGGTGGCTTTg                                                                                                                            >  7:138735/1‑151 (MQ=255)
                         gTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCCTGGTGGCTGGCACAAGTGCCGCAAGAGACCTCGGCGCTGCTGTCGATCAATCCGCAAAACGCTGCTGTTATGCCGCAGGGCGGTGGCTTTGAtt                                                                                                                         >  8:52696/1‑151 (MQ=255)
                           gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAttt                                                                                                                        <  4:1434671/150‑1 (MQ=255)
                           gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAttt                                                                                                                        <  3:991182/150‑1 (MQ=255)
                           gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTc                                                                                                                       <  4:1072791/151‑1 (MQ=255)
                            aCCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTc                                                                                                                       >  7:1286802/1‑150 (MQ=255)
                            aCCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCa                                                                                                                      <  4:587689/151‑1 (MQ=255)
                             ccGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCa                                                                                                                      >  6:1494665/1‑150 (MQ=255)
                              cGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAAt                                                                                                                    >  1:455891/1‑151 (MQ=255)
                                  gTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCaga                                                                                                                <  1:129279/151‑1 (MQ=255)
                                  gTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCag                                                                                                                 <  7:597388/150‑1 (MQ=255)
                                       gCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGc                                                                                                              <  5:173342/148‑1 (MQ=255)
                                        cAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGt                                                                                                          <  2:445191/151‑1 (MQ=255)
                                          aaCGGGTCAGCAAATCTGGTTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCCGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTCGATTGCATTTAAAACACAGCAAGttt                                                                                                        >  3:1003801/1‑151 (MQ=255)
                                                 cAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCGGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACcgcg                                                                                                 <  4:871429/151‑1 (MQ=255)
                                                 cAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACcgcg                                                                                                 <  6:1314283/151‑1 (MQ=255)
                                                  aGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCgcgc                                                                                                >  7:938421/1‑151 (MQ=255)
                                                    cAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGcc                                                                                               >  1:316210/1‑150 (MQ=255)
                                                        tCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCAccc                                                                                           >  7:1647085/1‑150 (MQ=255)
                                                           gggTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCa                                                                                          <  8:618545/148‑1 (MQ=255)
                                                           gggTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAg                                                                                         <  7:1459813/149‑1 (MQ=255)
                                                           gggTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAg                                                                                         <  8:1399210/149‑1 (MQ=255)
                                                             gTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGctc                                                                                     >  6:454833/1‑149 (MQ=255)
                                                               tCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACt                                                                                    <  1:222387/150‑1 (MQ=255)
                                                               tCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACt                                                                                    <  1:898769/150‑1 (MQ=255)
                                                               tCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACt                                                                                    <  4:1266058/150‑1 (MQ=255)
                                                                cGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGc                                                                                  <  4:1541230/151‑1 (MQ=255)
                                                                 gTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGAGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCg                                                                                 >  8:1441762/1‑151 (MQ=255)
                                                                    aGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTc                                                                               >  6:1215695/1‑150 (MQ=255)
                                                                    aGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTc                                                                               >  6:1115916/1‑150 (MQ=255)
                                                                    aGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTc                                                                               >  1:1226147/1‑150 (MQ=255)
                                                                    aGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCa                                                                              <  3:790627/151‑1 (MQ=255)
                                                                        tGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGt                                                                           <  8:27107/150‑1 (MQ=255)
                                                                        tGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTg                                                                          >  6:183306/1‑151 (MQ=255)
                                                                        tGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTg                                                                          >  3:877247/1‑151 (MQ=255)
                                                                          gCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGCTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTCGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTggg                                                                        >  8:1668040/1‑151 (MQ=255)
                                                                           cGATGCATGGTGGCTGGCACTAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGt                                                                       <  1:564804/151‑1 (MQ=255)
                                                                            gATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGtt                                                                      >  1:1453843/1‑151 (MQ=255)
                                                                                 aTGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAAc                                                                 >  6:1010599/1‑151 (MQ=255)
                                                                                 aTGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAAc                                                                 <  3:518126/151‑1 (MQ=255)
                                                                                 aTGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAAc                                                                 <  8:1049399/151‑1 (MQ=255)
                                                                                   ggtggCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACAt                                                               <  1:159068/151‑1 (MQ=255)
                                                                                   ggtggCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACAt                                                               >  2:357668/1‑151 (MQ=255)
                                                                                     tggctggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATc                                                              >  5:851746/1‑150 (MQ=255)
                                                                                      ggctggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCaa                                                            <  2:1507857/151‑1 (MQ=255)
                                                                                      ggctggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCa                                                             <  8:341716/150‑1 (MQ=255)
                                                                                         gggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAAcc                                                         <  1:134406/150‑1 (MQ=255)
                                                                                          ggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGCTCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCTCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCACACCg                                                        >  4:436088/1‑151 (MQ=255)
                                                                                          ggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAAcc                                                         <  4:1667847/150‑1 (MQ=255)
                                                                                            cacaAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGtt                                                      <  3:1092162/151‑1 (MQ=255)
                                                                                             acaAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGAACTGCGTCAGGTGGGTTCCAACATCAAACCGtt                                                      >  2:612499/1‑150 (MQ=255)
                                                                                                  tGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCctct                                                 >  7:1239010/1‑150 (MQ=255)
                                                                                                    ccGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTaca                                              <  8:1580265/151‑1 (MQ=255)
                                                                                                     cGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTacac                                             <  2:1682216/151‑1 (MQ=255)
                                                                                                        aaGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTCCCATTGATTTAAATCACAGCAAGTTTAACCTCGTCAACCAGGGACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTacaaccc                                          >  1:716004/1‑147 (MQ=255)
                                                                                                         aGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGTGTTCCAACATCAAACCGTTCCTCTACACcgcg                                         <  1:1079276/151‑1 (MQ=255)
                                                                                                         aGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACcgcg                                         >  3:1414346/1‑151 (MQ=255)
                                                                                                         aGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGGGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACcgcg                                         <  3:1191805/151‑1 (MQ=255)
                                                                                                            gAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCgcggc                                       <  4:753220/150‑1 (MQ=255)
                                                                                                              aCTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGAt                                    >  2:1112135/1‑151 (MQ=255)
                                                                                                                    cgcTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATa                               >  7:376178/1‑150 (MQ=255)
                                                                                                                    cgcTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATa                               >  6:285928/1‑150 (MQ=255)
                                                                                                                    cgcTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATCGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAACCAGTTCCTCTCCACGGCGGCGATGGaga                               >  2:1179514/1‑148 (MQ=255)
                                                                                                                    cgcTGGTGTCGATCAATCCGCAAAACGGTGCCGGTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                              <  4:1344420/151‑1 (MQ=255)
                                                                                                                     gcTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCgcgg                                        >  5:411120/1‑140 (MQ=255)
                                                                                                                     gcTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaaa                             <  7:1573997/151‑1 (MQ=255)
                                                                                                                     gcTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaaa                             <  3:1270687/151‑1 (MQ=255)
                                                                                                                        ggggtCGATCAATCCGCAAACCGGTGCGGGTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAg                            <  4:985998/146‑1 (MQ=255)
                                                                                                                         gtgtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTc                         <  3:1310687/151‑1 (MQ=255)
                                                                                                                           gtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                       <  4:636999/151‑1 (MQ=255)
                                                                                                                           gtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                       <  4:1500696/151‑1 (MQ=255)
                                                                                                                           gtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                       >  2:1382586/1‑151 (MQ=255)
                                                                                                                               aTCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGc                   <  7:237321/151‑1 (MQ=255)
                                                                                                                                 cAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCt                  <  7:751410/150‑1 (MQ=255)
                                                                                                                                 cAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCt                  <  6:1363289/150‑1 (MQ=255)
                                                                                                                                  aaTCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTgg                <  7:1668991/151‑1 (MQ=255)
                                                                                                                                     ccGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCa              <  7:321479/150‑1 (MQ=255)
                                                                                                                                     ccGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCa              <  2:1583283/150‑1 (MQ=255)
                                                                                                                                         aaaaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTAt         <  5:1461208/151‑1 (MQ=255)
                                                                                                                                              ggTGCTGTTAAGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTg     >  5:1529253/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAAC  >  NC_000913/3524037‑3524331

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