Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,948,951 C→T 5.0% M20I (ATG→ATA)  amiC ← N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,948,9510CT 5.0% 272.2 / 3.4 80M20I (ATG→ATA) amiCN‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base T (2/2);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

AGGGTTACTCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAATCTGTCA  >  NC_000913/2948807‑2949099
                                                                                                                                                |                                                                                                                                                    
aGGGTTACTCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATgg                                                                                                                                                  >  4:792806/1‑149 (MQ=255)
  ggTTACTCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCa                                                                                                                                                >  2:420778/1‑149 (MQ=255)
        tCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCcgcg                                                                                                                                          >  2:80142/1‑149 (MQ=255)
        tCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCcgcg                                                                                                                                          >  2:717199/1‑149 (MQ=255)
        tCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCcgcg                                                                                                                                          >  4:483265/1‑149 (MQ=255)
          aaCGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGcc                                                                                                                                        >  3:260370/1‑149 (MQ=255)
            cgcgAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCtt                                                                                                                                      <  2:200582/149‑1 (MQ=255)
             gcgAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTg                                                                                                                                     >  3:659947/1‑149 (MQ=255)
               gAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCa                                                                                                                                   <  2:635054/149‑1 (MQ=255)
                aaCTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAg                                                                                                                                  >  3:332960/1‑149 (MQ=255)
                 aCTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGt                                                                                                                                 <  2:7066/149‑1 (MQ=255)
                 aCTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGt                                                                                                                                 >  4:156551/1‑149 (MQ=255)
                   tGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTaa                                                                                                                               >  1:16581/1‑149 (MQ=255)
                     cTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAac                                                                                                                             >  3:497033/1‑149 (MQ=255)
                       tataTTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAacga                                                                                                                           <  2:677931/149‑1 (MQ=255)
                        atatTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAacgac                                                                                                                          <  4:659947/149‑1 (MQ=255)
                                  tGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGAt                                                                                                                >  2:516105/1‑149 (MQ=255)
                                   gACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGAtt                                                                                                               <  3:674125/149‑1 (MQ=255)
                                      gATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCa                                                                                                            >  3:575287/1‑149 (MQ=255)
                                      gATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCa                                                                                                            <  3:398782/149‑1 (MQ=255)
                                        tttGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAgt                                                                                                          >  2:394444/1‑149 (MQ=255)
                                         ttGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAgtg                                                                                                         >  3:632583/1‑149 (MQ=255)
                                           gATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGtt                                                                                                       <  1:160714/149‑1 (MQ=255)
                                             ttCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTgg                                                                                                     <  4:632583/149‑1 (MQ=255)
                                              tCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAacgacga                                                                                                                        <  4:747731/129‑1 (MQ=255)
                                              tCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAacgacga                                                                                                                        >  3:747731/1‑129 (MQ=255)
                                              tCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGa                                                                                                    <  4:12609/149‑1 (MQ=255)
                                              tCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGa                                                                                                    >  3:733119/1‑149 (MQ=255)
                                                cACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATc                                                                                                  <  4:733119/149‑1 (MQ=255)
                                                 aCCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCATATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATcc                                                                                                 >  4:684322/1‑149 (MQ=255)
                                                 aCCGTCACGCGGGTGCAGCTGTACCCAGGAAAGACTCGCACCGCCACGACCTGGCTGACCGCAGCCAGTCTGACTTGACGTACGCGCAACAGCCATATGGCACCCGCGCCTTACAGTAAGCGCCGACGGCTGATTGCAGTGTTGGATcc                                                                                                 >  4:496154/1‑149 (MQ=255)
                                                   cGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTg                                                                                               <  4:497033/149‑1 (MQ=255)
                                                    gTCACGCGGGGGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCAAATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGa                                                                                              <  2:94924/149‑1 (MQ=255)
                                                     tCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGAc                                                                                             <  4:606705/149‑1 (MQ=255)
                                                     tCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGAc                                                                                             <  1:717199/149‑1 (MQ=255)
                                                      cACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACa                                                                                            <  3:792806/149‑1 (MQ=255)
                                                        cgcgGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATg                                                                                          <  4:256038/149‑1 (MQ=255)
                                                         gcgGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGc                                                                                         <  1:555269/149‑1 (MQ=255)
                                                         gcgGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGc                                                                                         <  3:483265/149‑1 (MQ=255)
                                                           gggTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCct                                                                                       <  3:629581/149‑1 (MQ=255)
                                                                     ggACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCa                                                                             >  4:800393/1‑149 (MQ=255)
                                                                      gACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCaa                                                                            >  1:783347/1‑149 (MQ=255)
                                                                       aCGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCaaa                                                                           <  1:420778/149‑1 (MQ=255)
                                                                        cGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCaaaa                                                                          <  2:16581/149‑1 (MQ=255)
                                                                        cGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCaaaa                                                                          <  4:332960/149‑1 (MQ=255)
                                                                        cGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCaaaa                                                                          <  3:800393/149‑1 (MQ=255)
                                                                            ggCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGaaa                                                                      <  1:394444/149‑1 (MQ=255)
                                                                             gCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAAt                                                                     >  4:44622/1‑149 (MQ=255)
                                                                               cAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCt                                                                   >  3:192565/1‑149 (MQ=255)
                                                                                aGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTa                                                                  >  3:447442/1‑149 (MQ=255)
                                                                                   cgcgCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGc                                                               <  1:80142/149‑1 (MQ=255)
                                                                                   cgcgCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGc                                                               <  1:441053/149‑1 (MQ=255)
                                                                                      gcACCGCCGCGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCATATGGCACCCGCGCCTCGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCtgt                                                            <  3:496154/149‑1 (MQ=255)
                                                                                      gcACCGCCACGACGTGGCTGACCGCAGCCAGACTGACCTGACTTACCCTCAATAGCCACATGGCAc                                                                                                                                               <  4:679667/66‑1 (MQ=255)
                                                                                      gcACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCATATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCtgt                                                            <  3:684322/149‑1 (MQ=255)
                                                                                      gcACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCAc                                                                                                                                               >  3:679667/1‑66 (MQ=255)
                                                                                      gcACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCtgt                                                            <  4:110465/149‑1 (MQ=255)
                                                                                        aCCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTgtgt                                                          <  1:553953/149‑1 (MQ=255)
                                                                                             cACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCt                                                                                                >  1:313491/1‑106 (MQ=255)
                                                                                             cACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCt                                                                                                <  2:313491/106‑1 (MQ=255)
                                                                                             cACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACaa                                                     <  4:667243/149‑1 (MQ=255)
                                                                                              aCGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAacgacga                                                                                                                        >  4:426536/1‑81 (MQ=255)
                                                                                              aCGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAacgacga                                                                                                                        <  3:426536/81‑1 (MQ=255)
                                                                                              aCGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAg                                                    >  3:39592/1‑149 (MQ=255)
                                                                                                   cTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAAc                                               <  1:188625/149‑1 (MQ=255)
                                                                                                      gCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGAc                                            <  4:192565/149‑1 (MQ=255)
                                                                                                                gCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGAtt                                                                                                               <  1:249741/72‑1 (MQ=255)
                                                                                                                gCCAGACTGACCTGACTTACGCTCAATAGCAACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGAtt                                                                                                               >  2:249741/1‑72 (MQ=255)
                                                                                                                           ctgacTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAAc                                               >  2:286032/1‑125 (MQ=255)
                                                                                                                           ctgacTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAAc                                               <  1:286032/125‑1 (MQ=255)
                                                                                                                            tgacTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAAc                            <  3:544733/143‑1 (MQ=255)
                                                                                                                            tgacTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAAc                            >  4:544733/1‑143 (MQ=255)
                                                                                                                                     cTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACtt                          >  4:453411/1‑136 (MQ=255)
                                                                                                                                     cTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACtt                          <  3:453411/136‑1 (MQ=255)
                                                                                                                                      tCAATAGCCACATGGCACCTGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCa            >  4:552840/1‑149 (MQ=255)
                                                                                                                                        aaTAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATa          >  4:179907/1‑149 (MQ=255)
                                                                                                                                        aaTAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATa          >  2:60798/1‑149 (MQ=255)
                                                                                                                                          tAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAAt        <  4:560746/149‑1 (MQ=255)
                                                                                                                                          tAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAAt        >  4:319949/1‑149 (MQ=255)
                                                                                                                                             ccACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAATCTg     <  3:552840/149‑1 (MQ=255)
                                                                                                                                               acaTGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAATCTGTc   <  4:447442/149‑1 (MQ=255)
                                                                                                                                               acaTGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAATCTGTc   <  4:39592/149‑1 (MQ=255)
                                                                                                                                                caTGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAATCTGTCa  <  3:44622/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGGGTTACTCAACGCGAACTGCTTATATTTCAGCTGACGATTTGATTCCACCGTCACGCGGGTGTAGCTGGACGCAGGCCAGACGCGCACCGCCACGACCTGGCTGACCGCAGCCAGACTGACCTGACTTACGCTCAATAGCCACATGGCACCCGCGCCTTGCAGTAAACGACGACGGCTGATTGCAGTGTTGGATCCTGACATGCCTCTCCCGAGCAAAAGAAATCTAAGCTGTGTAACAAGTAAACGACTAATTTGACCGGAAACTTTAGCGAAAGACGCATAATCTGTCA  >  NC_000913/2948807‑2949099

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑