Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,654,891 T→C 7.1% intergenic (+43/‑215) aceA → / → PP_4118 isocitrate lyase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,654,8910TC 7.1% 163.6 / 2.5 56intergenic (+43/‑215)aceA/PP_4118isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base C (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGCC  >  NC_002947/4654743‑4655039
                                                                                                                                                     |                                                                                                                                                    
ggAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCt                                                                                                                                                      >  7:223910/1‑149 (MQ=255)
       ggCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGc                                                                                                                                                <  1:207788/148‑1 (MQ=255)
         cACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTAAAAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCt                                                                                                                                              <  6:110673/149‑1 (MQ=37)
         cACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCtt                                                                                                                                             >  2:277751/1‑149 (MQ=255)
              gCTACTTCGACGACATGACCACCGTGATCCAGGGGGGCGGTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                         <  2:103664/148‑1 (MQ=255)
              gCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCGCGTGGCTTGTGAAAGCTCCGGGCCTTTGttt                                                                                                                                        <  4:40916/149‑1 (MQ=255)
              gCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGttt                                                                                                                                        >  5:124253/1‑149 (MQ=255)
              gCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGttt                                                                                                                                        <  6:52969/149‑1 (MQ=255)
                    tCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAccc                                                                                                                                  >  4:221042/1‑149 (MQ=255)
                     cgacgaCATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAccc                                                                                                                                  >  1:16156/1‑148 (MQ=255)
                      gacgacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCga                                                                                                                                <  3:52698/149‑1 (MQ=255)
                      gacgacATGACCACCGTATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAccc                                                                                                                                  >  7:116980/1‑146 (MQ=255)
                       acgacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCga                                                                                                                                <  8:223910/148‑1 (MQ=255)
                           cATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAg                                                                                                                           <  4:118932/149‑1 (MQ=255)
                             tGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAggg                                                                                                                         <  1:277751/149‑1 (MQ=255)
                                   ccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgggc                                                                                                                   <  3:221042/149‑1 (MQ=255)
                                      tGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctg                                                                                                                <  7:261060/149‑1 (MQ=255)
                                        aTCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCTGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgct                                                                                                              >  5:35823/1‑149 (MQ=255)
                                         tCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgctg                                                                                                             <  6:124253/149‑1 (MQ=255)
                                                    gcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTa                                                                                                  <  8:116980/149‑1 (MQ=255)
                                                    gcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTa                                                                                                  <  2:16156/149‑1 (MQ=255)
                                                                aCTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAc                                                                                      >  2:19554/1‑149 (MQ=255)
                                                                  tGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCg                                                                                    >  8:81360/1‑149 (MQ=255)
                                                                       tGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTcc                                                                               >  4:13033/1‑149 (MQ=255)
                                                                                 gACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTAcc                                                                      <  7:81360/148‑1 (MQ=255)
                                                                                    cGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAt                                                                  >  7:25663/1‑149 (MQ=255)
                                                                                     gaagaaGAGCAGTTCCACTGATAGGGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAtt                                                                 <  5:165211/149‑1 (MQ=255)
                                                                                      aagaagAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTc                                                                >  4:240295/1‑149 (MQ=255)
                                                                                             gCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACACGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGt                                                          >  5:198006/1‑148 (MQ=255)
                                                                                                gTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccgc                                                                             <  4:88619/126‑1 (MQ=255)
                                                                                                gTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccgc                                                                             >  3:88619/1‑126 (MQ=255)
                                                                                                    cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTTGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                  <  8:113847/149‑1 (MQ=255)
                                                                                                    cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                   >  6:180494/1‑148 (MQ=255)
                                                                                                    cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                  <  4:108520/149‑1 (MQ=255)
                                                                                                       tGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCAcc                                               >  7:127658/1‑149 (MQ=255)
                                                                                                         aTAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTa                                             <  3:13033/149‑1 (MQ=255)
                                                                                                          tAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCTGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTac                                            <  6:35823/149‑1 (MQ=255)
                                                                                                          tAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTac                                            <  8:25663/149‑1 (MQ=255)
                                                                                                           aGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTaca                                           >  7:29227/1‑149 (MQ=255)
                                                                                                            gTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTacac                                          <  3:239741/149‑1 (MQ=255)
                                                                                                                gatGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                  <  6:33878/135‑1 (MQ=255)
                                                                                                                    ccGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCggg                                  >  6:175129/1‑149 (MQ=255)
                                                                                                                            aCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                            >  7:114398/1‑148 (MQ=255)
                                                                                                                              aaaGGCCCGGGGCTTGTGAAAGCTCCGAGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                          >  8:25392/1‑148 (MQ=255)
                                                                                                                              aaaGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                          >  8:78269/1‑148 (MQ=21)
                                                                                                                                aGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGacaca                        >  2:289502/1‑148 (MQ=255)
                                                                                                                                  gCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACAt                     >  5:41125/1‑149 (MQ=255)
                                                                                                                                  gCCCGGGGCTTGTGAAAGCTCCGAGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACAt                     <  7:25392/149‑1 (MQ=255)
                                                                                                                                    ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCg                   >  7:211450/1‑149 (MQ=255)
                                                                                                                                    ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCg                   >  3:268864/1‑149 (MQ=255)
                                                                                                                                     cGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg                  <  6:198006/149‑1 (MQ=255)
                                                                                                                                        ggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGa               <  8:114398/149‑1 (MQ=255)
                                                                                                                                          ctttcacaaGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGAtt             >  3:81322/8‑149 (MQ=255)
                                                                                                                                             tcacaaGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTgg           <  7:78269/144‑1 (MQ=255)
                                                                                                                                             gtgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAAcc                                                     <  3:152058/106‑1 (MQ=255)
                                                                                                                                             gtgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAAcc                                                     >  4:152058/1‑106 (MQ=255)
                                                                                                                                             gtgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTgg           <  3:240295/148‑1 (MQ=255)
                                                                                                                                                     tCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGcc  >  1:170992/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGCC  >  NC_002947/4654743‑4655039

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: