Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,614 G→C 21.4% intergenic (‑172/+46) pvdH ← / ← PP_4224 diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,796,6140GC21.4% 122.9 / 29.3 56intergenic (‑172/+46)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (15/29);  new base C (7/5);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

CTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGC  >  NC_002947/4796468‑4796759
                                                                                                                                                  |                                                                                                                                                     
cTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTa                                                                                                                                                     <  4:151828/149‑1 (MQ=255)
 tGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAAcccgcccg                                                                                                                                                       >  6:134979/1‑146 (MQ=255)
   gAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTa                                                                                                                                                     <  5:134979/146‑1 (MQ=255)
   gAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAg                                                                                                                                                    <  5:70503/147‑1 (MQ=255)
   gAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAg                                                                                                                                                    >  6:70503/1‑147 (MQ=255)
      gCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCgg                                                                                                                                              <  2:9046/149‑1 (MQ=255)
       cGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTGGAGTGGGGGTAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCgg                                                                                                                                              <  7:163644/149‑1 (MQ=255)
       cGGTGGAAATGGGTAAATCCTCGCGATCAAGCGGATGGGCTGTTATCGCTTTGTAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCgg                                                                                                                                              <  3:229188/149‑1 (MQ=255)
        ggTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTAcgggcggg                                                                                                                                             <  3:263998/149‑1 (MQ=38)
        ggTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCggg                                                                                                                                             <  6:312217/149‑1 (MQ=255)
           ggACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTc                                                                                                                                          <  4:266032/149‑1 (MQ=255)
             aattGGGTTATTCCTCGCGAGCAAGTGGATGTGCGGTTATGGCTTTGGTTACGCTTGAGTGGGGGGAGGATTAGCGGTTGTTTGGCGATTTTTTTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGTGCGTGTTCAc                                                                                                                                        <  6:26491/145‑1 (MQ=37)
              cATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAcc                                                                                                                                       <  3:102097/149‑1 (MQ=255)
               aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                      >  4:237659/1‑149 (MQ=38)
               aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                      <  3:237659/149‑1 (MQ=38)
               aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                                      <  3:64820/149‑1 (MQ=255)
                 gggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCggg                                                                                                                                             >  4:83627/1‑140 (MQ=255)
                 gggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCggg                                                                                                                                             <  3:83627/140‑1 (MQ=255)
                      aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGaa                                                                                                                                 >  4:9254/1‑147 (MQ=255)
                      aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               <  8:209238/149‑1 (MQ=255)
                          cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                            <  8:25034/148‑1 (MQ=255)
                            cgcgAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGGTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                         <  3:63646/149‑1 (MQ=255)
                                 gCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTAACCCGCGAAGAGGCCAGCACAg                                                                                                                    >  3:250334/1‑149 (MQ=255)
                                  cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                   >  7:237142/1‑149 (MQ=255)
                                     gCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGa                                                                                                               <  2:193440/149‑1 (MQ=255)
                                           gggCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAAc                                                                                                          <  8:237142/148‑1 (MQ=255)
                                                    tCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGccc                                                                                                <  7:20718/149‑1 (MQ=255)
                                                       cTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTCCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATa                                                                                                             <  1:124341/133‑1 (MQ=255)
                                                        tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGGTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCAGCCCGTAGGAGCGGGTTCACCCGCCAAGAGGCCAGCACCGGCGATCAACATCAATGCCCAGGa                                                                                            >  2:83755/1‑149 (MQ=39)
                                                            gataCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑cgc                                                                                       <  8:92369/145‑1 (MQ=255)
                                                               aCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCGCAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCg                                                                                   <  3:203624/149‑1 (MQ=255)
                                                                  cTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGCATAAACATCAATGCCCAGGA‑‑CGCGCCGGt                                                                                 <  1:22781/149‑1 (MQ=255)
                                                                          gggggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                         <  3:9254/148‑1 (MQ=255)
                                                                                            cgatctgcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCCAGCCGTCa                                                     <  5:238136/143‑1 (MQ=255)
                                                                                                   cgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAAAGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCa                                                     >  6:238136/1‑142 (MQ=255)
                                                                                                    gATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                             <  7:21660/149‑1 (MQ=255)
                                                                                                           tttttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTGACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGC‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAg                                      >  5:238140/1‑149 (MQ=255)
                                                                                                           ttatcgccTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGACGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                       <  5:122522/141‑1 (MQ=18)
                                                                                                           ttatcgccTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGACGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCa                                         >  6:122522/8‑146 (MQ=18)
                                                                                                              ttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAgcg                                    <  3:268080/148‑1 (MQ=255)
                                                                                                                   tgtgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGc                               >  1:247181/1‑148 (MQ=255)
                                                                                                                     tgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACCCGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAg                                 <  1:83755/146‑1 (MQ=255)
                                                                                                                        tGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCgg                         >  3:110106/1‑149 (MQ=255)
                                                                                                                               ttCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCg                  >  5:64110/1‑149 (MQ=37)
                                                                                                                                    gggTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa             >  8:267097/1‑149 (MQ=255)
                                                                                                                                     ggTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTc                                          >  8:10259/1‑119 (MQ=25)
                                                                                                                                     ggTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTc                                          <  7:10259/119‑1 (MQ=25)
                                                                                                                                      gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc            >  3:256212/1‑148 (MQ=255)
                                                                                                                                       tgaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa          >  3:258864/3‑149 (MQ=255)
                                                                                                                                       tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                         >  7:95748/1‑87 (MQ=255)
                                                                                                                                       tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                         <  8:95748/87‑1 (MQ=255)
                                                                                                                                       tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                         <  4:253706/87‑1 (MQ=255)
                                                                                                                                       tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                         >  3:253706/1‑87 (MQ=255)
                                                                                                                                        gaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGAt         >  3:219976/2‑149 (MQ=38)
                                                                                                                                        gaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc              <  7:267585/143‑1 (MQ=255)
                                                                                                                                        gaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc              >  8:267585/2‑144 (MQ=255)
                                                                                                                                          aCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCg       >  6:99531/1‑149 (MQ=255)
                                                                                                                                           cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc        <  4:169292/147‑1 (MQ=255)
                                                                                                                                           cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc        >  3:169292/1‑147 (MQ=255)
                                                                                                                                            ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCCACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                    <  1:196979/135‑1 (MQ=255)
                                                                                                                                            ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                    >  2:196979/1‑134 (MQ=255)
                                                                                                                                            ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc     >  5:267430/1‑149 (MQ=255)
                                                                                                                                               cccgTAGGAGCGGGTTCACCCGCTAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGc  <  7:244505/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGC  >  NC_002947/4796468‑4796759

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: