Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,257,593 A→G 6.1% D9D (GAT→GAC) ‡ aruR ← transcriptional regulator AruR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,257,5930AG 6.1% 196.2 / 2.3 65D9D (GAT→GAC) ‡aruRtranscriptional regulator AruR
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CCGGGCAGGCGGATGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCACCGG  >  NC_002947/4257447‑4257733
                                                                                                                                                   |                                                                                                                                             
ccGGGCAGGCGGATGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATc                                                                                                                                             >  2:118562/1‑149 (MQ=255)
      aGGCGGATGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGATAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGAt                                                                                                                                       >  6:104600/1‑149 (MQ=255)
        gCGGATGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCGAGTGACGGCGCAGTGTACGACGTAGCCTGCCCGGGCCAGGTAGGCATGCAGCAAGTCACGAACAAGCGGAACGTCATCGACGAt                                                                                                                                       >  1:110940/1‑148 (MQ=255)
          ggATGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCGATGACGATccgc                                                                                                                                   <  6:289430/149‑4 (MQ=255)
            aTGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCGATGACGATccgctt                                                                                                                                 <  6:302230/149‑6 (MQ=255)
                ccAGCAGCACGAGGTCGACGTCTTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATAccc                                                                                                                               >  5:365951/1‑147 (MQ=255)
                   gcagcaCGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATcgac                                                                                                                                          <  5:324599/133‑1 (MQ=255)
                   gcagcaCGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATcgac                                                                                                                                          >  6:324599/1‑133 (MQ=255)
                        aCGAGGGCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑t                                                                                                                    <  7:276239/149‑1 (MQ=255)
                             gTCGACGTCCTGTCTTACAAGCATGGCTTCGGCCCTTTCTGCT‑GTTTCGGAGAAGTGTACTTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCAGGTACTAGCGGTTCGTCATCGTCGATCAATTCCCGTTGCGTGT‑TTGCgt                                                                                                               <  8:186716/149‑1 (MQ=255)
                             gTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCgt                                                                                                               <  1:118562/149‑1 (MQ=255)
                             gTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCgt                                                                                                               <  5:366962/149‑1 (MQ=255)
                              tCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCgtg                                                                                                              >  4:305847/1‑149 (MQ=255)
                               cGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCgtgt                                                                                                             >  3:232286/1‑149 (MQ=255)
                                  cGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGtgttgt                                                                                                          <  3:305847/149‑1 (MQ=255)
                                      cTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTctct                                                                                                      >  6:50280/1‑149 (MQ=255)
                                       tGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCtt                                                                                                     >  3:309377/1‑149 (MQ=255)
                                         gCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCttt                                                                                                    <  4:81970/148‑1 (MQ=255)
                                           tGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTa                                                                                                 >  5:109814/1‑149 (MQ=255)
                                             ccAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTa                                                                                                 <  4:301072/147‑1 (MQ=255)
                                              cAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTg                                                                                              <  4:173222/149‑1 (MQ=255)
                                               aaGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTgg                                                                                             >  2:200078/1‑149 (MQ=255)
                                               aaGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTgg                                                                                             >  5:266026/1‑149 (MQ=255)
                                                aGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGgt                                                                                            <  7:42678/149‑1 (MQ=255)
                                                 gCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGgtg                                                                                           <  6:365951/149‑1 (MQ=255)
                                                   aGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGtgtt                                                                                         <  6:266026/149‑1 (MQ=255)
                                                     ggCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGtgttgt                                                                                       <  6:109814/149‑1 (MQ=255)
                                                            gCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCgt                                                                                                               <  6:282267/118‑1 (MQ=255)
                                                            gCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCgt                                                                                                               >  5:282267/1‑118 (MQ=255)
                                                            gCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTaca                                                                                 >  3:64854/1‑148 (MQ=255)
                                                            gCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTaca                                                                                 >  5:274334/1‑148 (MQ=255)
                                                             cccTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTacaca                                                                               >  1:148637/1‑149 (MQ=255)
                                                               cTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACt                                                                             <  8:71280/149‑1 (MQ=255)
                                                               cTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACt                                                                             >  3:108983/1‑149 (MQ=255)
                                                                     gcg‑gTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGa                                                                        >  6:298929/1‑148 (MQ=255)
                                                                            tCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACagcagc                                                                 <  2:110940/149‑1 (MQ=255)
                                                                                       tACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTAc                                                      <  5:179224/149‑1 (MQ=255)
                                                                                          gtcgtAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTg                                                    >  2:66839/1‑148 (MQ=255)
                                                                                                gCCTTCCCGGGACAGGTAGGCCTGCAGCTAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGc                                             <  5:100479/149‑1 (MQ=255)
                                                                                                gCCTTCCCGGGACAGATAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGc                                             <  5:104600/149‑1 (MQ=255)
                                                                                                   ttCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGAt                                          >  4:94073/1‑149 (MQ=255)
                                                                                                    tCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGAt                                          >  5:252048/1‑148 (MQ=255)
                                                                                                         ggACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAg                                    <  4:232286/149‑1 (MQ=255)
                                                                                                              ggTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATcgac                                                                                                                                          >  8:133796/1‑43 (MQ=255)
                                                                                                              ggTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATcgac                                                                                                                                          <  7:133796/43‑1 (MQ=255)
                                                                                                               gTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCg                              <  4:309377/149‑1 (MQ=255)
                                                                                                               gTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGGGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCg                              <  4:108983/149‑1 (MQ=255)
                                                                                                                     cTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAgg                        >  5:132740/1‑149 (MQ=255)
                                                                                                                        cagcaAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGc                       <  2:46606/147‑1 (MQ=255)
                                                                                                                            aaGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAgc                 <  5:298929/149‑1 (MQ=255)
                                                                                                                             aGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAgcg                >  4:198498/1‑149 (MQ=255)
                                                                                                                                   gAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGt                                                         >  1:20120/1‑102 (MQ=255)
                                                                                                                                   gAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGt                                                         <  2:20120/102‑1 (MQ=255)
                                                                                                                                   gAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCg          >  8:152665/1‑149 (MQ=255)
                                                                                                                                    aaCAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCgg         >  1:172865/1‑149 (MQ=255)
                                                                                                                                       gtattgATCGTCGATGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCa      >  7:34308/6‑149 (MQ=255)
                                                                                                                                       aaGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCa      <  6:274334/149‑1 (MQ=255)
                                                                                                                                       aaGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCa      <  4:64854/149‑1 (MQ=255)
                                                                                                                                         gCGGATCGTCATCGACGATCAATACCCGAGGCGTTATTTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACAAACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATAAACAGCCCTCGTACAGGCTGAAGCGCCTGCGGCCAc     >  1:67122/1‑149 (MQ=255)
                                                                                                                                          ttgATCGTCGATGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCACCg   >  5:227325/3‑149 (MQ=255)
                                                                                                                                           ggATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCACCgg  <  1:200078/149‑1 (MQ=255)
                                                                                                                                              tCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCt                                <  3:368638/116‑1 (MQ=255)
                                                                                                                                              tCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCt                                >  4:368638/1‑116 (MQ=255)
                                                                                                                                                gTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAgcgc               >  6:111468/1‑131 (MQ=255)
                                                                                                                                                gTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAgcgc               <  5:111468/131‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CCGGGCAGGCGGATGTCCAGCAGCACGAGGTCGACGTCCTGGCTGCCAAGCAGGGCTTCGGCCCTTTCGGCG‑GTGTCGGCGCAGTGTACGTCGTAGCCTTCCCGGGACAGGTAGGCCTGCAGCAAGTCACGAACAAGCGGATCGTCATCGACGATCAATACCCGAGGCGTCA‑TTGCGTGTTGTCTCTTTTTATTGGTGTTGTTCTACACACTCGCGACAGCAGCGTACTCGTTACTGAACAGGCGATCAACAGCCCTCGGACAGGCTGAAGCGCCTGCGGCCACCGG  >  NC_002947/4257447‑4257733

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