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,1460GA97.5% 118.2 / NA 40A138A (GCG→GCAsltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/22);  minor base C (1/0);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630913‑4631387
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGcc                                                                                                                                                                                                                                  >  1:131715/1‑250 (MQ=255)
      ccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAAGTACTAGTATGCAAAATGGATCACCGTGCAGAGTg                                                                                                                                                                                                                             <  2:27585/247‑1 (MQ=255)
       cccTTTATCCCTACCTGGAATACCGTCAGATTACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAGGGAACACCGGGCAGAGTgaa                                                                                                                                                                                                                           <  2:26910/249‑1 (MQ=255)
               cccTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCt                                                                                                                                                                                                                    <  1:72570/250‑1 (MQ=255)
               cccTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTg                                                                                                                                                                                                                   <  2:215121/251‑1 (MQ=255)
                  tACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGCCTGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCACTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGGCCGGAAAAGCCCGGAACTACCGGAGCCCAAACTAATTACTACTTAGCCAAATGGAAAAACGGGCAAGGGGAAGAAGCCTGGcc                                                                                                                                                                                                                >  1:115941/1‑250 (MQ=255)
                     cTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGAAAAGCCCGGAACTAACGAAGCGCAATGTAATTTATACTATGCAAAAAGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgg                                                                                                                                                                                                             >  1:115434/1‑251 (MQ=255)
                     cTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGGACTTCCGAAGCCCCATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGCAGAAGACTgggcagg                                                                                                                                                                                                             >  2:96552/1‑246 (MQ=255)
                         aaTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGAGTGATGACTGACGTGGCTTGTTAGGCTTTAGGCCCGGAAAGCCCCGAACTCCCGAAGCGCAAAGCAATTACCACTCTGCAAAATGGAACACAGGGAAGAGTTAAGAAGCCTGGCCAGGGGCg                                                                                                                                                                                                         >  1:23524/1‑251 (MQ=255)
                                 cAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAATGGGCGAAAGAGCt                                                                                                                                                                                                 <  2:164754/251‑1 (MQ=255)
                                                          aaCCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCt                                                                                                                                                                        <  1:202591/251‑1 (MQ=255)
                                                              ggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAAACGGCAAGAGCCAGCCTAACg                                                                                                                                                                    >  1:193259/1‑251 (MQ=255)
                                                              ggcggTGACTGTCACTAACTTTGTTCGCAATAACCCCACGCTTCCTCCCGCTCGCAAGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAACAACCCGACCTACCGAAGCGCAATGTAATTACTCCTATGCAACATGGAACACCGGGCACAGTGAAGAAGACTGGCAAGGGGGGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACg                                                                                                                                                                    >  2:12473/1‑251 (MQ=255)
                                                                  gTGACTGTCACTAACGTTGTTCGCGCTAACCCCACGCTTCCTCTCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGAAAAGCCCGGATCTACCGAAGCGCAATGTAATGCCTACTATGCAAAATGTAACACCGGGCAGAGGGAAGAAGCCTGGCAAGGGGCGAAAGCGCTGTGGCTAACCGGCAAGAGACAGTCTAACGCCtg                                                                                                                                                                <  2:207464/251‑1 (MQ=255)
                                                                     actgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCtgtgg                                                                                                                                                             >  2:61646/4‑250 (MQ=255)
                                                                         tCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCACTGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGACTTTATCCCGGAAAAGCCCGGAACTACCGAAACGCAATGTGATTACTACTATGCAAAATGGAACACCCGGGCGAGCGAAGAAGCCTGGCAAGGGGCGAAAGGGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCtgtgcca                                                                                                                                                           >  1:1679/1‑246 (MQ=255)
                                                                         tCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACa                                                                                                                                                           >  1:92138/1‑249 (MQ=255)
                                                                                        gcgcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCTTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAGGCGCAATTTAATTACTACTAGGCAAAATGGAACACCGGGAAGGGGGAAGAAGCCTGGCAAGGGGCCAAAGAGCTGTGGGTAACCGGCAAGAGCCAGCCTAAAGGCTGTGACAAGGTATTTAGGGTCTgg                                                                                                                                          >  2:200962/1‑251 (MQ=255)
                                                                                          gcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGc                                                                                                                                                                                    <  1:5987/207‑1 (MQ=255)
                                                                                                     gCTTCCTCCCGCTCGAACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGCAGCCTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGCACTACCGAAGCGCAATGTTATTACTACTATGCAAAATGGACACCCGGGCAGCGTTGAGAAGCCTGGCCAGGGGCGAAAGAGCTGTGGCTAACCGGCACGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGGCTGGCGTGCGTCGGGt                                                                                                                              >  2:134505/1‑250 (MQ=255)
                                                                                                            cccGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGa                                                                                                                      <  1:114333/251‑1 (MQ=255)
                                                                                                                         tGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAGTGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGGCAGGGAAACAAGATCCGCTGGCGTAt                                                                                                         >  1:228635/1‑251 (MQ=255)
                                                                                                                          gCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTAt                                                                                                         <  1:141140/250‑1 (MQ=255)
                                                                                                                          gCAATCGCGGTTCGGCAATGAACGGGCGGGGAGTGAAGACTGGCGTGGCTTGTGAGCCTTTAGCCCGTAAAAGCCCTGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGTGCAGAGTTAAGAAGCCTGGCAATGGGCGAAAGAGCTGTTGCTAACCGGCAAGAGCCAGCTTAACGTCTGTGACTATTTAGTTATTGGGTGGCGTGCGTCAGGTAAACAAGATCCGCTTGAGTAtt                                                                                                        <  2:243837/251‑1 (MQ=255)
                                                                                                                           cAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTAttt                                                                                                       <  2:115434/251‑1 (MQ=255)
                                                                                                                                             gAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGc                                                                                     <  1:12473/251‑1 (MQ=255)
                                                                                                                                                     gcggcgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAg                                                                             <  1:181235/251‑1 (MQ=255)
                                                                                                                                                         cgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACCCCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGGAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCCGTGACAAGTTATTTAGAGGCAGGCGTGCGGCAGGTAAACCAGATACGCTGGCGAATTTAGCGCGTATCCGACTGGCGGTGAAAGAggg                                                                         >  2:192684/1‑251 (MQ=255)
                                                                                                                                                                     cGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCt                                                             <  1:200962/251‑1 (MQ=255)
                                                                                                                                                                            tGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAAc                                                       >  1:155125/1‑250 (MQ=255)
                                                                                                                                                                                     ttAGCCCGGAAAAGCACGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGTTAACAGTGCTgg                                               <  1:62728/249‑1 (MQ=255)
                                                                                                                                                                                       aGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTTCTggcagg                                           <  2:131715/251‑1 (MQ=255)
                                                                                                                                                                                              gAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAt                                     <  1:62177/251‑1 (MQ=255)
                                                                                                                                                                                                aaaGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGAAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgc                                   >  1:238221/1‑251 (MQ=255)
                                                                                                                                                                                                       cGGAACTACCGAAGCGCAATGTAATTACTCCTATGCAAAATGGAACACCGGGTAGAGTGAAGAAGCCTGCCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGTCTGTGACTAGTGATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTACCACGGGCCTGGTAAGCGTGCTGGCAAGCCAGATGTCTGCCg                             <  2:109958/251‑1 (MQ=255)
                                                                                                                                                                                                                      gcAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCg              >  1:5711/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaTGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGcc            >  1:87816/1‑251 (MQ=255)
                                                                                                                                                                                                                           gtcaTTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGTCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGTCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGTTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTc          <  1:119124/247‑1 (MQ=255)
                                                                                                                                                                                                                           gTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTc          <  1:47773/250‑1 (MQ=255)
                                                                                                                                                                                                                                    actaTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCa  <  2:193259/249‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
GATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630913‑4631387

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