Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,281 G→A intergenic (+271/‑208) 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,2810GA91.3% 105.5 / 3.7 46intergenic (+271/‑208)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (23/19);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01

CCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGA  >  NC_000913/17042‑17520
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ccccaaccccccGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCt                                                                                                                                                                                                                                           <  1:100885/242‑1 (MQ=255)
 cccAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCtt                                                                                                                                                                                                                                          <  1:50099/249‑1 (MQ=255)
    aaGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTACTTTTTTAATGATTTCATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGCTTATTTTTGAACTATTGGCCTAACAAGCATCGCCGCCTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCtaccctta                                                                                                                                                                                                                                         >  2:301432/1‑242 (MQ=255)
     aGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGc                                                                                                                                                                                                                                      >  1:427373/1‑251 (MQ=255)
              cgcTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTAt                                                                                                                                                                                                                            <  1:203224/251‑1 (MQ=255)
                          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)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
CCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGA  >  NC_000913/17042‑17520

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