Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 566,620 C→T 6.6% intergenic (+32/+41) bfr‑I → / ← uvrA bacterioferritin/uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947566,6200CT 6.6% 175.3 / 2.1 60intergenic (+32/+41)bfr‑I/uvrAbacterioferritin/uvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/2);  total (31/29)
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.27e-01

ACATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTT  >  NC_002947/566473‑566768
                                                                                                                                                   |                                                                                                                                                    
aCATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGc                                                                                                                                                      <  7:97482/148‑1 (MQ=255)
   tCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCa                                                                                                                                                  >  2:107666/1‑149 (MQ=255)
     cTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCAt                                                                                                                                                 <  2:300905/148‑1 (MQ=255)
     cTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACAACTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCAACTGACTCATTTAAAAGCCCATGGCACCATc                                                                                                                                                >  3:234047/1‑149 (MQ=255)
        cgtgcgCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCg                                                                                                                                               <  5:320958/147‑1 (MQ=255)
           gcgcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAc                                                                                                                                           >  7:146381/1‑148 (MQ=255)
            ctcAGCTGGCCGAAACCGAATAAGATCACACCTACTGGCTTGAGAAGCAGCAGGTCCTGATCAATGCCATCGGGCTGTAAAACTACCTTCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAcc                                                                                                                                          <  4:234047/146‑1 (MQ=255)
             gcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAg                                                                                                                                        <  7:32493/149‑1 (MQ=255)
                  tGGCCGACACCGAAGAAGCTCACACCTACTCGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAc                                                                                                                                           >  8:7688/1‑141 (MQ=255)
                  tGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAc                                                                                                                                           <  7:7688/141‑1 (MQ=255)
                     ccGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCtttt                                                                                                                                >  8:193082/1‑149 (MQ=255)
                        acacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCa                                                                                                                             <  1:335254/149‑1 (MQ=255)
                          acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCCCCAGCGGTCCCAGGGGCTTTTTCAtt                                                                                                                           >  8:295270/1‑149 (MQ=255)
                           ccGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCGTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAATTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAACAGGGGCTTTTTCAtt                                                                                                                           <  2:342703/148‑1 (MQ=255)
                           ccGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTg                                                                                                                          <  5:349519/149‑1 (MQ=255)
                                     tCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAAt                                                                                                                <  8:146381/149‑1 (MQ=255)
                                           cTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAg                                                                                                          >  6:163116/1‑149 (MQ=255)
                                            tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCgg                                                                                                                                              <  7:45100/112‑1 (MQ=255)
                                            tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCgg                                                                                                                                              >  8:45100/1‑112 (MQ=255)
                                            tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAgg                                                                                                         >  5:331297/1‑149 (MQ=255)
                                               tGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGccc                                                                                                      >  5:14490/1‑149 (MQ=255)
                                                 gCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTa                                                                                                    >  4:202247/1‑149 (MQ=255)
                                                    tGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAgg                                                                                                                                       <  1:317057/111‑1 (MQ=39)
                                                    tGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAgg                                                                                                                                       >  2:317057/1‑111 (MQ=39)
                                                    tGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATcc                                                                                                 >  6:113108/1‑149 (MQ=255)
                                                       gaAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCg                                                                                               >  4:334872/1‑148 (MQ=255)
                                                       gaAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGc                                                                                              <  5:128505/149‑1 (MQ=255)
                                                                   ccTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgg                                                                                                                                       >  8:315602/1‑96 (MQ=255)
                                                                   ccTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgg                                                                                                                                       <  7:315602/96‑1 (MQ=255)
                                                                     tGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGtt                                                                                 >  8:46339/1‑148 (MQ=255)
                                                                       gatcaGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                               <  3:186970/144‑1 (MQ=255)
                                                                               cATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGAACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTc                                                                                            <  1:48312/127‑1 (MQ=255)
                                                                               cATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGAACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTc                                                                                            >  2:48312/1‑127 (MQ=255)
                                                                                   gggCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCgg                                                                  <  6:331297/149‑1 (MQ=255)
                                                                                   gggCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCgg                                                                  <  5:163116/149‑1 (MQ=255)
                                                                                   gggCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                   <  6:14490/148‑1 (MQ=255)
                                                                                         gAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                             <  1:334468/148‑1 (MQ=255)
                                                                                                     cAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCAtct                                                <  7:295270/149‑1 (MQ=255)
                                                                                                       gTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCg                                              <  7:215270/149‑1 (MQ=255)
                                                                                                       gTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATAGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCg                                              >  1:108039/1‑149 (MQ=255)
                                                                                                          gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCt                                            >  1:148135/1‑148 (MQ=255)
                                                                                                                gTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                             >  6:12655/1‑109 (MQ=255)
                                                                                                                gTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                             <  5:12655/109‑1 (MQ=255)
                                                                                                                   aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttc                                   <  2:148135/148‑1 (MQ=255)
                                                                                                                           acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGCCTGCTTCATGTCGCTCAGCTCTTCTGGTGTAcc                          >  2:245951/1‑149 (MQ=255)
                                                                                                                             aGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGc                        <  5:113108/149‑1 (MQ=255)
                                                                                                                                cTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAgg                     <  3:202247/149‑1 (MQ=255)
                                                                                                                                   atgaaaaaGCCCCTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGtata                                                              <  6:256324/101‑1 (MQ=37)
                                                                                                                                   atgaaaaaGCCCCTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGtata                                                              >  5:256324/5‑105 (MQ=37)
                                                                                                                                   aaGTAAAAGCCCCTGGCACCATCGGTACCAGGTGCTTTTTCAt                                                                                                                            <  4:333008/43‑1 (MQ=255)
                                                                                                                                   aaGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAt                                                                                                                            >  3:333008/1‑43 (MQ=255)
                                                                                                                                    aGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAAt                 >  6:234367/1‑149 (MQ=255)
                                                                                                                                           gCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGc          <  7:46339/149‑1 (MQ=255)
                                                                                                                                            ccccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGcc         >  4:310190/1‑149 (MQ=255)
                                                                                                                                               cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAc      >  4:37750/1‑149 (MQ=255)
                                                                                                                                               cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAc      >  8:60669/1‑149 (MQ=255)
                                                                                                                                                tGGCACCATCGGTACCAGGTGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAcc     >  6:39855/1‑149 (MQ=255)
                                                                                                                                                tGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTACCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                 >  3:109139/1‑121 (MQ=255)
                                                                                                                                                tGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTACCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                 <  4:109139/121‑1 (MQ=255)
                                                                                                                                                 ggCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAccc    >  5:146270/1‑149 (MQ=255)
                                                                                                                                                  gCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCt   >  1:237628/1‑149 (MQ=255)
                                                                                                                                                   cACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCtt  <  2:237628/149‑1 (MQ=255)
                                                                                                                                                   cACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCtt  <  3:334872/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTT  >  NC_002947/566473‑566768

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