Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,309 A→G intergenic (+299/‑180) sokC → / → nhaA small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,3090AG93.3% 114.6 / 0.4 46intergenic (+299/‑180)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (23/19);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

GCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATC  >  NC_000913/17068‑17545
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
gCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAggg                                                                                                                                                                                                                                         >  1:60793/1‑251 (MQ=255)
       ttAACCCATTCTGTTTTATCGATTCTAATCCTTAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATTATTTAATTATTTAACTTGAATTTATCTCATTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑AATATCTGGGCTTAGGGTCAc                                                                                                                                                                                                                                     <  2:396938/248‑1 (MQ=255)
         aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTg                                                                                                                                                                                                                                >  1:417908/1‑251 (MQ=255)
          aCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGa                                                                                                                                                                                                                               <  1:79782/251‑1 (MQ=255)
              attctGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAgcgc                                                                                                                                                                                                                           >  1:37814/5‑251 (MQ=255)
                     tttATCGATTCTACTCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGc                                                                                                                                                                                                                    >  2:330748/1‑251 (MQ=255)
                     tttATCGATTCTAATCCTGACGCCGCCTCGCATTTTTGTGGCGTAATTTTTTACTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTATTTAATCCCTCAGAAATCCTAACCTTAAGCTATG‑ATTATCTTGCCTTAGGGTCACTCGTGAGCGCTTAccgc                                                                                                                                                                                                                    >  2:91218/1‑248 (MQ=255)
                       tATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTGTCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCg                                                                                                                                                                                                                  >  2:406917/1‑251 (MQ=255)
                                        aaGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCaa‑aa‑aactcatctc                                                                                                                                                                                                >  1:505923/1‑241 (MQ=255)
                                             gCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAAGCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCt                                                                                                                                                                                         <  2:417908/251‑1 (MQ=255)
                                                                    tttttAATTATTTAATTATTTAACTTTAATTTAACTCTTTAACGTAATTATTGACGACAATCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGATTGACAATAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCTCCTTAATCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGTTTACAGCCGTCAA‑AA‑AAC‑TCATCCCACCGTTGATGTCGCAAATTCTTCAAtgg                                                                                                                                                                   <  2:494893/250‑3 (MQ=255)
                                                                          aTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTATAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGATGATGGCGCACATTATTCAATGGATCGTa                                                                                                                                                             >  2:397388/1‑250 (MQ=255)
                                                                                     attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTTGATTATTTTTGAAATATTGTCCAAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCGTAAGCTATG‑ATTATCGGGGCTTAGGGTCACTCGTGAGCGCTGTCAGCCGTCAA‑AA‑AAC‑ACATCTCTCCGCGGATGGCGGAAATTCTTCAATGGATCGTAAAAAACGagtt                                                                                                                                                  <  2:308009/250‑4 (MQ=255)
                                                                                     attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAACCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGagtt                                                                                                                                                  <  2:151532/250‑4 (MQ=255)
                                                                                          aCTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTACTTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AACA‑C‑TAATATCACCGCTGATGGCACAAATTCTTCAATGGATCGTAAAAAACgtgtaattcct                                                                                                                                            >  2:278499/1‑241 (MQ=255)
                                                                                                  ttATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGTAATCTTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGACGCAAATTCTTCAATGGCTCGTAAAAAACGAGTTATTCATACACtat                                                                                                                                     >  2:426814/1‑250 (MQ=255)
                                                                                                      ctctTTATCGCACATATTGACGACAAGCTGGATTATTTTTGAAATATTGTCATAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCCCCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTAGTGAGCGCTTACAGCCGTCAA‑AA‑A‑CCTCATATCACCGCTGATGGCGCAAATTCTTCAATGGTTCGAAAACAAAGAGTTATTCCTACACTATCATCt                                                                                                                                >  2:378540/1‑251 (MQ=255)
                                                                                                            agcGCAAATATTGACGACAAGCTGGATTATTTGTGTAAGATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGTGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtgt                                                                                                                           <  2:204937/248‑3 (MQ=255)
                                                                                                            aTCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtttt                                                                                                                          >  2:323167/1‑251 (MQ=255)
                                                                                                                   ttattGACGACAAGCTGGATTATTGTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGACATCCTCACATTAAGCTATG‑ATTATCTGGACTTAGGGTCACTCGTGAGCGCTTACATCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACTAGTTATTCCTACACTATAATCTGATTTTAACgatg                                                                                                                   >  2:378945/1‑251 (MQ=255)
                                                                                                                             cAAGCTGGATTATTTTTGAAATATTGCCCTACCAATTACCGCCGCCTGCCAAGAAATTAAT‑ATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCggg                                                                                                        <  1:223063/251‑1 (MQ=255)
                                                                                                                             aaaGCTGGATTATTTTGTAAAAATTTTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTGCCTAATTAATCCCTCAGAAATCCTCACCTTAAGC‑ATGGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGTCGTCAA‑AA‑AAC‑GCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCgg                                                                                                         <  2:303611/250‑1 (MQ=255)
                                                                                                                                  tGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGag                                                                                                                                                    >  1:56401/1‑202 (MQ=255)
                                                                                                                                  tGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAAGCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGag                                                                                                                                                    <  2:56401/203‑2 (MQ=255)
                                                                                                                                     attattTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAAt                                                                                                 <  1:341175/251‑1 (MQ=255)
                                                                                                                                      ttattTTTGAAATCTTGTCCTAACAAATAACGCCGAATGAAAAGAACTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGATTAGGGTCACTCGTGAGCGCTTACAGCCGTCAACAA‑A‑C‑TCATCTCACCGCTGATGGCGCAACTTCTTCAATGGATCGTAAAAAACGTGTTATTCCTACACTATAACATGATTTTAACgatgctttattcggggaacaat                                                                                                 >  1:365850/1‑232 (MQ=255)
                                                                                                                                           tttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGAAAAATAGTaa                                                                                            >  1:293822/1‑250 (MQ=255)
                                                                                                                                             tGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTaaaaa                                                                                         <  1:190485/250‑1 (MQ=255)
                                                                                                                                                   atTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGTCAA‑AATA‑C‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAGACTAGTAAACACGATCt                                                                                   >  1:244306/1‑251 (MQ=255)
                                                                                                                                                            acaaatcatgCCGACTGACAAGAAATTACTTATTACCCTTCCTAATTAACCCCTCACAAATACTCACCTTAAGCTATG‑ATTTTCTGGGCTTAGGCTCACTAGTGACCGCTTACAGCCGTCAA‑AA‑A‑CATCATCTCACCGCTGACGGCGCAAATTCTTCAATGGATCGTAAAAAAAGAATTATTCCTACACTATAATCTGATTTTAACAATGATTTATGCGGGGTAAACTAGTAAAAAAGCTCTATCCACCTg                                                                          >  2:236474/10‑251 (MQ=255)
                                                                                                                                                            aaaaataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCAACTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTg                                                                          <  2:427193/243‑1 (MQ=255)
                                                                                                                                                             cAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATACTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCa                                                                                                                                                                                                                                      <  1:53687/97‑1 (MQ=255)
                                                                                                                                                             cAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATACTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCa                                                                                                                                                                                                                                      >  2:53687/1‑97 (MQ=255)
                                                                                                                                                                   acGCCGACTGACAAGAAAGTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAgaga                                                                   <  1:483827/250‑1 (MQ=255)
                                                                                                                                                                               aagaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGa                                                       <  2:37814/248‑1 (MQ=255)
                                                                                                                                                                                 aaaaTTAATTATTACTTTACCTAATTAATCCCTCAGAAAACCTCACCTTAAGCTATT‑ATTATCTGGGCTTAGGGTCACTCGTGATCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCATATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAATAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGaaa                                                     <  2:62093/250‑1 (MQ=255)
                                                                                                                                                                                        attattACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGATATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAACCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTCACGATGATTTGTGCGGGGTAAATTAGTAAAAACGATCTTTTCACCTGAAAGAGAA‑Taaaaacgacaaacatatgc                                             >  2:444916/1‑237 (MQ=255)
                                                                                                                                                                                              aCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATACAAAGTGAAACATCTGCATCGAt                                        >  1:270606/1‑251 (MQ=255)
                                                                                                                                                                                                       aattaatCCCTCAAAAATCCTCACCTTAAGCTATG‑ATTAACTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATGTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCa                               <  2:99928/251‑1 (MQ=255)
                                                                                                                                                                                                         ttaatCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGt                             >  2:262957/1‑251 (MQ=255)
                                                                                                                                                                                                           aatCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTATTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAACTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACCTCTGCATCGATTCTTTAGCAGTgg                           >  2:123799/1‑250 (MQ=255)
                                                                                                                                                                                                                  cagaaaTCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                     <  1:406917/246‑1 (MQ=255)
                                                                                                                                                                                                                      aaTCCTTACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCTTGAGCGCTTACAGCCGTAAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCTCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAg                <  2:401456/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaTCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                  >  1:217932/1‑249 (MQ=255)
                                                                                                                                                                                                                               cctAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCGCGGGAGGCATTATTc       <  2:302391/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc  >  2:367861/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
GCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATC  >  NC_000913/17068‑17545

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