Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,820,992 T→G 5.9% M297L (ATG→CTG)  ptxD ← phosphonate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,820,9920TG 5.9% 196.0 / 2.2 62M297L (ATG→CTG) ptxDphosphonate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base G (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTGGCCCGCAGCGCTTCG  >  NC_002947/3820844‑3821140
                                                                                                                                                     |                                                                                                                                                    
cGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCAt                                                                                                                                                      >  7:239043/1‑149 (MQ=255)
        cGCAAGGCAGGCAACTGCCGCTGACCGTGTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTgg                                                                                                                                               >  1:182609/1‑149 (MQ=255)
         gCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGc                                                                                                                                              <  6:114332/148‑1 (MQ=255)
          cAAGGCAGGCAACTGCCGCTGACCGTTTGAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCgg                                                                                                                                            >  1:182956/1‑149 (MQ=255)
          cAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCgg                                                                                                                                            <  8:239043/149‑1 (MQ=255)
            aggcaggcAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGt                                                                                                                                          >  1:98627/1‑149 (MQ=255)
             ggcaggcaACTGGCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGGGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGtt                                                                                                                                         <  2:182609/149‑1 (MQ=255)
                   caACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCg                                                                                                                                   <  5:115072/149‑1 (MQ=255)
                   caACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGACCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCg                                                                                                                                   <  7:60097/149‑1 (MQ=255)
                         ccGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACAAGGTTTACCTGCCGCTCACCGGCCATCGCGCTGAGCAGGTTGTCCACCGCGCAAAGCGCCATGGCCTGGCGGCTTTCCTCGGTCgcc                                                                                                                             >  8:219693/1‑149 (MQ=255)
                                ccGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCGCCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCgccgagccga                                                                                                                      >  4:7028/1‑149 (MQ=255)
                                 cGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCgccgagccga                                                                                                                      >  1:211708/1‑148 (MQ=255)
                                 cGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCgccgagccga                                                                                                                      >  4:304315/1‑148 (MQ=255)
                                                 gcgGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTAc                                                                                                      <  2:98627/149‑1 (MQ=255)
                                                   ggCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCa                                                                                                    <  2:320030/149‑1 (MQ=255)
                                                    gCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCAc                                                                                                   >  5:186014/1‑149 (MQ=255)
                                                    gCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGGCGCTACCAc                                                                                                   <  3:189833/149‑1 (MQ=255)
                                                     ccccTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACg                                                                                                  <  2:182956/149‑1 (MQ=255)
                                                            gggAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATc                                                                                             >  4:239168/1‑147 (MQ=255)
                                                            gggAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATc                                                                                             <  3:239168/147‑1 (MQ=255)
                                                                      gcgcGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAg                                                                                 >  3:27386/1‑149 (MQ=255)
                                                                          gTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAgcggc                                                                             <  2:211708/149‑1 (MQ=255)
                                                                            cAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCgg                                                                                                               >  2:208416/1‑113 (MQ=255)
                                                                            cAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCgg                                                                                                               <  1:208416/113‑1 (MQ=255)
                                                                             aGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGa                                                                           <  3:215905/148‑1 (MQ=255)
                                                                             aGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGGCGCTACCACGTTATCCAGCTGCAACAGCGGCGAg                                                                          <  8:47315/149‑1 (MQ=255)
                                                                             aGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCAGGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAg                                                                          <  8:43229/149‑1 (MQ=255)
                                                                              gCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGt                                                                         <  3:304315/149‑1 (MQ=255)
                                                                                 gTTGACCAGTTTTACCGTCCGCTCACCGGCCAGTTCGATGAGCAGGTTGTCAACCGCTCAACGCGCCATGGCCTGGCTGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGa                                                                      <  7:219693/149‑1 (MQ=255)
                                                                                     ctgaGGTTTAACGGGCTGTCACGGGCCTGCGCGCTGAGCTGGTTGTCAACCGCTCAACGCGCCATGGCCTGGTGGGTTTGCTCGGTCGGCGAGCGGATATGGGGGGTCGCTACCAGTTTATCCTGCTGCTTGAGCGGGGAGTCGATCg                                                                   <  4:38795/145‑1 (MQ=255)
                                                                                     aCCAGGTTTACCGGCCGCTCACAGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCgg                                                                  >  4:311277/1‑149 (MQ=255)
                                                                                      ccAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCt                                                                                        >  5:52131/1‑126 (MQ=255)
                                                                                      ccAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCt                                                                                        <  6:52131/126‑1 (MQ=255)
                                                                                        aGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGCGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCg                                                                   >  3:38795/1‑145 (MQ=255)
                                                                                               ccGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCg                                                        >  1:157745/1‑149 (MQ=255)
                                                                                                cGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAgc                                                                                >  6:63297/1‑124 (MQ=255)
                                                                                                cGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAgc                                                                                <  5:63297/124‑1 (MQ=255)
                                                                                                   ccGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCGCCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCGACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCa                                                    <  3:7028/149‑1 (MQ=255)
                                                                                                     gCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACg                                                  <  2:157745/149‑1 (MQ=255)
                                                                                                           cGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACg                                            >  4:78423/1‑149 (MQ=255)
                                                                                                           cGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACg                                            >  7:45491/1‑149 (MQ=255)
                                                                                                             gCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTc                                          <  3:78423/149‑1 (MQ=255)
                                                                                                              ccAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGTATCGTGCACGAACACGTcc                                         <  4:27386/149‑1 (MQ=255)
                                                                                                               cAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATCCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCa                                        >  4:245225/1‑149 (MQ=255)
                                                                                                                   gcgcTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGccc                                    >  2:174850/1‑149 (MQ=255)
                                                                                                                     gcTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCg                                  >  3:142277/1‑149 (MQ=255)
                                                                                                                       tGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCa                                >  3:25709/1‑149 (MQ=255)
                                                                                                                              gTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGa                         >  4:141469/1‑149 (MQ=255)
                                                                                                                              gTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGa                         <  1:174850/149‑1 (MQ=255)
                                                                                                                               ttGTCCTCCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTAGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGAt                        <  1:233934/149‑1 (MQ=255)
                                                                                                                                  tCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATgc                      <  6:186014/148‑1 (MQ=255)
                                                                                                                                       gaggaaaCCCGCCAGGCCATGGCGCGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTgg                >  2:18945/6‑149 (MQ=255)
                                                                                                                                       cgcgCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCt                  >  6:124077/1‑147 (MQ=255)
                                                                                                                                       cgcgCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCt                  <  5:124077/147‑1 (MQ=255)
                                                                                                                                        gcgcAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGt                                                       <  4:192040/109‑1 (MQ=255)
                                                                                                                                        gcgcAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGt                                                       >  3:192040/1‑109 (MQ=255)
                                                                                                                                         ctcagcgcgCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTGGc               <  6:56583/143‑1 (MQ=255)
                                                                                                                                          gcAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATgcgc                    >  6:285288/1‑142 (MQ=255)
                                                                                                                                          gcAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATgcgc                    <  5:285288/142‑1 (MQ=255)
                                                                                                                                            aaCCCGCCAGGCCATGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACcgc                            <  6:180918/132‑1 (MQ=255)
                                                                                                                                            aaCCCGCCAGGCCATGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACcgc                            >  5:180918/1‑132 (MQ=255)
                                                                                                                                               gcgcCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTGGCCCGCAgc        <  3:311277/149‑1 (MQ=255)
                                                                                                                                                cgcCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTGGCCCGCAgcg       <  4:25709/149‑1 (MQ=255)
                                                                                                                                                     tGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTGGCCCGCAGCGCTTCg  <  7:20747/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTCGATCGGCAGTGGCTCGTGCACGAACACGTCCAGCCCCGCACCGCGGATGCGCCTGGCCCGCAGCGCTTCG  >  NC_002947/3820844‑3821140

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