Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,204,350 A→G Q344R (CAG→CGG)  zraR → fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,204,3500AG85.2% 142.5 / ‑0.6 54Q344R (CAG→CGG) zraRfused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/23);  minor base C (8/0);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATCC  >  NC_000913/4204201‑4204493
                                                                                                                                                          |                                                                                                                                               
tCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCa                                                                                                                                                 <  2:862367/150‑1 (MQ=255)
  gCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGc                                                                                                                                              <  2:831023/151‑1 (MQ=255)
    cAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATTCCT‑CTGCTGGCTGGCCATTTTTTGCAGCGCTTT‑GCCGAGCGGAATCGC‑CA‑GG‑CGGGAAAAAGGTTTATGCCCCCGGCg                                                                                                                                             >  1:30519/1‑150 (MQ=255)
      ggATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATg                                                                                                                                          <  1:487561/151‑1 (MQ=255)
          ctctACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGAAATTCCT‑CTGCTTGCTGGCCATTTTTTGCAGCGCCTT‑TCCGAGGGGAATCGC‑AAGGG‑GGGGAAAAGGTTTTAGCCCCCGGGGGAgggt                                                                                                                                        >  2:552935/1‑149 (MQ=255)
            ctactaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGa                                                                                                                                        <  1:703285/147‑1 (MQ=255)
            ctactaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATc                                                                                                                                      <  2:22862/149‑1 (MQ=255)
            ctactaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGAAATCCCC‑CTTCTGGGCGGCCACTTTTTTCCGCGCTTT‑TGCGGGGGGTATCTC‑AA‑GGAGGGGAAAAGGTTTT‑CTCCCCCGGCg                                                                                                                                             >  2:1001042/1‑142 (MQ=255)
              actaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGGTGGCCATTTTTTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTggt                                                                                                                                  >  1:440416/1‑149 (MQ=255)
              actaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCt                                                                                                                                     <  1:1393/148‑1 (MQ=255)
                  tCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGGTGGGCCTTTTCTGCAGCGCGTT‑TCCGCGCGGAAACCC‑AA‑AG‑CGGGAAAAGGGTTTACGCCCCGGGCGATGGGTCTGTTGAtt                                                                                                                              >  2:251653/1‑151 (MQ=255)
                   cGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTc                                                                                                                             <  1:271312/151‑1 (MQ=255)
                    gCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCa                                                                                                                            <  1:40073/151‑1 (MQ=255)
                     ccTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGGTATTCCT‑CTGCTGGGTGGCCCATTTTTGCCGCGCCTT‑TCCGGGCGGAATCGC‑AA‑AG‑CGGGCAAAGGGTTTTCGCCCCCGGCGGTGGGTCTGGTGATTCAt                                                                                                                           >  1:763389/1‑151 (MQ=255)
                        gAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTa                                                                                                                         <  1:124266/150‑1 (MQ=255)
                        gAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCC‑CTGCTGGGTGGCCATTTTCTTCAGCGCTTT‑TGCGAGCGTAATCGC‑CA‑GG‑CGGTAAAAGGGTTTTCGCCCCCGGCGATGGAGCTGTTGTTTCAtt                                                                                                                          >  2:854049/1‑149 (MQ=255)
                        gAATGTGGTGGCGATTGAAGTACCATCGCTGCGGAAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTAc                                                                                                                        <  1:141584/151‑1 (MQ=255)
                            gtggtgGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGAAATCACC‑CCGCCGGCTGGCCATTTTTTTCAAGGCTTT‑TGCCAGCGTAAACGCAAA‑GG‑GGGAAAAAGGTTTTA‑GCCCCCGGGGAGGGGTCTGTTTGTTCAtt                                                                                                                          >  1:377898/1‑145 (MQ=255)
                              ggtggCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑AG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTg                                                                                                                   <  1:985678/150‑1 (MQ=255)
                                  gCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATTCCC‑CCGCTTGCTGGGCCCTTTTTTCCGCGCTTTGGCCGAG‑GGAATCCC‑AA‑AG‑CGGGAAAAGGGTTTAAGCCCCGGGCGGTGGGTCTTTTGTTTTAATTCGACTTGCCgg                                                                                                              >  2:552569/1‑151 (MQ=255)
                                    gATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGa                                                                                                            <  1:303584/151‑1 (MQ=255)
                                      ttGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCCTTTTTTTCACCGCTTT‑GCCGAGCGTAATCGC‑AA‑AG‑CGGGAAAAGGTTTTTCGCCCCCGGCGGTGGGTCTGTTGTTTCATTTAGACTGGCCgggga                                                                                                           >  1:605954/1‑148 (MQ=255)
                                          aGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATAtt                                                                                                      <  2:841144/151‑1 (MQ=255)
                                          aGTACCATCGCTGCGGCAACGGCGGAAAGATATACCTCCTGC‑GGCTGGCCATTCTCTTCTGCACCTT‑GGCCACCCCAATCCG‑GA‑CG‑CGGGGAAAGGATTTTCGCCCCCGGCGGTGGTTCTTTTTGTTTATTACGAcacgcc                                                                                                                >  2:384046/1‑136 (MQ=255)
                                               cATCGCTGCGGCAACGGCGGGAAGACATCCCT‑CCGCTTGCTGGCCCTTTTTTTCACCGCTTT‑TCCGAGCGGAAACCC‑AA‑GG‑CGGGAAAAGGGTTTATGCCCCGGGGGATGGGTCTGTTTATTCATTACGACTGGCCGGGGAATATTCGTg                                                                                                  >  1:430443/1‑150 (MQ=255)
                                                aTCGCTGCGGCAACGGCGGGGAGAGAACCCT‑CTTCCGGCTGGCCCATTTTTGCAGCGGCTT‑GCCGAGAGGAAACGC‑AA‑AG‑GGGGAAAAGGGTTTACGCCCCGGGCGGTGGGGCTGTTGATTCATTTAGAACGGCCGGGGAATATTtgtggg                                                                                                >  2:96128/1‑149 (MQ=255)
                                                   gCTGCGGCAACGGCGGGAAGAAATCCCC‑CTCCAGGCCGGCCCTTTTCTGCCGACCTCT‑GTCCGGCGCAAAACC‑AA‑AG‑CGGGAAAAAGTTTTTAGCCCCCGGCGAGGGGTCTGTTGTTTCATTACTAACGGCCCGGAAAGATTTGGGAGc                                                                                               >  2:972412/1‑149 (MQ=255)
                                                     tGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGGTGGCCAATTTTTGCAGCGCTTT‑GCCGGGCGGAAACGC‑AA‑AG‑CGGGAAAAGGGTTTACGCCCCGGGCGATGGATCTGTTGGTTTATTTAGACTGGGCGGGGAATAATTGTGGGCTgg                                                                                            >  2:358393/1‑150 (MQ=255)
                                                       cGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGAc                                                                                                                     <  1:414407/123‑1 (MQ=255)
                                                        ggCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGaaaa                                                                                        <  1:131787/151‑1 (MQ=255)
                                                             cggcggGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCgg                                                                                   >  1:94605/1‑151 (MQ=255)
                                                                tgggAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTgg                                                                                >  2:116350/2‑151 (MQ=255)
                                                                       atatCCCT‑CTCCTGGCTGGCCATTTTCTGCACCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGCAATATTCGTGAGCTGGAAACCGCGGTGGAACggg                                                                          >  2:1049548/1‑150 (MQ=255)
                                                                        tatCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGc                                                                         >  2:968249/1‑150 (MQ=255)
                                                                         atCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAg                                                                       <  1:879052/151‑1 (MQ=255)
                                                                                     ggctggcCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTgctgct                                                             <  1:738652/150‑1 (MQ=255)
                                                                                     ggctggcCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGctgctg                                                            >  1:958418/1‑151 (MQ=255)
                                                                                       ctggcCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGAc                                                          >  1:941393/1‑151 (MQ=255)
                                                                                         agcgatttgatgaaGCGCGTT‑GCCGAGCGGAATCAC‑AA‑GG‑CGGGAAACGGCTTTACGCCCAGGGCGATGGATCTGCTGATTCACTACGACTGGCCGGGAAATATTCGCGAGCTGGAAAACGCGGTGGAATGGACAGTGGGGCTGCTGACCg                                                        >  2:1048366/14‑151 (MQ=255)
                                                                                            cATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCgggg                                                     <  1:94445/151‑1 (MQ=255)
                                                                                               tttCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAat                                                  >  1:233553/1‑151 (MQ=255)
                                                                                                ttCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAat                                                  <  2:263097/150‑1 (MQ=255)
                                                                                                  cTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatat                                                >  2:232788/1‑150 (MQ=255)
                                                                                                   tGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatatat                                              >  1:984900/1‑151 (MQ=255)
                                                                                                               gCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAcgc                                    <  2:877814/150‑1 (MQ=255)
                                                                                                                gcGAGCGTAATCGC‑AA‑GG‑CGGTGAAAGGGTTTACGCCCCGGGCGATGGATCTGTGGTTTCATTACGACTGGCAGGGAAATATTCGTGAGCTGGAAAAACCGGTGGAACGGGCAGTGGGGCTGCTGACCGGGGAATATATTTCCGAAcgc                                    >  2:217603/2‑149 (MQ=255)
                                                                                                                ccGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGa                                  <  1:461906/151‑1 (MQ=255)
                                                                                                                    gCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCt                               <  2:457859/150‑1 (MQ=255)
                                                                                                                    gCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTg                              <  1:700773/151‑1 (MQ=255)
                                                                                                                       tAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCg                           >  2:560354/1‑151 (MQ=255)
                                                                                                                         aTCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCt                         <  1:25770/151‑1 (MQ=255)
                                                                                                                             c‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCg                     >  1:738910/1‑151 (MQ=255)
                                                                                                                                         aaaaGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGt            >  1:552816/1‑151 (MQ=255)
                                                                                                                                                  tACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATc   >  2:459923/1‑151 (MQ=255)
                                                                                                                                                   aCGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATcc  >  1:429233/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                               
TCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCT‑CTGCTGGCTGGCCATTTTCTGCAGCGCTTT‑GCCGAGCGTAATCGC‑AA‑GG‑CGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATCC  >  NC_000913/4204201‑4204493

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