Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,622,2500AG36.0% 77.1 / 13.2 51E161G (GAG→GGG) yjjJtoxin YjjJ
Reads supporting (aligned to +/- strand):  ref base A (4/28);  new base G (18/0);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.05e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATT  >  NC_000913/4622011‑4622497
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                               
gtgCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGAtttat                                                                                                                                                                                                                                                      >  1:444228/1‑248 (MQ=255)
  gCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGTAGGGGAATTATCAg                                                                                                                                                                                                                                                    <  2:670440/251‑1 (MQ=255)
  gCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTGGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCTGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATTATATACGTCTCGGTCAGGAAGAAGATGTGCTTTACTCGCTGACCGTATTTATCGGTGAATATACTGGCTGTTGGTGGGTTGGGGAGGGGAAGTATCa                                                                                                                                                                                                                                                     <  2:348261/250‑1 (MQ=255)
       ccGCAGGGAAGTTGTCTGGTTACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGACGATGTGCTCTACGCCCTGACCGTATTTAACGGAGAATATACTGGCGGTGGGTTGGTTGGGGTGGGGAGTTATAAGGGATg                                                                                                                                                                                                                                               >  2:519566/1‑251 (MQ=255)
         gCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTTGTCGGGGGGGGGGAtttttcaggatgtt                                                                                                                                                                                                                                             >  1:120129/1‑239 (MQ=255)
               aaGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTGTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTTCTCTACGCCCTGACCGTATTTAACGGTGAATATTCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGTA‑TTACtc                                                                                                                                                                                                                                      <  2:42471/251‑2 (MQ=255)
                    gTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGcaca                                                                                                                                                                                                                                  <  1:513581/250‑1 (MQ=255)
                      ccggTAACAGGCGGAGATGGCGTCGAACAGTGGTTTGATGGTTTGCCCTGGTAGTTGACCGATATCCGAACGTAGGGCTTGTTAGGGGGGGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTGTGGCAGGAAGAAGATGTGCGCTGCGGCCGGACCGTAGTTAACGGTGAAGATACTGGCGGTTGGTTGGGCGGGGAGGGGAATTATCAGCGATGGA‑TTTCTGCacaac                                                                                                                                                                                                                                <  2:329528/248‑1 (MQ=255)
                        gggaaCAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGGGAATATACTGGCGGTTGGTGGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCacaac                                                                                                                                                                                                                                <  2:644905/245‑1 (MQ=255)
                        ggTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTGTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACAcc                                                                                                                                                                                                                             <  2:429188/251‑1 (MQ=255)
                            aCAGGCGCAGATGGTGGCGAACAGTGGTTTGGTGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGc                                                                                                                                                                                                                         <  2:653167/251‑1 (MQ=255)
                              aGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGTCGGGGAGGGGGGTTATCAGCGATGGG‑TTACTTCCCAACACCCTGCgg                                                                                                                                                                                                                       >  1:608646/1‑251 (MQ=255)
                                gcgcAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTTTCAGCGATGGA‑TTACTGCacaaacacctgcggg                                                                                                                                                                                                                      >  1:322559/1‑239 (MQ=255)
                                   cAGATGGCGACGAACAGTGGGTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGGCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGAGGGA‑TTACTGCGCAACGCCCTGCGGAAA‑tt                                                                                                                                                                                                                 <  2:131903/251‑1 (MQ=255)
                                    aGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGAGCGTATTTAACGGTGAATATACTGGCGGGTGGGTGGTCGGGGAGGGGAATGATCAGCGATGGA‑TTACTGCacaaca                                                                                                                                                                                                                               <  2:414459/237‑1 (MQ=255)
                                    aGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGGTTATCAGCGATGGA‑TTACTGCacaacaagatcggaaaagg                                                                                                                                                                                                                  >  1:414459/1‑237 (MQ=255)
                                     gATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGGCGGGGGGGGGGATTTTTTGCGAGGGA‑TTACTTCacaaacccccgcggagattac                                                                                                                                                                                                                >  2:101780/1‑234 (MQ=255)
                                        ggCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGGCCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCctct                                                                                                                                                                                                            <  2:395081/251‑1 (MQ=255)
                                        ggCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTTGGGGGGGGGAAATATTAACGATGGA‑TTACTGCACAAAACCCTGCGGAAA‑TTTctct                                                                                                                                                                                                            >  1:513789/1‑251 (MQ=255)
                                                              gTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGAATATCAGCGAAGGG‑GTATTGCACAACCACCTTCGGGAA‑TACTCCTGGGTCAAAAAATCAccccatac                                                                                                                                                                                      >  2:49756/1‑246 (MQ=255)
                                                                         tATTTGACCGATCTCCGACCGCAGGGCTTTTTAGAGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGTGCGGGGGGGGGTATTTTCTACGATGGAGTTA‑TGCCCCACACCCTGCCGAAA‑ATACTCTGGATCAAAAAATCAcccccttccacacaccggg                                                                                                                                                                           >  2:542166/1‑235 (MQ=255)
                                                                                        cGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGc                                                                                                                                                            <  1:4608/251‑1 (MQ=255)
                                                                                              cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGCCCGTATTTAACGGTGAATATACTGGCGGCTGGTTGGGCGGGGGGGGGAATTATCCGCGATGGG‑TTTCTGCCCAAAACCCTGCGGGAA‑ATCCTCTTGATCAAAAACTCACCCACTACCAAAAGCCTGCAAGTGGTGCAATGGCAgggg                                                                                                                                                       >  1:289928/1‑248 (MQ=255)
                                                                                              cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGA‑TTACTGCACAAACCCCTGCGGGAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAggggg                                                                                                                                                      >  1:322113/1‑248 (MQ=255)
                                                                                                      tttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGCTGGG‑TTTCTGCCCAAACCCCCGCGGAAA‑TTTCTCTTGATCCAAAAATCACCCATTACGAACCGCTGGCAAGGGATGCAATGGCAGG‑A‑G‑GAATTTTgg                                                                                                                                           >  1:36607/1‑251 (MQ=255)
                                                                                                          gggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTAGTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCGGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCAGTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGttc                                                                                                                                       <  1:75377/251‑1 (MQ=255)
                                                                                                           ggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGCTGGTCGGGGGGGGGGATTATCAGAGATGGA‑TTAATGCCCCACACCCTGCGGAAA‑TTCCTCTGGGTCAAAAACTCACCCCTTACGAACAGCTGGCAAGTGATTCACTGGCAGG‑T‑G‑AAATTGTGGGttct                                                                                                                                      >  1:82024/1‑251 (MQ=255)
                                                                                                              gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTtcttct                                                                                                                                   <  1:613376/251‑1 (MQ=255)
                                                                                                              gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTtcttct                                                                                                                                   <  1:689867/251‑1 (MQ=255)
                                                                                                                                     gcgcTACTGAAGTTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTTCGCCATGTCCGTATTTAACGGTGAATATACTGGCGGTTGGTTGATCGGGGAGGGGAATTATCAGCGATGTA‑TTACTGCACAACACCCTGCTGAGA‑TTCCTCTGGATCAAAAACTCACCCATTGCGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑T‑ACATTTTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCaaaacc                                                                                                           <  2:20977/251‑3 (MQ=255)
                                                                                                                                          acggAATCTGACTGTGGATATCCGT‑TTTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGTTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑GTTACcta                                                                                                    <  2:459075/248‑1 (MQ=255)
                                                                                                                                            tGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATCAGCGATGGA‑TTAATGCACAACACCCTGCGGGAA‑TTCCTCTGGAACAAAAAATCCCCCATTACGAACAGCTGGCAAGTGATGCACCTGCCGG‑A‑G‑AAATTGTGGGGTCTTCTGCGGGGGGGGGGC‑AGCCAAAA‑ATTACctac                                                                                                   >  1:477323/1‑251 (MQ=255)
                                                                                                                                              aaTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGGTGGG‑TTACTGCACAACACCCTGCGGGAA‑TTTCTCTGGATCAAAAACTCACCCACTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑GGAAATT‑TGGGTTCTTCTTCGGGCGGCGGGC‑AGCCAAAA‑TTTACctaat                                                                                                  >  1:30738/1‑248 (MQ=255)
                                                                                                                                               aTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTTGTTGGTCGGGGGGGGGGATTATCAGGGGTGGG‑TTTCTGCACAACACCCTGCGGAAAATTCC‑CTGGATCAAAAAAACCCCCACTACACAAAGATGGCAAGGGAGACACTGGCAGGGA‑G‑AAAT‑GTGGGGTCTTTTGCGGGGGGGGAGG‑AGACCAAA‑TTTAcaatctat                                                                                                >  1:523326/1‑244 (MQ=255)
                                                                                                                                                   gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATg                                                                                               <  2:30738/248‑1 (MQ=255)
                                                                                                                                                    actgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTCTTTAACGGTGAATATACTGGCGGATGGTTGGTCGGGGGGGGGAATTAACAGCGATCGA‑TTACTGCACAACACCCTGCGGGAA‑TTCCTCTGGATCAAAAACCCACCCATTACGAACAGCTGGCAAGTGATGCACCGGCAGG‑A‑G‑GAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATTcac                                                                                            >  1:108616/1‑250 (MQ=255)
                                                                                                                                                                cGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTAGGCACAAACGCCGTCAgg                                                                               <  1:519566/251‑1 (MQ=255)
                                                                                                                                                                     ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑GTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGTCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCa                                                                             <  1:590129/248‑1 (MQ=255)
                                                                                                                                                                       ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGGTCTGCTGCTGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAATCGCCGTCAGGCGATaaa                                                                        <  2:305345/251‑1 (MQ=255)
                                                                                                                                                                               gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGCTGGGCGGGGGGGGGAATTATCAGCGATGGT‑TTACTGCCAAACACCCCGCGGGAA‑TTCCTCTGGATCAAAAAATCAACCATTACGAACAGCTGGCAAGGGATGCACTCGCAGG‑AAG‑AAATT‑TGGGTTCTTCTGCGGGCGGCGAGC‑AGACCAAA‑TTTACCTACTATGGACAAACGACGCCACGCAAAAAACAtaggtt                                                                >  2:578556/1‑246 (MQ=255)
                                                                                                                                                                                            aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCa                                                    <  2:7234/249‑1 (MQ=255)
                                                                                                                                                                                             cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGGTGGTCGGGGAGGGGTAATATCCGCGATGGA‑TTACTGGACAACACCCTGCGGGAA‑ATCCTCTGGGTTCAAAACTCCCCCCTTACGAACAGCCGGCAAGAGATGCACTGGCAGG‑A‑G‑AAAATGTGGGTTTTTCTTCGGGGGGGGGGG‑CGCCAAAAGTTTAC‑TACTATTCCCCAACGACGTCAGGCAATAAACCTGTTTTGGGG‑AAAttacccc                                                 >  2:525681/1‑246 (MQ=255)
                                                                                                                                                                                             cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTTTTAGCGGTGGA‑TTAATGCACAACACCCCGCGGAAA‑TTTCTCTGGGTCAAAAAATCACCCATTACGAACAGCCGGCCAGAGATTCACTGGGAGG‑G‑G‑AAAATTTGGGGTTTTCTTCGGGCGGCGAGCCAACAAAAA‑TTT‑CCCACTATGAAAAAACGCCGGAAGGCAATAAACATGTGTTTGTGGAAATTCAcc                                                  >  1:112726/1‑251 (MQ=255)
                                                                                                                                                                                                   gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACcac                                           <  1:467288/251‑1 (MQ=255)
                                                                                                                                                                                                           ttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGAATATCAACGATTGC‑TTACTGCCCAACAACCCGCGGAAA‑ATCCTCTGGAACCAAAACCCACCCAATACCGACAGCTGGCAAGGGATGCACCGGAAGG‑A‑GGAAATT‑TGGGGTTTTCTGCGGGCGGCGAGG‑CGGCAAAA‑ATTACCTACCTTTCACAAAAGCCATCAAGCAATAAACAATTGTTGGTG‑GAATTCACCggacacccccaaacc                                   >  1:185134/1‑237 (MQ=255)
                                                                                                                                                                                                                         gcTGGCGGTGGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCGGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAAGTCACCGTACCACAGCAAACCGCGGTCAGCCAACg                     <  2:120129/250‑1 (MQ=255)
                                                                                                                                                                                                                          cTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑C‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCaa                                                                            >  1:367286/1‑196 (MQ=255)
                                                                                                                                                                                                                          cTGGCGGTTGGTTGGGCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑GTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGAGGCACTGGCAGT‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTAGGCACAAACGCCGTCAGGCaa                                                                            <  2:367286/196‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGt                    <  1:392595/249‑1 (MQ=255)
                                                                                                                                                                                                                                gttggttgGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGt              <  1:49756/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggtCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TGCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCGTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGc       <  1:364119/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gggaggggaATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAAtt  <  1:251009/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                               
GTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGA‑TTACTGCACAACACCCTGCGGAAA‑TTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGG‑A‑G‑AAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAA‑TTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTG‑AAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATT  >  NC_000913/4622011‑4622497

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