Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 167,295 T→C 6.1% intergenic (+31/‑189) mrcB → / → fhuA peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,2950TC 6.1% 180.0 / 2.8 66intergenic (+31/‑189)mrcB/fhuApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base T (21/41);  new base C (2/2);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-02

ATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAA  >  NC_000913/167155‑167437
                                                                                                                                            |                                                                                                                                              
aTCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGATCAGCAAGAGCAGAAAGACCGCGACGGCGAAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGc                                                                                                                                              <  1:42553/143‑1 (MQ=255)
  cAGTCTTCGCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCa                                                                                                                                            <  1:213574/143‑1 (MQ=255)
      cgactcAGACGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCg                                                                                                                                        <  1:69188/140‑1 (MQ=255)
       ttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCgg                                                                                                                                       <  1:223805/143‑1 (MQ=255)
       ttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCgg                                                                                                                                       <  2:147125/143‑1 (MQ=255)
            aGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaaa                                                                                                                                  <  1:589236/143‑1 (MQ=255)
              ccGCAGCTATTGCCTCAACCACAACCGGCTCAGAAACAGCAGAAAGACGGCGACGGTGTAGCCGGTGGGATCAAGGAGATGTTGGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATa                                                                                                                                <  1:149653/143‑1 (MQ=255)
                 cagcaACAGCCGCAACAGCAACCTGCTCAGCACGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACa                                                                                                                             >  2:243245/1‑143 (MQ=255)
                    caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATcc                                                                                                                          >  1:904614/1‑143 (MQ=255)
                      aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCa                                                                                                                        >  2:255522/1‑143 (MQ=255)
                      aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGGTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCa                                                                                                                        >  2:873406/1‑143 (MQ=255)
                        aGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGGAAGATAACATCCCATc                                                                                                                      >  1:932690/1‑143 (MQ=255)
                           cGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGGCAGGGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGGAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTaa                                                                                                                   <  2:178876/143‑1 (MQ=255)
                            ggaaCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                  <  1:870296/141‑1 (MQ=255)
                            gCAACAGGAACCTTCTCAGCAAGAGCAGAAAGGCATCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                  <  2:205873/143‑1 (MQ=255)
                             cAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGa                                                                                                                 <  1:550964/143‑1 (MQ=255)
                              aaCAGCAACCTGCGCAGCAAGAGGAGAAAGACAGCGACGGTGTAGTCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGat                                                                                                                <  2:287430/143‑1 (MQ=255)
                                cAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGatat                                                                                                              >  1:34042/1‑143 (MQ=255)
                                      caagCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAAccc                                                                                                        <  1:939810/140‑1 (MQ=255)
                                         gCTCGGCGAGAGCGGAGAGAGAGCGAGGGTGTGGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCcttt                                                                                                     <  2:493592/143‑1 (MQ=255)
                                          cTCAGCACGAGCAGAAAGACAGCGCCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttc                                                                                                    <  1:608543/143‑1 (MQ=255)
                                            cAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATCTTAACCctgtctt                                                                                                  >  2:1111644/1‑143 (MQ=255)
                                                aaGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCa                                                                                              <  1:255522/143‑1 (MQ=255)
                                                     cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTgg                                                                                         >  2:613195/1‑143 (MQ=255)
                                                      aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgt                                                                                        >  1:178317/1‑143 (MQ=255)
                                                        aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttg                                                                                      >  1:850296/1‑143 (MQ=255)
                                                         aaGACAGCGACGATGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTCATTACCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAATTAACATCCCATCTACGATATAAACCCTTTCTTTTCATCTGgttgt                                                                                     <  2:990527/143‑1 (MQ=255)
                                                          aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                    <  1:177986/143‑1 (MQ=255)
                                                          aGACAGCGACGGGGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTATCCCTTTCTTTTCATCTGgttgtt                                                                                    <  1:317362/143‑1 (MQ=255)
                                                            aCAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                  >  1:290025/1‑143 (MQ=255)
                                                              aGCGACGGTGTAGCGGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                <  2:342642/143‑1 (MQ=255)
                                                               gcgtttgtgtAGCCGTTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatta                                                                               <  1:499810/137‑1 (MQ=255)
                                                                cGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                              <  1:123675/143‑1 (MQ=255)
                                                                cGAAGGTGGAGCTGGTTGGATCAAGGAGATGGTGGGTAGTAATTAACATCTGAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCGCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                              <  2:78438/143‑1 (MQ=255)
                                                                  aCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAAcc                                                                            <  1:894520/143‑1 (MQ=255)
                                                                      tgtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAATATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                        <  2:738821/143‑1 (MQ=255)
                                                                           ccGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGaa                                                                   <  2:813892/143‑1 (MQ=255)
                                                                            cGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGATTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAAc                                                                  <  2:17199/143‑1 (MQ=255)
                                                                              gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGGTGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCATTCAGGAACGc                                                                >  1:214091/1‑143 (MQ=255)
                                                                                   aTCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGa                                                           >  2:106810/1‑143 (MQ=255)
                                                                                   aTCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTATTCATCTGGTTGCTTATTAACCCTTCAGGAACGCTCAGa                                                           >  2:772193/1‑143 (MQ=255)
                                                                                    tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGACCGCTCAGAt                                                          >  2:574644/1‑143 (MQ=255)
                                                                                        ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCg                                                      >  1:282838/1‑143 (MQ=255)
                                                                                        ggATATGTTGGGGAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTGCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCg                                                      <  1:915500/143‑1 (MQ=255)
                                                                                           tatGTTTGGTAGTAATTAACATCTAAGCGTGAAACACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTAc                                                   <  1:364480/143‑1 (MQ=255)
                                                                                             tGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCg                                                 <  2:339991/143‑1 (MQ=255)
                                                                                             tGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCg                                                 <  2:34042/143‑1 (MQ=255)
                                                                                             tGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGCTGTTTATTATCCCTTCAGGACCGCTCAGCTTGCCTACCg                                                 >  2:1117153/1‑143 (MQ=255)
                                                                                                ttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCGATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCtt                                              <  2:498199/143‑1 (MQ=255)
                                                                                                 tGGTAGTAATTAACATCTAAGCGTGACCTACCGGATGGCCATTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTg                                             <  2:124713/143‑1 (MQ=255)
                                                                                                     agtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGaa                                         >  2:765461/1‑143 (MQ=255)
                                                                                                     agtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGaa                                         <  2:464012/143‑1 (MQ=255)
                                                                                                      gtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAc                                        <  2:942835/143‑1 (MQ=255)
                                                                                                              cATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGc                                >  1:409115/1‑143 (MQ=255)
                                                                                                                 cTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGtt                             >  1:678634/1‑143 (MQ=255)
                                                                                                                  tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCTCCTGCGttt                            <  2:1111338/143‑1 (MQ=255)
                                                                                                                  tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTCTTGCTCTGGTTGTTTCTTAACCCTGCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                            >  2:237851/1‑143 (MQ=255)
                                                                                                                   aaGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                           <  1:1032006/143‑1 (MQ=255)
                                                                                                                      cGTGAAATACCGGATGGCGAGTTGCCACCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTATCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGaa                        <  1:588119/143‑1 (MQ=255)
                                                                                                                          aatgaCCGGATGGCGAGTTGCGATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAAtatt                    <  1:326029/139‑1 (MQ=255)
                                                                                                                          aaaTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCGTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAAtatt                    <  2:521534/143‑1 (MQ=255)
                                                                                                                             tACCGGATGGCGATTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGACCGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                 <  2:923088/143‑1 (MQ=255)
                                                                                                                             caCCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTCTCTTTTCATCTGGTTGGTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                 <  1:449506/142‑1 (MQ=255)
                                                                                                                                cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGCACGCTCTGCTTGCGTACCGCTTGCTAACCCTCCAACGTAACTAATAttatctta              >  1:34390/1‑143 (MQ=255)
                                                                                                                                cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttctctta              >  1:145367/1‑143 (MQ=255)
                                                                                                                                      gCGAGTTGCCATCCGTTAAAATATCATCACATCTAAGATATTAACCATTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTta        <  2:99753/143‑1 (MQ=255)
                                                                                                                                            tGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGTAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTAtaataa  <  1:884921/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
ATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAA  >  NC_000913/167155‑167437

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