Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,978,164 T→A 23.4% Y245F (TAT→TTT)  lysA ← diaminopimelate decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,978,1640TA23.4% 95.7 / 6.0 53Y245F (TAT→TTT) lysAdiaminopimelate decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (26/14);  new base A (2/11);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 3.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAGAT  >  NC_000913/2977935‑2978402
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
ggCCCGGGCTCTGATAGGGTGGTAACCCACGTACCTTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCCCAAAGTGGCGGCCCCCCATTTGTTTGACGCTCCGCCCCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGtt                                                                                                                                                                                                                                       <  2:930691/238‑1 (MQ=255)
 tgtgtataaaaggcagggtggtAACTTCCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTGGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGtt                                                                                                                                                                                                                                       <  1:675618/222‑1 (MQ=255)
          ttgttttgggtttAACCACCCTACATTGCCGGGCCCATCGATTCGTTGACCCCGGCATCAACCCGCAAAAGTGGGCGGCCCCCCATTTGTTGGCCGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAAc                                                                                                                                                                                                                             <  2:1060647/227‑1 (MQ=255)
          attattttgttgtAACCACCGTATATTCGCGGGCCCTCCAGTTCGTTGAACCCGGCCACCACCAGCCAAAGGGGGCGGCCCCCCCTTTGTTTGACGCTCCGCCACGGAGTAATTAATACGCCAGACGGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACGGGGTGGCCAAAATGGCGGGCTATTTGCTCACCCGCGCCATCCCACAGACCATAATAATGTTCGGTATCAAc                                                                                                                                                                                                                             <  2:472045/228‑1 (MQ=255)
           ttattgtggggtAACCACCGGTCATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTc                                                                                                                                                                                                                                      <  2:1208501/220‑1 (MQ=255)
                  gtggtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGGTTTGCTCCCGCGCGGCATTCCACAGACCCTAATAAAGTTCGGTATCAACCGCccct                                                                                                                                                                                                                      >  1:463565/1‑238 (MQ=255)
                   tgtaacacaCCGTACATTGCGGGGCCCAGCAGTTCTTTGACCCCGCCCACACCCAACAAAAGGGGGCGTCTCCCATTTTGTTTGACGCCCCGCCCCTGATTAATTAATACGCAAGACTGTGCTACCAGGAAGGGCCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCAAAAGGGCGGGCTATTTGCTCCCGCGCGGCATTCCACAGACCAAAAAAATGTTGGGTACCAACCGCCTCtt                                                                                                                                                                                                                     <  2:1141459/230‑1 (MQ=255)
                     gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGGGATTTGCTCACACCCGGCATTCCCCAGACACAAATAAAGTTCGGGTTCAACCGctcgtcctt                                                                                                                                                                                                                  >  1:1208501/1‑231 (MQ=255)
                     gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAAGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAAAGTTCGGTATCAACCGCCTCTTCa                                                                                                                                                                                                                   >  1:1060647/1‑238 (MQ=255)
                        aCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTTACGCACCGCACCCGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGGTCAATTTTCAGTTTCCACGGGTGGCCCAAAAAGGCGGCACTTTGCTCACCCGCGGCCTTCCACAGCACATAATAAAGGTCgggttacaccaccccttaccccc                                                                                                                                                                                                               >  2:750115/1‑218 (MQ=255)
                        aCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTTCCAGGGTGGCCCAAATGGCGGGCGATTTGCTCTCGCGCGGCATTCCCCAGACCCTAATAATGTTCGGGATCAACCGCCCCCTCAc                                                                                                                                                                                                                  >  1:234078/1‑236 (MQ=255)
                            ccGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCCTTTGTTTGACGCTCCCCCCCCGAGAAATTAATACGCCCGACTGCGCTACCCGGAAGCGGCCCCGGTCCATTTCCCATTTTACAGGGGGGCCCAAAAAGGGGGCGGCTTTCTCTCCCGGCGCCTTTCACAAAACAAAATAAAATTCTGTGTCAAACCCCTCTTctttaaaac                                                                                                                                                                                                            >  2:675618/1‑230 (MQ=255)
                                                        tGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGGTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCCCAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCCGTTGATAAGGAACAGAAAGCCCAccgcccccc                                                                                                                                                                               >  1:186774/1‑236 (MQ=255)
                                                             ccGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCCCCCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCCGTTCAATTTCCAGTTTCACAGGGTGGCCCCAAAGGCGGGCGATTTGCTCACGCGCGGCGTTTCCCAGACCATAATAATGTTCGGTGTCCACCCCCCCTCCACCCTCTTGTTAAGGAAAAAAAAACACACCGCCCCCCGaaaa                                                                                                                                                                          >  2:1219285/1‑238 (MQ=255)
                                                                                                      tgACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAACCCACCGCCCCCAGAAATAGCCTGTAAATCCTGACCGAAATCGATGACCTGACGAAcc                                                                                                                                  >  2:1208694/1‑239 (MQ=255)
                                                                                                          gCTCCGACCCTGAGTATTTAATACGCCAGACTGCGCTACCAGGAGGCGACCCGGTTCAATTTCCAGTTCCACAGGGTGGCCCAAAGGGCGGGCGTTTTGCTCACGCGCGGCATCCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  1:1219285/239‑1 (MQ=255)
                                                                                                          ccTCCGCACCGAGTAAATTAATCCGCAAGACTGCGCTCCAAGGAGGCGCCCCGGTTCAATTCCCAGTTCCACGGGGTGGCCAAAAGGGCGGGCAATTTGCTCACGCGCGGCATTCCACAGCCAATAAAAATGTCGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  1:1208694/238‑1 (MQ=255)
                                                                                                             ttgtttaaagggaaagtaatCCGCCGGACGGCGCTCCCAGGAAGCGCCCGGGTTCAATTTCCATTTGCACAGGGTGGCCAAAATGGGGGGCGTTTTGCCCACGCGCGGCATCCCACAGACAAAAATAATGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  2:8517/221‑1 (MQ=255)
                                                                                                                gccttggcgttatttacaccccAGACTGCGCGACGAGGAGGCGACCCTTTTCAATTTCCGATGTCACAGGGTGCCCCAATGTGGGGGCGTTTTGCTCACTCGCGCTATTCCCAAGCCCATAATAGTTTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGCCCTGACGCACCATAGcacca                                                                                                                         <  2:337628/219‑1 (MQ=255)
                                                                                                                      aGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCt                                                                                                                  >  1:381135/1‑239 (MQ=255)
                                                                                                                      aGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATCCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGACAGCCCACCGCCCGCAGAAATAGCATGTAAATCCTGACCGAATTCGAGGACCTGACGCACCATAGcacctgaccctt                                                                                                                  >  1:1336188/1‑231 (MQ=255)
                                                                                                                        taattaatACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGTTAAGGAACAGAAAGCCCCCCCCCCCCAGAAAAAGCCTGTAAACCCTGACCGAATTCGGTGACCtcaacacccaaaacccccaaca                                                                                                                     >  2:1120278/1‑213 (MQ=255)
                                                                                                                               agccGCAAGGCTGCGCAACCAGGAAGCCACCCGGTTAAATTCCCAGTTTCCCGGGGGGGCCAAAATGGCGGGCAATTTGCTCACGCGCGGAATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                          <  1:1008958/236‑1 (MQ=255)
                                                                                                                                 aCGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCCTAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTAGATGAACTGACGCCCCCTACCACCACCCACCCGTTCCaacga                                                                                                        >  2:1008958/1‑235 (MQ=255)
                                                                                                                                     caagctgccttACCAGGAAGGCACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCAAAAGGGCGGGCAATTTGCCCACGCGCGGCATTCCACAGACAATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                                                                    <  1:1179127/229‑1 (MQ=255)
                                                                                                                                        aCTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACAccctgtcccttttacacat                                                                                                  >  1:8517/1‑221 (MQ=255)
                                                                                                                                            cgcTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCCGTTCCAGATGGGCATAATCa                                                                                              >  1:209642/1‑238 (MQ=255)
                                                                                                                                            cgcTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAAAAGCCTGTAAATCCTCACCGAATTCGATGACCCGACCCACCACAACCACACCACCCTGTTCCCGCTGGGCATAAtaa                                                                                              >  2:938049/1‑236 (MQ=255)
                                                                                                                                             gcTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACCCCTGTTCCAAATGGGCATAATCaa                                                                                             >  1:1270711/1‑238 (MQ=255)
                                                                                                                                                 ccgggagggcccccGTTTCAATTTCGAGTTTCCAGGGGGGGCCCAAAGGGGGGGGTGTTTGCGCCCGGGCGGATTTCAAAAGACAAAAATAATGTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCa                                                                                        <  2:755384/228‑1 (MQ=255)
                                                                                                                                                 ccAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACCCGCGGCATTCCACCGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAAAAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGCCCGAATTCGCTGACATGACCGCCCAAAGCACCACCCACCCGTTCCAAGTGGGCCAAAACAACCCCa                                                                                        >  2:1307120/1‑239 (MQ=255)
                                                                                                                                                      aaGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTACATCCTGACCGAATTCGATGACCTTACGCACCATAGCACCACACCCCCGTCCCAACTGGGCTTAATCAACGCCAAAAcc                                                                                   >  1:128334/1‑239 (MQ=255)
                                                                                                                                                         gcgcccGGTTTTATTTCCCATTTTCAGGGGGGGCCCAAAATGGCGGGGAATTTCCCCCCGCGCGGCATCCCACAGACCAAAATAATGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  2:381135/235‑1 (MQ=255)
                                                                                                                                                         gcccccGGTTAATTTTCCATTTCCCAAGGGGGCCCCAAAGGGGGGGGTTTTTGCCCACCCGCGGCATCCCACAACCAAAAAAAAGGTGCGGCATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  2:1270711/234‑1 (MQ=255)
                                                                                                                                                         agACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTTCACAGACCATAATAATGTTCGGTATCCACCGCCTCTTTACCCTGTTGATAGGGAACAGAAAGCCCCCCGCCCGCAGAAATAGCCAGTAAATCCTGACCGACTTCGCTGCCCTTACCCACCATTGCACCACCCACCTTTTCCAACTGGGGTTATAAAACGCCCGAACCacc                                                                                >  2:1036991/2‑237 (MQ=255)
                                                                                                                                                          gACCCGGTTTCAATTTCCAGTTTAACAGGGTGGCCAAAATGGCGGGCGATTTGCTCACGCGGGGCATTCCCAGGCACATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  1:938049/237‑1 (MQ=255)
                                                                                                                                                           aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCCCACACCTGTTCCAGATGGGCATAAt                                                                                                >  1:16505/1‑221 (MQ=255)
                                                                                                                                                              cGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCCCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCCCACACCGGTTCCAGCTGGGCATAATAAACGCCAGAACCAAAGTGCa                                                                           >  1:552063/1‑239 (MQ=255)
                                                                                                                                                                        ttcccatgtccacaGGGTGGCCAAAAGTGCTGGCGATTTTCTCACGCCACGCATCCCACAGACACAAATGATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATg                                                                   <  1:672938/229‑1 (MQ=255)
                                                                                                                                                                                           gCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCacacac                                                                                                                    >  1:272719/1‑171 (MQ=255)
                                                                                                                                                                                                                  cACGCGCGGCATGCCACAGCCCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTgcg                         <  1:305909/239‑1 (MQ=255)
                                                                                                                                                                                                                  cACGCGCGGCATGCCAAAGACCATAATAAGTTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTgcg                         <  1:997651/239‑1 (MQ=255)
                                                                                                                                                                                                                            ctaccACAGACATAAAGAATGTTCGGTACAAACCGCCTCTTCACCCTGTTGAAAAGGAACAGAAAGCCCCCCGCCCGAAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  2:683485/236‑1 (MQ=255)
                                                                                                                                                                                                                            catCCACAGACCATAAAAATGTCCGGTAACAACCGCCTCTTCACCCTGTTGAAAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  1:113186/237‑1 (MQ=255)
                                                                                                                                                                                                                             agaaaagggccaataaaaATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  1:470486/227‑1 (MQ=255)
                                                                                                                                                                                                                                cccaGACAAAAATATTGTGCGGTATCACCCCCCTCTTCCCCCTTTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  2:16508/237‑1 (MQ=255)
                                                                                                                                                                                                                                ccaaGACAAAAATAATGTTCGTTAACAACCGCCTTTTCACCCTGGTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:1052968/236‑1 (MQ=255)
                                                                                                                                                                                                                                agacagccATAATAATGTTCGGTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:912065/233‑1 (MQ=255)
                                                                                                                                                                                                                                acaaGACAAAAATAATTTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:750115/236‑1 (MQ=255)
                                                                                                                                                                                                                                     aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCACTCTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGGGCGGCGGGAACATCGGGGTAcc      >  1:1354119/1‑239 (MQ=255)
                                                                                                                                                                                                                                     aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGAACTAGTATCACGTCCAATGCGCCGGCCaaggccgcgtacc      >  1:442342/1‑227 (MQ=255)
                                                                                                                                                                                                                                     aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACCACAGAACCAATGTGCATGGGAATGCCGACCAGCTTCAGATGATGACGTTGTATCACGTCCAGAGCGGCgggacaaaagggggacc      >  2:499725/1‑225 (MQ=255)
                                                                                                                                                                                                                                     aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGGTCCAGATGGGCATAATCAACGCCAGAACCAATTGGCATGTGAATCCCGACCAGCTGCAGCGGATTACGTTGTATAACGCCCAGGGAGGGGGGCGGagagggggagc      >  2:96286/1‑229 (MQ=255)
                                                                                                                                                                                                                                        cATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCAAACACCTGTTCCCATTGGGCATAATCAACGCACGAACCAATTTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCCGATCGGGGTACCaa    >  2:1060460/1‑238 (MQ=255)
                                                                                                                                                                                                                                          taataaTGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCGGTGCGGCGGGAAGATCGGTGTACCaccg  >  1:683485/1‑236 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
TGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAGAT  >  NC_000913/2977935‑2978402

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