Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,146 G→A A138A (GCG→GCA)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,1460GA94.9% 156.9 / ‑0.2 61A138A (GCG→GCA) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (27/29);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

TGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTG  >  NC_000913/4630899‑4631395
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
ggCCTGGACTGAAAGATTAGCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCTCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACGTGCGCGGCGTGAATACTGGCGTGGCTTGTTAGTCTTGAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCaaaa                                                                                                                                                                                                                                                        <  1:62715/250‑1 (MQ=255)
  ccTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑ACCCGGAACTACCGGAGCGCCAAT‑TAATTACTACTATGCAAAATg                                                                                                                                                                                                                                                      >  1:63547/1‑251 (MQ=255)
   cTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGGCTTTAGCCCGGAAAA‑GGCCAGCACTACCGGAGCGC‑ACTGTAATTACTACTATGCCAAATgg                                                                                                                                                                                                                                                     >  2:40670/1‑251 (MQ=255)
      gACTGAAAGATTATCTACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTACCCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGTCTGGCGTGGCTTGTTAGCCTTTAGCTCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAac                                                                                                                                                                                                                                                  <  2:455632/251‑1 (MQ=255)
      gACTGAAAGATTATCCCCTTTATCCCTACCTGGAATTCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAAGGGAac                                                                                                                                                                                                                                                  <  2:455641/251‑1 (MQ=255)
      gACTGAAAGATTATCCACTTTATCCCTACCTTGAATACCGTCAGATCCCCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACCTTGTTCGCGCTAACCCCACGCTTCCGCTCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAac                                                                                                                                                                                                                                                  <  1:364007/251‑1 (MQ=255)
      gACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTTGCCCGGAAAAACCCCGGAA‑TACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAac                                                                                                                                                                                                                                                  >  2:177446/1‑251 (MQ=255)
         ggaaagATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGTCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACAcc                                                                                                                                                                                                                                               <  2:471795/246‑1 (MQ=255)
               aTTATCCACTTTATCCCTACCTGGCATACAGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGAGAAGACTGGCGTGGCTTGTTAGACTTTAGACCGGAAAA‑GCCCGGGACTACCGAAGCGG‑AAAGTAATTAATACTATGCAAAATGGAAAACCgggg                                                                                                                                                                                                                                           >  1:357877/1‑248 (MQ=255)
                         ttATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTgaagaa                                                                                                                                                                                                                                <  1:262683/250‑1 (MQ=255)
                          tATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGaagaag                                                                                                                                                                                                                               <  1:21120/250‑1 (MQ=255)
                            cccctccctGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACGACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCt                                                                                                                                                                                                                            <  2:452649/250‑1 (MQ=255)
                                   ccggaaacacttcAGATCACCTATGACCTGAGGAACCAACCGGCGGTGCCTGTCACTAACTTTGTTCGCGCTAACCCACCGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATTAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgg                                                                                                                                                                                                                     <  2:119721/239‑1 (MQ=255)
                                                aGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGGAAA‑ACCCGGAACTACCGAAGCGC‑AATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAAAAGCCTGGCCAGGGGCGAAagaggtg                                                                                                                                                                                                        >  1:270689/1‑248 (MQ=255)
                                                          tgaccTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTAACGAAGCGC‑AAAGTAATTACTACTATGCAAAATGGAACACCGGGACGAGGGAAAAAGCCTGGCACGGGGCGAAAGAGCAGTGGCCAACCg                                                                                                                                                                                              >  1:494586/5‑251 (MQ=255)
                                                               tgatgaATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACACTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAAATTGCGCGGCGTGAAGACTGGCGTTGCTTGTTAGCCTTTAACCCGGAAAA‑GCCCGGAACTACCCAAGCGC‑AATGTAATTTCAACTTTGTAAAATGGAACACCGGAGAGAGTGAAAAAGCCCTGCAAGGGGGGaaaaacattgggcaaacctgcca                                                                                                                                                                                          >  1:191645/1‑230 (MQ=255)
                                                                gatgaATCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAA‑GCCCGGAACTAACGAAGCGC‑AATGTAATTACTAATCAGCACAATGGAACACCGGGCAGAGTGAAGACACCTGGCAAGGGGCGAAAGAGCTGTGGGTAACCgggcaa                                                                                                                                                                                         >  1:196000/1‑246 (MQ=255)
                                                                 atgaATCAACCGGCGGTGACTGTCACTAACGTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAagag                                                                                                                                                                                       <  1:32720/251‑1 (MQ=255)
                                                                           cggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAAc                                                                                                                                                                             <  1:179769/251‑1 (MQ=255)
                                                                              cggTGACTGTCACTATCTTTGTTCGCTCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCTTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCATGGGGCGAAAGAGCTGTGGCTAACCGGCAAGATCCCGCCTAACGcc                                                                                                                                                                          <  2:499609/251‑1 (MQ=255)
                                                                                  gactgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCtgtg                                                                                                                                                                      <  1:261614/247‑1 (MQ=255)
                                                                                      gTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGAc                                                                                                                                                                    <  2:396692/249‑1 (MQ=255)
                                                                                               tttGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGGGACAAGTTAttt                                                                                                                                                           >  1:34337/1‑249 (MQ=255)
                                                                                                 tGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAgt                                                                                                                                                        <  1:118450/250‑2 (MQ=255)
                                                                                                  gTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTATCCTTTAGCCCGGCAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAAAGCTGTGGCTAACCGGCAAGAGCCAGCCGAACGCCTGGGACAAGTTATTTAgggt                                                                                                                                                      >  2:123835/1‑248 (MQ=255)
                                                                                                          tAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGgcgtg                                                                                                                                              >  1:58784/1‑251 (MQ=255)
                                                                                                             ccccACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGCACA‑GCCAGGAACTCCCGCAACGC‑AATGTAAATACTACTATGCAAAATGGAAAACCGGGCAGAGTGCAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCCGCCTAACGCCTGTGACAAGTTATTTAGGGTCTGgcgtgcgt                                                                                                                                           >  1:326873/1‑251 (MQ=255)
                                                                                                                  cGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                                      <  1:462864/251‑1 (MQ=255)
                                                                                                                   gCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGACGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGGAAA‑GCCCCGAACTACCGCAGCAC‑CATGTAATTACTACTATGCAAAAAGGGAAACCGGGCAGAGTGAAGAAGCCCGGAATGGGGCGAAAGAGCTGTGGCAAACCAGCAAATGCCCACCTAACGCCTGTGTCAAGTGATTTAGAGACTGgcgtgcggccggg                                                                                                                                      >  2:192812/1‑244 (MQ=255)
                                                                                                                      tcctccCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCCAAAAGGAACACCGGGCCGAGTGAAGAAGACTGGCACGGGGAGAAAGAGCTGTGGCTAACCGGCAAGAGCCCGCCCAACGCCTGTGACAAGTTATTTAGTGTCTGgcgggggtaaggaaaac                                                                                                                                  >  2:457390/1‑237 (MQ=255)
                                                                                                                      tcctccCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTaaa                                                                                                                                   <  1:166143/250‑1 (MQ=255)
                                                                                                                       cctccCGCACGCACGCTGCAATCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACa                                                                                                                                 <  1:229383/251‑1 (MQ=255)
                                                                                                                        ctccCGCTCGCACGCTGCAATATCGTTTCGTCAATGAACTGGCGCGGCGTGAAGTCTGGCGAGGCTTGTTAGCCTTTAGCCCGGCAAA‑GGCCGGGACTACCGAAGCGC‑CATGTAATTTCTACTATGTAAAATGGAACCCCGGGCAGGGTGAAGAAACCTGGCAAGGGGCGCAAGGGCTTGGGCTAACCGGCCAGAGCCCGCCAAAAGCCCTTGAATAGTTATTTAGTGTCTgggcggccgccagtcacca                                                                                                                                 >  2:230921/1‑233 (MQ=255)
                                                                                                                              cTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCg                                                                                                                          >  1:303836/1‑251 (MQ=255)
                                                                                                                                  cACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTCGCCTTTAGCCCGGAAAA‑GCCCGGCACAACCCGAGCGC‑AATGTAATTAATACTCTGCAACATGGCACCCCGGGCAGAGGGAAGAAGCCTGGCAAGGCGCGAAAGAACTGTGGGTACCCCGGAAGAGCCAGCCTCACGACTGTGACACGTTATTTAATGTCTTGGGTGCGTCAGGTCAACAAGTTcccgtgg                                                                                                                      >  1:34970/1‑246 (MQ=255)
                                                                                                                                      cTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCCTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGCAAGGGGCGAAGAAACTGTGGCTAACCGGCAAGAGCCAGCCTACCGACTGGGACAAGTTTTTTAGGGGCTGGCGGGCGGCAGGGAACAAAGATCCAGTGGCGtg                                                                                                                  >  2:346872/1‑250 (MQ=255)
                                                                                                                                         cAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAACCTTTAGCCCGGAAAA‑GCCCGGAACTACCCAAGCGC‑AATGTAATTACTACTATGCAAAATTGGACACCGGGCAGAGGGAAATAGCCTGGCAAGGAGCGAAAGAGATGTTGGTAACAGGCAAGAGCCAGCCCAAAGACTGTGAAAAGTTATTTAGTGTATGGGGTGCGGCAGGTAAACAAGATCAGCTGGCGtcatt                                                                                                               >  1:363077/1‑247 (MQ=255)
                                                                                                                                          aaTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGGGC‑AATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGGCTGGCCAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCCGCCTAACGCCTTGGACAAGTTATTTAGTGTCTGGCTGGCGTCAGGGAAACAAGGTCCGCTGGccgtttt                                                                                                               >  2:359457/1‑244 (MQ=255)
                                                                                                                                               cGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGc                                                                                                          <  1:322175/250‑1 (MQ=255)
                                                                                                                                               cGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGGGCGTCAGGTAAACAAGATCCGCTGGCGTATTTaga                                                                                                            <  2:462772/248‑1 (MQ=255)
                                                                                                                                                       cAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCGGGCCAGGGGCGAAAGAGCTGTGGCTAACCAGCCAGAGCCAGCCTAAAGCCTGGGACAAGTTATTTAGTGTCTGGCGGGCGGCCGGGAAACAAGATCCGCGGGCGTATTTTGAGCGGATCCggc                                                                                                 >  2:362586/1‑249 (MQ=255)
                                                                                                                                                        aaTGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTAGGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCt                                                                                                <  2:211951/251‑1 (MQ=255)
                                                                                                                                                                 gcgcGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGaa                                                                                       <  1:270817/251‑1 (MQ=255)
                                                                                                                                                                  cgcgGCGTGAAGCCTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGAAAGGGGCGGAAGAGCTGTGGCTAACCGGCAAGAGCCAGGCTAACCCCGGGGGCAAGGTATTTAATGTTCGGCGTGGGTCCGGTAAACAAGGTACGCTGGCGTCTTTTGATCGTAGCCgcgcagcgattaa                                                                                       >  2:85779/1‑238 (MQ=255)
                                                                                                                                                                                  gCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGcc                                                                      <  2:155982/251‑1 (MQ=255)
                                                                                                                                                                                    gTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTg                                                                    <  1:360772/251‑1 (MQ=255)
                                                                                                                                                                                      ggCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAATCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGt                                                                  <  1:109584/251‑1 (MQ=255)
                                                                                                                                                                                        cTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGCGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTa                                                                 >  1:310795/1‑250 (MQ=255)
                                                                                                                                                                                           gTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACa                                                              <  2:235764/250‑1 (MQ=255)
                                                                                                                                                                                            ttAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGt                                                            >  1:361080/1‑251 (MQ=255)
                                                                                                                                                                                                   ttAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggc                                                      <  2:361080/250‑1 (MQ=255)
                                                                                                                                                                                                     aGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagg                                                   >  1:154005/1‑251 (MQ=255)
                                                                                                                                                                                                         cGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagggc                                                 >  1:404282/1‑249 (MQ=255)
                                                                                                                                                                                                           gAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAt                                             <  1:177446/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccctgAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCGGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcctgc                                       <  2:34337/247‑1 (MQ=255)
                                                                                                                                                                                                                         cTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTAc                                >  1:319224/1‑251 (MQ=255)
                                                                                                                                                                                                                            ccGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAGGCCTGCCGATTAca                               >  2:458308/1‑248 (MQ=255)
                                                                                                                                                                                                                                 gcgc‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTAt                        <  1:40670/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGCTTACCAGACTATCGCCt                   >  1:31547/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ttACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCa              >  1:50377/1‑250 (MQ=255)
                                                                                                                                                                                                                                               actactAGGCAAAATGGAACCCCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAAt            <  2:58784/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        cGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTg  >  1:62487/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
TGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTG  >  NC_000913/4630899‑4631395

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