Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,960 T→C D76D (GAT→GAC slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9600TC98.2% 178.8 / NA 54D76D (GAT→GACsltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/30);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGA  >  NC_000913/4630734‑4631193
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
gggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGAct                                                                                                                                                                                                                      <  2:346478/250‑2 (MQ=255)
  gAAAAAGCCAAACAATGTTCCTGGGGGTTGTTGGCGGCCGGGGGCTGTCTGCTGGAGGTCACCAGCGTGGCTCGAGCCGGCCCACTGGATTAGGAGCGTAGCCGTTTCTCGCAAAACAAGCAGGCCTGGGATAATCGACAAATGGATGTTGTCGAAAAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGActgt                                                                                                                                                                                                                    <  2:369971/251‑4 (MQ=255)
          cAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACt                                                                                                                                                                                                            >  2:103036/1‑251 (MQ=255)
           aaaCAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACtt                                                                                                                                                                                                           >  1:267966/1‑251 (MQ=255)
            aaCAATTTCCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACttt                                                                                                                                                                                                          <  1:292907/251‑1 (MQ=255)
                        gggcggcTGTGGGCTGCCGGTGTCTGTCTGCGGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTaa                                                                                                                                                                                               <  2:84010/250‑1 (MQ=255)
                           cggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAccc                                                                                                                                                                                            <  1:265696/251‑1 (MQ=255)
                                gTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCAc                                                                                                                                                                                         <  1:33570/249‑1 (MQ=255)
                                     cTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTACCCCCACGCTtcct                                                                                                                                                                                  >  2:250964/1‑251 (MQ=255)
                                       gCCGGTGTGTGTCCGCGGACGGTCAGCTGCTTGGCGCGTGCCGACTTACGGTAGTAGCAGCGTAGGCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAACGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTtcctcc                                                                                                                                                                                <  2:333/251‑1 (MQ=255)
                                         cGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCAATAACTTTGTTCGCGCTAACCCCACGCTTCCTcccc                                                                                                                                                                              >  2:473292/1‑250 (MQ=255)
                                           gtgtCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCt                                                                                                                                                                            <  1:41574/251‑1 (MQ=255)
                                            tgtctgtcTGCTGACGGTCAGCAGCGTGGCGCGCGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGACTACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCCACCGCTc                                                                                                                                                                           >  2:454108/1‑251 (MQ=255)
                                            tgtctgtcTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCCAGCAGGCATGGGATAATCGACAAATGGATGTGGTCGCACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACCTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTc                                                                                                                                                                           >  2:173658/1‑251 (MQ=255)
                                                 gtctgCTGACGGTCAGCAGCGTGGCGCGAGCCTACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCAc                                                                                                                                                                       <  1:365673/250‑1 (MQ=255)
                                                                 cagcGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGCTTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGACTCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCCCCCGCCCGCTGCCATCTCGTTTCg                                                                                                                                                      >  2:301964/1‑251 (MQ=255)
                                                                             aGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAACTCAAGCAGGCCTGGGCTAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCCCGCTAACCCCACGCATCCTCCCGCCCGCCCGCTGCAATCTCGTTTCGTCAATggaatgg                                                                                                                                          >  2:482203/1‑245 (MQ=255)
                                                                                      aCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTg                                                                                                                                 <  2:275164/251‑1 (MQ=255)
                                                                                         ggATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAg                                                                                                                              <  1:243338/251‑1 (MQ=255)
                                                                                                       gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                                                <  2:214617/249‑1 (MQ=255)
                                                                                                       gccGTTACGCGCAAATCAAGCAGGCCCGGGACAATCGACAACTGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                                                >  1:374354/3‑251 (MQ=255)
                                                                                                         cggtACGCGCAAATCAAGCAGGCCTGGGATATCTGATAAATGGATGTGGTCTAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTGCCTTGAATACCGTCAGATCACCGATGACCTGATGAATCATCCGGCGGTGACTGTCACTAACTCTGTTCGCGCTAACCCCACTCTTACTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGtt                                                                                                              <  1:386477/248‑1 (MQ=255)
                                                                                                            tACGCGTAAATCAAGTAGGTCTGGGAGAATAGACAAATGGATGTGGTCGAACAAATGATGCCTGGATTGAAAGATTATACACTTTATCCATTCCTGGAATATAGTCAGATAAACGATGATCTGATGAATCAACCGGTGGTGACTGTCACTAACTTTGTGCGTGCTAAACCCACGCTTCCGCCCGCCCGCACGCTGGAATCTAGTTTTGTCAATGAATGGGAGCGGCGTGAAGACTGGCGTGGCTTggtgg                                                                                                            <  1:303914/250‑5 (MQ=255)
                                                                                                             aCGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGcc                                                                                                         >  1:79578/1‑251 (MQ=255)
                                                                                                                   aaaTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTTGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGc                                                                                                   <  2:474894/251‑1 (MQ=255)
                                                                                                                          gCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGaaa                                                                                           <  2:420723/251‑1 (MQ=255)
                                                                                                                          gCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAAGGAACTGGCGCGGCGTTAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGaaa                                                                                           <  1:283020/251‑1 (MQ=255)
                                                                                                                               ccTGGTATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTTGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATTAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGccc                                                                                      <  1:175944/251‑1 (MQ=255)
                                                                                                                                 tGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAAACCCgg                                                                                    >  1:428999/1‑251 (MQ=255)
                                                                                                                                    gATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGAATGGCTGGGCTTTTTTGCCTTTAGCCCGGGAAAAGccggaa                                                                                   >  2:369171/1‑246 (MQ=255)
                                                                                                                                        aTCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTAcc                                                                            <  2:231723/251‑1 (MQ=255)
                                                                                                                                            gACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGc                                                                                                                   <  1:138952/212‑1 (MQ=255)
                                                                                                                                                 aTGGATGTGGTCGAACAAATGATTCCTGGACTTAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGACCTCCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCa                                                                    <  2:334330/251‑1 (MQ=255)
                                                                                                                                                   ggATGTGGTCGACCAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCaa                                                                   >  1:354689/1‑250 (MQ=255)
                                                                                                                                                         ggTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCCATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAggccaatgat                                                               >  1:473761/1‑239 (MQ=255)
                                                                                                                                                         ggTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCGGGAAGACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTa                                                               <  2:473761/248‑1 (MQ=255)
                                                                                                                                                             gAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATtacta                                                        <  1:418181/251‑1 (MQ=255)
                                                                                                                                                                 aaaTGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGGATGGCGTGGCTTGTTACACCTTTAGCCGGGAAAGCCCCGAGCTACCGCAGAGCACTGTAATTTCTACTAt                                                    >  2:77886/1‑251 (MQ=255)
                                                                                                                                                                     gatgCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAACCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCCTTTTAGCCTTTAACCCGGGAAAGCCCGGAACCAACGAAACGCAAATTAATTACTACTAttca                                                 >  1:436007/1‑247 (MQ=255)
                                                                                                                                                                            ggACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGGACTACCGAAACGCAATGTAATTACTACTATGCAAAATGGa                                          >  1:364418/1‑250 (MQ=255)
                                                                                                                                                                               cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTAGCCTTTAGCCCGGGAAAGCCCGGAACTACCGAAGCGCCATGTAATtagatcgggagagcaccagtctg                                      >  1:254520/1‑230 (MQ=255)
                                                                                                                                                                               cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAAtt                                                            <  2:254520/229‑1 (MQ=255)
                                                                                                                                                                                tGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGGGTGGGTTGTTAGCCTTTTGCCCGGAAAAGCCCGGGAACACAGAAGCGCAATGTAATtaactataagcaaaaaggagcacc                                     >  2:156723/1‑229 (MQ=255)
                                                                                                                                                                                     gATTATCCACTTTATCCCTACCTGGAATACCGTCAGTTCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGtt                                                                                                              >  1:171970/1‑176 (MQ=255)
                                                                                                                                                                                     gATTATCCACTTTATCCCTACCTGGAATACCGTCAGTTCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGtt                                                                                                              <  2:171970/176‑1 (MQ=255)
                                                                                                                                                                                           ccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGTAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGTAACACCGGGCAGAGTg                           <  1:118937/247‑1 (MQ=255)
                                                                                                                                                                                           ccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTGGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGTGTga                          <  2:273105/248‑1 (MQ=255)
                                                                                                                                                                                            cacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTgaa                         <  2:354689/249‑1 (MQ=255)
                                                                                                                                                                                                 tATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCCTTAGCCCGGCAAAGCCCGGGACTACCGGAGCGCCATGGAATTTCTACTATGCCAAATGGGAAACCGGGcgacgtgaaga                       >  2:88515/1‑238 (MQ=255)
                                                                                                                                                                                                    cccTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCt                  <  1:369171/250‑1 (MQ=255)
                                                                                                                                                                                                            ggAAAACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgg           <  2:267966/249‑1 (MQ=255)
                                                                                                                                                                                                              aaTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAAGTAATTACTACTATTCAAAATGGAACACCGGGAAGAGTGAAGAAGCCTGGCAAGGGGCg       >  1:325486/1‑251 (MQ=255)
                                                                                                                                                                                                                   cgtcAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGTGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAaga  <  1:103036/248‑1 (MQ=255)
                                                                                                                                                                                                                           cACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAAGCCCGGGACTACCGAAGCGCAATGt                                                                >  1:159356/1‑181 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
TGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGA  >  NC_000913/4630734‑4631193

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