Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,868,804 A→T 6.5% intergenic (‑51/+89) nlpD ← / ← pcm activator of AmiC murein hydrolase activity, lipoprotein/L‑isoaspartate protein carboxylmethyltransferase type II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,868,8040AT 6.5% 199.6 / 4.9 62intergenic (‑51/+89)nlpD/pcmactivator of AmiC murein hydrolase activity, lipoprotein/L‑isoaspartate protein carboxylmethyltransferase type II
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base T (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

TGGATTTGAAGTGTCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGA  >  NC_000913/2868658‑2868948
                                                                                                                                                  |                                                                                                                                                
tgGATTTGAAGTGTCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAgg                                                                                                                                                <  4:404240/149‑1 (MQ=255)
    tttGAAGTGTCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCTGCGAACGGTGAATGTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCAcc                                                                                                                                            >  3:384655/1‑149 (MQ=255)
       gAAGTGTCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTa                                                                                                                                         >  4:336968/1‑149 (MQ=255)
            gtCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGaa                                                                                                                                    <  3:266325/149‑1 (MQ=255)
            gtCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGaa                                                                                                                                    <  1:56360/149‑1 (MQ=255)
             tCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAac                                                                                                                                   >  4:58411/1‑149 (MQ=255)
                 aaCAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTg                                                                                                                               <  3:509531/149‑1 (MQ=255)
                    aGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGaa                                                                                                                            >  1:128255/1‑149 (MQ=255)
                          ccagccaTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATg                                                                                                                      >  1:442838/1‑149 (MQ=255)
                           cagccaTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGt                                                                                                                     >  2:147216/1‑149 (MQ=255)
                              ccaTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCAcc                                                                                                                                            <  3:169279/123‑1 (MQ=255)
                              ccaTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCAcc                                                                                                                                            >  4:169279/1‑123 (MQ=255)
                              ccaTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCa                                                                                                                  <  2:287514/149‑1 (MQ=255)
                               caTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCAc                                                                                                                 <  3:439926/149‑1 (MQ=255)
                                aTAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACg                                                                                                                >  2:268460/1‑149 (MQ=255)
                                 tAGCGAAACCAGTGAAAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGa                                                                                                               <  4:240986/149‑1 (MQ=255)
                                 gaGCGAAAAGAGTGTCAAAGCGGCAGTGCGGCGAATGGTGTATTTTGCGCCTCCCCCGCTCATTCCTCCCACAGGAAAAATTGGTTAATAAGCAGAGACATAATTACCGTGGAAGGCAGCCTACTGAACACTGGAAAAGATGTTCAAGa                                                                                                               >  4:135395/2‑149 (MQ=255)
                                           aGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACGGAACACTGGAAAAGATGTTCACGATACGCTGAcc                                                                                                     >  3:339972/1‑149 (MQ=255)
                                           aGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTCGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGAcc                                                                                                     >  1:324674/1‑149 (MQ=255)
                                              gACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGc                                                                                                  <  2:545108/149‑1 (MQ=255)
                                               cgttcgCCGCATTGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACAGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCg                                                                                                 >  3:182713/6‑149 (MQ=255)
                                               aCAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCg                                                                                                 <  1:556821/149‑1 (MQ=255)
                                                cAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCgg                                                                                                <  2:396484/149‑1 (MQ=255)
                                                 aaaGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGc                                                                                               >  1:3834/1‑149 (MQ=255)
                                                 aaaGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGc                                                                                               >  2:185143/1‑149 (MQ=255)
                                                 aaaGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGc                                                                                               >  1:389753/1‑149 (MQ=255)
                                                       gCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATa                                                                                         >  4:89493/1‑149 (MQ=255)
                                                        cAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATaa                                                                                        <  3:41011/149‑1 (MQ=255)
                                                         aaTGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCATGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGc                                                                                                  >  4:246275/1‑138 (MQ=255)
                                                         aaTGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCATGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGc                                                                                                  >  4:163579/1‑138 (MQ=255)
                                                         aaTGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCATGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGc                                                                                                  <  3:163579/138‑1 (MQ=255)
                                                         aaTGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCATGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGc                                                                                                  <  3:246275/138‑1 (MQ=255)
                                                            gcggcgAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAg                                                                                    <  1:57786/149‑1 (MQ=255)
                                                            gcggcgAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAg                                                                                    <  4:129040/149‑1 (MQ=255)
                                                                    cGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATc                                                                            <  2:128255/149‑1 (MQ=255)
                                                                                 gCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTc                                                               >  2:328849/1‑149 (MQ=255)
                                                                                   ttCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGTAAAAATCCAGGTATTTCCTCAc                                                             >  3:347756/1‑149 (MQ=255)
                                                                                   ttCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCAc                                                             <  4:289312/149‑1 (MQ=255)
                                                                                        gcgcTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGtttt                                                        <  2:384874/149‑1 (MQ=255)
                                                                                          gcTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTaa                                                      <  4:9859/149‑1 (MQ=255)
                                                                                            tCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGc                                                    <  2:133665/149‑1 (MQ=255)
                                                                                               tttATCGCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAGCACTGGAAAAGATGTTTACGATGCGCTGACCTGCTGCAAAATAACCAGGAGAAATCCAGGTATTTGCTCACGTTTTAAGCCAg                                                 >  4:51966/1‑149 (MQ=255)
                                                                                                     ccccAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCAcc                                           >  3:183052/1‑149 (MQ=255)
                                                                                                      cccAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCAccc                                          >  4:310599/1‑149 (MQ=255)
                                                                                                         aGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTc                                       >  2:454496/1‑149 (MQ=255)
                                                                                                            aaaaaTTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACt                                    >  1:209500/1‑149 (MQ=255)
                                                                                                               aaTTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTaaa                                 <  3:89493/149‑1 (MQ=255)
                                                                                                                aTTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTaagg                                >  3:113486/1‑147 (MQ=255)
                                                                                                                 ttGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACAGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAgg                               <  4:182713/149‑1 (MQ=255)
                                                                                                                    gTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTAAAGCCAGCTCACCCTTCACTAAAGGGAc                            <  4:347756/149‑1 (MQ=255)
                                                                                                                     ttAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACa                           >  3:478404/1‑149 (MQ=255)
                                                                                                                       aataaCCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACaaa                         >  3:306310/1‑149 (MQ=255)
                                                                                                                        ataaCCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCCCGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAg                        <  3:71972/149‑1 (MQ=255)
                                                                                                                             cAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCAc                   <  2:404677/149‑1 (MQ=255)
                                                                                                                               gTGACATAATTACGGTGCCAGGCACCCTACTGAACACTGGCAAAGATGTTCACGATACGCTGACCTGCCGCAAAATAACCAGGAAAAATCCATGTATTTCCTCACGTTTTACGCCAGCTTACCCTTCACTAAAGGGACAAAGCGCACgc                 >  4:279930/1‑148 (MQ=255)
                                                                                                                               gTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACgg                 >  2:412452/1‑149 (MQ=255)
                                                                                                                                     tAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCa           <  2:3834/149‑1 (MQ=255)
                                                                                                                                      aaTTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCAc          <  4:209968/149‑1 (MQ=255)
                                                                                                                                       aTTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACg         <  4:384655/149‑1 (MQ=255)
                                                                                                                                              tGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGa  <  2:175913/149‑1 (MQ=255)
                                                                                                                                               gCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTaaa                                 <  2:243405/117‑1 (MQ=255)
                                                                                                                                               gCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTaaa                                 >  1:243405/1‑117 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGATTTGAAGTGTCAGAACAGCCTGCCAGCCATAGCGAAACCAGTGACAAAGCCGCAATGCGGCGAACGGTGAATTTTGGGCTTCCCGCGCTCATTTATCCCCCAGGAAAAATTGGTTAATAACCAGTGACATAATTACCGTGCAAGGCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGA  >  NC_000913/2868658‑2868948

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