Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,244,996 G→A 8.2% V135M (GTG→ATG)  hmgR → Hmg transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,244,9960GA 8.2% 172.3 / 5.1 61V135M (GTG→ATG) hmgRHmg transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (32/24);  new base A (3/2);  total (35/26)
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.74e-01

GCTGGGCAGCATCGATGTGCTGAAGGTGGCGGCGATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGCCAGTGACTATCACA  >  NC_002947/5244851‑5245140
                                                                                                                                                 |                                                                                                                                                
gcTGGGCAGCATCGATGTGCTGAAGGTGGCGGCGATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCgtgg                                                                                                                                               >  8:208713/1‑149 (MQ=255)
   gggCAGCATCGATGTGCTGAAGGTGGCGGCGATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGa                                                                                                                                            >  1:237832/1‑149 (MQ=255)
               tgtgCTGAAGGTGGCGGCGATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGAtcggtt                                                                                                                                <  8:201830/149‑1 (MQ=255)
                   cTGAAGCTGGCGGCGCTGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGACAGCTGCTTCATCGCCGGGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGAtcggtt                                                                                                                                <  6:233756/145‑1 (MQ=255)
                            gcggcgATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGAtcggtt                                                                                                                                <  8:56130/136‑1 (MQ=255)
                            gcggcgATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGAtcggtt                                                                                                                                >  7:56130/1‑136 (MQ=255)
                            gcggcgATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGtt                                                                                                                   >  6:309009/1‑149 (MQ=255)
                             cggcgATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTg                                                                                                                  <  3:187303/149‑1 (MQ=255)
                             cggcgATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTg                                                                                                                  >  3:191960/1‑149 (MQ=255)
                                   tGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGt                                                                                                            >  7:137352/1‑149 (MQ=255)
                                    gCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGtt                                                                                                           >  7:40513/1‑149 (MQ=255)
                                           tCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCAc                                                                                                    <  4:191960/149‑1 (MQ=255)
                                             gCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCg                                                                                                  <  7:11791/149‑1 (MQ=255)
                                              cAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGa                                                                                                                                      >  4:85805/1‑112 (MQ=255)
                                              cAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGa                                                                                                                                      <  3:85805/112‑1 (MQ=255)
                                                   gCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTgg                                                                                            >  4:201474/1‑149 (MQ=255)
                                                    cGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGt                                                                                           <  5:309009/149‑1 (MQ=255)
                                                        aTGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTc                                                                                       >  1:137748/1‑149 (MQ=255)
                                                         tGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCATGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTcg                                                                                      >  5:215980/1‑149 (MQ=255)
                                                         tGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCATGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTcg                                                                                      >  5:125522/1‑149 (MQ=255)
                                                              tCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCAcc                                                                                                   <  8:228050/131‑1 (MQ=255)
                                                              tCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCAcc                                                                                                   >  7:228050/1‑131 (MQ=255)
                                                                aaCGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCg                                                                                                  >  6:152934/1‑130 (MQ=255)
                                                                aaCGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCg                                                                                                  <  5:152934/130‑1 (MQ=255)
                                                                aaCGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTa                                                                               <  7:208713/149‑1 (MQ=255)
                                                                aaCGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTa                                                                               >  7:157603/1‑149 (MQ=255)
                                                                          gctTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCcgag                                                                     >  5:324287/1‑149 (MQ=255)
                                                                          gctTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCcgag                                                                     >  4:86377/1‑149 (MQ=255)
                                                                           ctTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCcgagc                                                                    <  3:201474/149‑1 (MQ=255)
                                                                             tCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTAtt                                                                                        >  3:218852/1‑127 (MQ=255)
                                                                             tCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTAtt                                                                                        <  4:218852/127‑1 (MQ=255)
                                                                                      tgtgGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAg                                                         <  8:40513/149‑1 (MQ=255)
                                                                                        tgGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAgg                                                        <  1:133299/148‑1 (MQ=255)
                                                                                         ggggCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAg                                                      >  4:262476/1‑149 (MQ=255)
                                                                                          gggCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAg                                                      <  8:157603/148‑1 (MQ=255)
                                                                                          gggCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGc                                                     >  7:132573/1‑149 (MQ=255)
                                                                                           ggCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGAGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTAGTCGCCGGCTCCCTTCCCGAGCGAGAGACCCAGGAGCt                                                    >  4:274358/1‑149 (MQ=255)
                                                                                            gCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTg                                                   >  1:290268/1‑149 (MQ=255)
                                                                                                                  aaGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTc                                                                                                       <  5:136685/75‑1 (MQ=255)
                                                                                                                  aaGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTc                                                                                                       >  6:136685/1‑75 (MQ=255)
                                                                                                                  aaGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCg                                <  4:38821/146‑1 (MQ=255)
                                                                                                                  aaGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCg                                >  3:38821/1‑146 (MQ=255)
                                                                                                                  aaGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCtt                             >  5:43657/1‑149 (MQ=255)
                                                                                                                  aaGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCtt                             >  1:16186/1‑149 (MQ=255)
                                                                                                                   aGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCtt                             <  8:132573/148‑1 (MQ=255)
                                                                                                                     cATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGc                          <  2:237832/149‑1 (MQ=255)
                                                                                                                      aTCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGc                          >  2:46805/1‑148 (MQ=255)
                                                                                                                       tCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGcaa                        >  2:253285/1‑149 (MQ=255)
                                                                                                                       tCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGcaa                        <  6:324287/149‑1 (MQ=255)
                                                                                                                        cGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGcaa                        >  7:312501/1‑148 (MQ=255)
                                                                                                                         gAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGcaac                       >  5:79303/1‑148 (MQ=255)
                                                                                                                          aaCCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGa                                                       <  4:120786/115‑1 (MQ=255)
                                                                                                                          aaCCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGa                                                       >  3:120786/1‑115 (MQ=255)
                                                                                                                          aaCCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGcaaca                      >  2:32236/1‑148 (MQ=255)
                                                                                                                           aCCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACaaa                    >  4:19066/1‑149 (MQ=255)
                                                                                                                            ccTGCCGTGCGCGCGGTCACCATGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAAc                   <  6:215980/149‑1 (MQ=255)
                                                                                                                            ccTGCCGTGCGCGCGGTCACCATGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAAc                   <  6:125522/149‑1 (MQ=255)
                                                                                                                                ccGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGc                <  2:16186/148‑1 (MQ=255)
                                                                                                                                   tGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGCCAGt            <  2:290268/149‑1 (MQ=255)
                                                                                                                                     atctgggTCACCACGGTGACCCAGATCGGTTCGGTGCTGCCGTTGGTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGCCAGTGa          >  1:242430/6‑149 (MQ=255)
                                                                                                                                           gTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGTGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGCCAGTGACTATca    >  7:264665/1‑149 (MQ=255)
                                                                                                                                              aCCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCAACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGCCAGTGACTATcaca  >  1:253966/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCTGGGCAGCATCGATGTGCTGAAGGTGGCGGCGATGCCGTTGTCGCAGTTGCGAGATGAGCTCAACGAAAGCTGCTTCATCGCCGTGTGGGGCAACCAGGGGGCCACGGTGGTAAGCATCGAACCTGCCGTGCGCGCGGTCACCGTGGTGACCCAGATCGGTTCGGTGCTGCCGTTGCTCAGTTCGTCCACCGGGCTGGTATTCGCCGCCTACCTGCCCGAGCGAGAAACCCAGGAGCTGCGCGACAGGGAACTGGCCGCTTTGCAACAAACCGCCAGTGACTATCACA  >  NC_002947/5244851‑5245140

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: ‑