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→GCAsltsoluble 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

TGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTG  >  NC_000913/4630899‑4631395
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
ggCCTGGACTGAAAGATTAGCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCTCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACGTGCGCGGCGTGAATACTGGCGTGGCTTGTTAGTCTTGAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCaaaa                                                                                                                                                                                                                                                        <  1:62715/250‑1 (MQ=255)
  ccTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGAACTACCGGAGCGCCAATTAATTACTACTATGCAAAATg                                                                                                                                                                                                                                                      >  1:63547/1‑251 (MQ=255)
   cTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGGCTTTAGCCCGGAAAAGGCCAGCACTACCGGAGCGCACTGTAATTACTACTATGCCAAATgg                                                                                                                                                                                                                                                     >  2:40670/1‑251 (MQ=255)
      gACTGAAAGATTATCTACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTACCCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGTCTGGCGTGGCTTGTTAGCCTTTAGCTCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAac                                                                                                                                                                                                                                                  <  2:455632/251‑1 (MQ=255)
      gACTGAAAGATTATCCCCTTTATCCCTACCTGGAATTCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAGGGAac                                                                                                                                                                                                                                                  <  2:455641/251‑1 (MQ=255)
      gACTGAAAGATTATCCACTTTATCCCTACCTTGAATACCGTCAGATCCCCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACCTTGTTCGCGCTAACCCCACGCTTCCGCTCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAac                                                                                                                                                                                                                                                  <  1:364007/251‑1 (MQ=255)
      gACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTTGCCCGGAAAAACCCCGGAATACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAac                                                                                                                                                                                                                                                  >  2:177446/1‑251 (MQ=255)
         ggaaagATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGTCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACAcc                                                                                                                                                                                                                                               <  2:471795/246‑1 (MQ=255)
               aTTATCCACTTTATCCCTACCTGGCATACAGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGAGAAGACTGGCGTGGCTTGTTAGACTTTAGACCGGAAAAGCCCGGGACTACCGAAGCGGAAAGTAATTAATACTATGCAAAATGGAAAACCgggg                                                                                                                                                                                                                                           >  1:357877/1‑248 (MQ=255)
                         ttATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTgaagaa                                                                                                                                                                                                                                <  1:262683/250‑1 (MQ=255)
                          tATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGaagaag                                                                                                                                                                                                                               <  1:21120/250‑1 (MQ=255)
                            cccctccctGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACGACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCt                                                                                                                                                                                                                            <  2:452649/250‑1 (MQ=255)
                                   ccggaaacacttcAGATCACCTATGACCTGAGGAACCAACCGGCGGTGCCTGTCACTAACTTTGTTCGCGCTAACCCACCGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATTAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgg                                                                                                                                                                                                                     <  2:119721/239‑1 (MQ=255)
                                                aGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGGAAAACCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAAAAGCCTGGCCAGGGGCGAAagaggtg                                                                                                                                                                                                        >  1:270689/1‑248 (MQ=255)
                                                          tgaccTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTAACGAAGCGCAAAGTAATTACTACTATGCAAAATGGAACACCGGGACGAGGGAAAAAGCCTGGCACGGGGCGAAAGAGCAGTGGCCAACCg                                                                                                                                                                                              >  1:494586/5‑251 (MQ=255)
                                                               tgatgaATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACACTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAAATTGCGCGGCGTGAAGACTGGCGTTGCTTGTTAGCCTTTAACCCGGAAAAGCCCGGAACTACCCAAGCGCAATGTAATTTCAACTTTGTAAAATGGAACACCGGAGAGAGTGAAAAAGCCCTGCAAGGGGGGaaaaacattgggcaaacctgcca                                                                                                                                                                                          >  1:191645/1‑230 (MQ=255)
                                                                gatgaATCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAAGCCCGGAACTAACGAAGCGCAATGTAATTACTAATCAGCACAATGGAACACCGGGCAGAGTGAAGACACCTGGCAAGGGGCGAAAGAGCTGTGGGTAACCgggcaa                                                                                                                                                                                         >  1:196000/1‑246 (MQ=255)
                                                                 atgaATCAACCGGCGGTGACTGTCACTAACGTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAagag                                                                                                                                                                                       <  1:32720/251‑1 (MQ=255)
                                                                           cggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAAc                                                                                                                                                                             <  1:179769/251‑1 (MQ=255)
                                                                              cggTGACTGTCACTATCTTTGTTCGCTCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCTTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCATGGGGCGAAAGAGCTGTGGCTAACCGGCAAGATCCCGCCTAACGcc                                                                                                                                                                          <  2:499609/251‑1 (MQ=255)
                                                                                  gactgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCtgtg                                                                                                                                                                      <  1:261614/247‑1 (MQ=255)
                                                                                      gTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGAc                                                                                                                                                                    <  2:396692/249‑1 (MQ=255)
                                                                                               tttGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGGGACAAGTTAttt                                                                                                                                                           >  1:34337/1‑249 (MQ=255)
                                                                                                 tGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAgt                                                                                                                                                        <  1:118450/250‑2 (MQ=255)
                                                                                                  gTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTATCCTTTAGCCCGGCAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAAAGCTGTGGCTAACCGGCAAGAGCCAGCCGAACGCCTGGGACAAGTTATTTAgggt                                                                                                                                                      >  2:123835/1‑248 (MQ=255)
                                                                                                          tAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGgcgtg                                                                                                                                              >  1:58784/1‑251 (MQ=255)
                                                                                                             ccccACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGCACAGCCAGGAACTCCCGCAACGCAATGTAAATACTACTATGCAAAATGGAAAACCGGGCAGAGTGCAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCCGCCTAACGCCTGTGACAAGTTATTTAGGGTCTGgcgtgcgt                                                                                                                                           >  1:326873/1‑251 (MQ=255)
                                                                                                                  cGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                                      <  1:462864/251‑1 (MQ=255)
                                                                                                                   gCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGACGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGGAAAGCCCCGAACTACCGCAGCACCATGTAATTACTACTATGCAAAAAGGGAAACCGGGCAGAGTGAAGAAGCCCGGAATGGGGCGAAAGAGCTGTGGCAAACCAGCAAATGCCCACCTAACGCCTGTGTCAAGTGATTTAGAGACTGgcgtgcggccggg                                                                                                                                      >  2:192812/1‑244 (MQ=255)
                                                                                                                      tcctccCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCCAAAAGGAACACCGGGCCGAGTGAAGAAGACTGGCACGGGGAGAAAGAGCTGTGGCTAACCGGCAAGAGCCCGCCCAACGCCTGTGACAAGTTATTTAGTGTCTGgcgggggtaaggaaaac                                                                                                                                  >  2:457390/1‑237 (MQ=255)
                                                                                                                      tcctccCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTaaa                                                                                                                                   <  1:166143/250‑1 (MQ=255)
                                                                                                                       cctccCGCACGCACGCTGCAATCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACa                                                                                                                                 <  1:229383/251‑1 (MQ=255)
                                                                                                                        ctccCGCTCGCACGCTGCAATATCGTTTCGTCAATGAACTGGCGCGGCGTGAAGTCTGGCGAGGCTTGTTAGCCTTTAGCCCGGCAAAGGCCGGGACTACCGAAGCGCCATGTAATTTCTACTATGTAAAATGGAACCCCGGGCAGGGTGAAGAAACCTGGCAAGGGGCGCAAGGGCTTGGGCTAACCGGCCAGAGCCCGCCAAAAGCCCTTGAATAGTTATTTAGTGTCTgggcggccgccagtcacca                                                                                                                                 >  2:230921/1‑233 (MQ=255)
                                                                                                                              cTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCg                                                                                                                          >  1:303836/1‑251 (MQ=255)
                                                                                                                                  cACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTCGCCTTTAGCCCGGAAAAGCCCGGCACAACCCGAGCGCAATGTAATTAATACTCTGCAACATGGCACCCCGGGCAGAGGGAAGAAGCCTGGCAAGGCGCGAAAGAACTGTGGGTACCCCGGAAGAGCCAGCCTCACGACTGTGACACGTTATTTAATGTCTTGGGTGCGTCAGGTCAACAAGTTcccgtgg                                                                                                                      >  1:34970/1‑246 (MQ=255)
                                                                                                                                      cTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCCTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGCAAGGGGCGAAGAAACTGTGGCTAACCGGCAAGAGCCAGCCTACCGACTGGGACAAGTTTTTTAGGGGCTGGCGGGCGGCAGGGAACAAAGATCCAGTGGCGtg                                                                                                                  >  2:346872/1‑250 (MQ=255)
                                                                                                                                         cAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAACCTTTAGCCCGGAAAAGCCCGGAACTACCCAAGCGCAATGTAATTACTACTATGCAAAATTGGACACCGGGCAGAGGGAAATAGCCTGGCAAGGAGCGAAAGAGATGTTGGTAACAGGCAAGAGCCAGCCCAAAGACTGTGAAAAGTTATTTAGTGTATGGGGTGCGGCAGGTAAACAAGATCAGCTGGCGtcatt                                                                                                               >  1:363077/1‑247 (MQ=255)
                                                                                                                                          aaTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGGGCAATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGGCTGGCCAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCCGCCTAACGCCTTGGACAAGTTATTTAGTGTCTGGCTGGCGTCAGGGAAACAAGGTCCGCTGGccgtttt                                                                                                               >  2:359457/1‑244 (MQ=255)
                                                                                                                                               cGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGc                                                                                                          <  1:322175/250‑1 (MQ=255)
                                                                                                                                               cGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGGGCGTCAGGTAAACAAGATCCGCTGGCGTATTTaga                                                                                                            <  2:462772/248‑1 (MQ=255)
                                                                                                                                                       cAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCGGGCCAGGGGCGAAAGAGCTGTGGCTAACCAGCCAGAGCCAGCCTAAAGCCTGGGACAAGTTATTTAGTGTCTGGCGGGCGGCCGGGAAACAAGATCCGCGGGCGTATTTTGAGCGGATCCggc                                                                                                 >  2:362586/1‑249 (MQ=255)
                                                                                                                                                        aaTGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTAGGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCt                                                                                                <  2:211951/251‑1 (MQ=255)
                                                                                                                                                                 gcgcGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGaa                                                                                       <  1:270817/251‑1 (MQ=255)
                                                                                                                                                                  cgcgGCGTGAAGCCTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGAAAGGGGCGGAAGAGCTGTGGCTAACCGGCAAGAGCCAGGCTAACCCCGGGGGCAAGGTATTTAATGTTCGGCGTGGGTCCGGTAAACAAGGTACGCTGGCGTCTTTTGATCGTAGCCgcgcagcgattaa                                                                                       >  2:85779/1‑238 (MQ=255)
                                                                                                                                                                                  gCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGcc                                                                      <  2:155982/251‑1 (MQ=255)
                                                                                                                                                                                    gTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTg                                                                    <  1:360772/251‑1 (MQ=255)
                                                                                                                                                                                      ggCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAATCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGt                                                                  <  1:109584/251‑1 (MQ=255)
                                                                                                                                                                                        cTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGCGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTa                                                                 >  1:310795/1‑250 (MQ=255)
                                                                                                                                                                                           gTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACa                                                              <  2:235764/250‑1 (MQ=255)
                                                                                                                                                                                            ttAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGt                                                            >  1:361080/1‑251 (MQ=255)
                                                                                                                                                                                                   ttAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggc                                                      <  2:361080/250‑1 (MQ=255)
                                                                                                                                                                                                     aGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagg                                                   >  1:154005/1‑251 (MQ=255)
                                                                                                                                                                                                         cGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagggc                                                 >  1:404282/1‑249 (MQ=255)
                                                                                                                                                                                                           gAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAt                                             <  1:177446/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccctgAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCGGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcctgc                                       <  2:34337/247‑1 (MQ=255)
                                                                                                                                                                                                                         cTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTAc                                >  1:319224/1‑251 (MQ=255)
                                                                                                                                                                                                                            ccGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAGGCCTGCCGATTAca                               >  2:458308/1‑248 (MQ=255)
                                                                                                                                                                                                                                 gcgcAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTAt                        <  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)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
TGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTG  >  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: