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

TGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAGAT  >  NC_000913/2977935‑2978402
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
ggCCCGGGCTCTGATAGGGTGGTAACCCACGTACCTTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCCCAAAGT‑GGCGGCCCCCCATTTGTTTGACGCTCCGCCCCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGtt                                                                                                                                                                                                                                       <  2:930691/238‑1 (MQ=255)
 tgtgtataaaaggcagggtggtAACTTCCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTGGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGtt                                                                                                                                                                                                                                       <  1:675618/222‑1 (MQ=255)
          ttgttttgggtttAACCACCCTACATTGCCGGGCCCATCGATTCGTTGACCCCGGCATCAACCCGCA‑AAAGTGGGCGGCCCCCCATTTGTTGGCCGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAAc                                                                                                                                                                                                                             <  2:1060647/227‑1 (MQ=255)
          attattttgttgtAACCACCGTATATTCGCGGGCCCTCCAGTTCGTTGAACCCGGCCACCACCAGCCAAAGGG‑GGCGGCCCCCCCTTTGTTTGACGCTCCGCCACGGAGT‑AATTAATACGCCAGACGGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACGGGGTGGCCAAAATGGCGGGCTATTTGCTCACCCGCGCCATCCCACAGAC‑CATAATAATGTTCGGTATCAAc                                                                                                                                                                                                                             <  2:472045/228‑1 (MQ=255)
           ttattgtggggtAACCACCGGTCATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTc                                                                                                                                                                                                                                      <  2:1208501/220‑1 (MQ=255)
                  gtggtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGGTTTGCTCCCGCGCGGCATTCCACAGAC‑CCTAATAAAGTTCGGTATCAACCGCccct                                                                                                                                                                                                                      >  1:463565/1‑238 (MQ=255)
                   tgtaacacaCCGTACATTGCGGGGCCCAGCAGTTCTTTGACCCCGCCCACACCCAACAAAAGGG‑GGCGTCTCCCATTTTGTTTGACGCCCCGCCCCTGATT‑AATTAATACGCAAGACTGTGCTACCAGGAAGGGCCCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCAAAAGGGCGGGCTATTTGCTCCCGCGCGGCATTCCACAGAC‑CAAAAAAATGTTGGGTACCAACCGCCTCtt                                                                                                                                                                                                                     <  2:1141459/230‑1 (MQ=255)
                     gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGGGATTTGCTCACACCCGGCATTCCCCAGAC‑ACAAATAAAGTTCGGGTTCAACCGctcgtcctt                                                                                                                                                                                                                  >  1:1208501/1‑231 (MQ=255)
                     gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAAGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAAAGTTCGGTATCAACCGCCTCTTCa                                                                                                                                                                                                                   >  1:1060647/1‑238 (MQ=255)
                        aCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTTACGCACCGCACCCGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑GTCAATTT‑TCAGTTTCCACGGGTGGCCCAAAAAGGCGGCACTTTGCTCACCCGCGGCCTTCCACAGCA‑CATAATAAAGGTCgggttacaccaccccttaccccc                                                                                                                                                                                                               >  2:750115/1‑218 (MQ=255)
                        aCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCCCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTTCCAGGGTGGCCCAAATGGCGGGCGATTTGCTCTCGCGCGGCATTCCCCAGAC‑CCTAATAATGTTCGGGATCAACCGCCCCCTCAc                                                                                                                                                                                                                  >  1:234078/1‑236 (MQ=255)
                            ccGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCCTTTGTTTGACGCTCCCCCCCCGAGA‑AATTAATACGCCCGACTGCGCTACCCGGAAGCGGCCCCG‑GTCCATTT‑CCCATTTTACAGGGGGGCCCAAAAAGGGGGCGGCTTTCTCTCCCGGCGCCTTTCACAAAA‑CAAAATAAAATTCTGTGTCAAACCCCTCTTctttaaaac                                                                                                                                                                                                            >  2:675618/1‑230 (MQ=255)
                                                        tGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑GTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCCCAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCCGTTGATAAGGAACAGAAAGCCCAccgcccccc                                                                                                                                                                               >  1:186774/1‑236 (MQ=255)
                                                             ccGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCCCCCCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCCG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCCAAAGGCGGGCGATTTGCTCACGCGCGGCGTTTCCCAGAC‑CATAATAATGTTCGGTGTCCACCCCCCCTCCACCCTCTTGTTAAGGAAAAAAAAACACACCGCCCCCCGaaaa                                                                                                                                                                          >  2:1219285/1‑238 (MQ=255)
                                                                                                      tgACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAACCCACCGCCCCCAGAAATAGCCTGTAAATCCTGACCGAAATCGATGACCTGACGAAcc                                                                                                                                  >  2:1208694/1‑239 (MQ=255)
                                                                                                          gCTCCGACCCTGAGT‑ATTTAATACGCCAGACTGCGCTACCAGGAGGCGACCCGG‑TTCAATTT‑CCAGTTCCACAGGGTGGCCCAAAGGGCGGGCGTTTTGCTCACGCGCGGCATCCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  1:1219285/239‑1 (MQ=255)
                                                                                                          ccTCCGCACC‑GAGTAAATTAATCCGCAAGACTGCGCTCCAAGGAGGCGCCCCGG‑TTCAATTC‑CCAGTTCCACGGGGTGGCCAAAAGGGCGGGCAATTTGCTCACGCGCGGCATTCCACAGCC‑AATAAAAATGTCGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  1:1208694/238‑1 (MQ=255)
                                                                                                             ttgtttaaagggaaagtaatCCGCCGGACGGCGCTCCCAGGAAGCGCCCGGG‑TTCAATTT‑CCATTTGCACAGGGTGGCCAAAATGGGGGGCGTTTTGCCCACGCGCGGCATCCCACAGAC‑AAAAATAATGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  2:8517/221‑1 (MQ=255)
                                                                                                                gccttggcgttatttacaccccAGACTGCGCGACGAGGAGGCGACCCTT‑TTCAATTT‑CCGATGTCACAGGGTGCCCCAATGTGGGGGCGTTTTGCTCACTCGCGCTATTCCCAAGCC‑CATAATAGTTTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGCCCTGACGCACCATAGcacca                                                                                                                         <  2:337628/219‑1 (MQ=255)
                                                                                                                      aGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCt                                                                                                                  >  1:381135/1‑239 (MQ=255)
                                                                                                                      aGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATCCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGACAGCCCACCGCCCGCAGAAATAGCATGTAAATCCTGACCGAATTCGAGGACCTGACGCACCATAGcacctgaccctt                                                                                                                  >  1:1336188/1‑231 (MQ=255)
                                                                                                                        t‑aattaatACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGTTAAGGAACAGAAAGCCCCCCCCCCCCAGAAAAAGCCTGTAAACCCTGACCGAATTCGGTGACCtcaacacccaaaacccccaaca                                                                                                                     >  2:1120278/1‑213 (MQ=255)
                                                                                                                               agccGCAAGGCTGCGCAACCAGGAAGCCACCCGG‑TTAAATTC‑CCAGTTTCCCGGGGGGGCCAAAATGGCGGGCAATTTGCTCACGCGCGGAATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                          <  1:1008958/236‑1 (MQ=255)
                                                                                                                                 aCGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CCTAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTAGATGAACTGACGCCCCCTACCACCACCCACCCGTTCCaacga                                                                                                        >  2:1008958/1‑235 (MQ=255)
                                                                                                                                     caagctgccttACCAGGAAGGCACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCAAAAGGGCGGGCAATTTGCCCACGCGCGGCATTCCACAGAC‑AATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                                                                    <  1:1179127/229‑1 (MQ=255)
                                                                                                                                        aCTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACAccctgtcccttttacacat                                                                                                  >  1:8517/1‑221 (MQ=255)
                                                                                                                                            cgcTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCCGTTCCAGATGGGCATAATCa                                                                                              >  1:209642/1‑238 (MQ=255)
                                                                                                                                            cgcTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAAAAGCCTGTAAATCCTCACCGAATTCGATGACCCGACCCACCACAACCACACCACCCTGTTCCCGCTGGGCATAAtaa                                                                                              >  2:938049/1‑236 (MQ=255)
                                                                                                                                             gcTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACCCCTGTTCCAAATGGGCATAATCaa                                                                                             >  1:1270711/1‑238 (MQ=255)
                                                                                                                                                 ccgggagggcccccGT‑TTCAATTT‑CGAGTTTCCAGGGGGGGCCCAAAGGGGGGGGTGTTTGCGCCCGGGCGGATTTCAAAAGAC‑AAAAATAATGTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCa                                                                                        <  2:755384/228‑1 (MQ=255)
                                                                                                                                                 ccAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACCCGCGGCATTCCACCGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAAAAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGCCCGAATTCGCTGACATGACCGCCCAAAGCACCACCCACCCGTTCCAAGTGGGCCAAAACAACCCCa                                                                                        >  2:1307120/1‑239 (MQ=255)
                                                                                                                                                      aaGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTACATCCTGACCGAATTCGATGACCTTACGCACCATAGCACCACACCCCCGTCCCAACTGGGCTTAATCAACGCCAAAAcc                                                                                   >  1:128334/1‑239 (MQ=255)
                                                                                                                                                         gcgcccGG‑TTTTATTTCCCATTTTCAGGGGGGGCCCAAAATGGCGGGGAATTTCCCCCCGCGCGGCATCCCACAGAC‑CAAAATAATGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  2:381135/235‑1 (MQ=255)
                                                                                                                                                         gcccccGG‑TTAATTTT‑CCATTTCCCAAGGGGGCCCCAAAGGGGGGGGTTTTTGCCCACCCGCGGCATCCCACAACC‑AAAAAAAAGGTGCGGCATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  2:1270711/234‑1 (MQ=255)
                                                                                                                                                         agACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTTCACAGAC‑CATAATAATGTTCGGTATCCACCGCCTCTTTACCCTGTTGATAGGGAACAGAAAGCCCCCCGCCCGCAGAAATAGCCAGTAAATCCTGACCGACTTCGCTGCCCTTACCCACCATTGCACCACCCACCTTTTCCAACTGGGGTTATAAAACGCCCGAACCacc                                                                                >  2:1036991/2‑237 (MQ=255)
                                                                                                                                                          gACCCGGTTTCAATTT‑CCAGTTTAACAGGGTGGCCAAAATGGCGGGCGATTTGCTCACGCGGGGCATTCC‑CAGGCACATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  1:938049/237‑1 (MQ=255)
                                                                                                                                                           aCCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCCCACACCTGTTCCAGATGGGCATAAt                                                                                                >  1:16505/1‑221 (MQ=255)
                                                                                                                                                              cGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCCCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCCCACACCGGTTCCAGCTGGGCATAATAAACGCCAGAACCAAAGTGCa                                                                           >  1:552063/1‑239 (MQ=255)
                                                                                                                                                                        ttcccatgtccacaGGGTGGCCAAAAGTGCTGGCGATTTTCTCACGCCACGCATCCCACAGAC‑ACAAATGATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATg                                                                   <  1:672938/229‑1 (MQ=255)
                                                                                                                                                                                           gCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCacacac                                                                                                                    >  1:272719/1‑171 (MQ=255)
                                                                                                                                                                                                                  cACGCGCGGCATGCCACAGCC‑CAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTgcg                         <  1:305909/239‑1 (MQ=255)
                                                                                                                                                                                                                  cACGCGCGGCATGCCAAAGAC‑CATAATAAGTTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTgcg                         <  1:997651/239‑1 (MQ=255)
                                                                                                                                                                                                                            ctaccACAGAC‑ATAAAGAATGTTCGGTACAAACCGCCTCTTCACCCTGTTGAAAAGGAACAGAAAGCCCCCCGCCCGAAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  2:683485/236‑1 (MQ=255)
                                                                                                                                                                                                                            catCCACAGAC‑CATAAAAATGTCCGGTAACAACCGCCTCTTCACCCTGTTGAAAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  1:113186/237‑1 (MQ=255)
                                                                                                                                                                                                                             agaaaagggccaataaaaATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  1:470486/227‑1 (MQ=255)
                                                                                                                                                                                                                                cccaGAC‑AAAAATATTGTGCGGTATCACCCCCCTCTTCCCCCTTTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  2:16508/237‑1 (MQ=255)
                                                                                                                                                                                                                                ccaaGAC‑AAAAATAATGTTCGTTAACAACCGCCTTTTCACCCTGGTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:1052968/236‑1 (MQ=255)
                                                                                                                                                                                                                                agacagc‑cATAATAATGTTCGGTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:912065/233‑1 (MQ=255)
                                                                                                                                                                                                                                acaaGAC‑AAAAATAATTTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:750115/236‑1 (MQ=255)
                                                                                                                                                                                                                                     aC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCACTCTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGGGCGGCGGGAACATCGGGGTAcc      >  1:1354119/1‑239 (MQ=255)
                                                                                                                                                                                                                                     aC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGAACTAGTATCACGTCCAATGCGCCGGCCaaggccgcgtacc      >  1:442342/1‑227 (MQ=255)
                                                                                                                                                                                                                                     aC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACCACAGAACCAATGTGCATGGGAATGCCGACCAGCTTCAGATGATGACGTTGTATCACGTCCAGAGCGGCgggacaaaagggggacc      >  2:499725/1‑225 (MQ=255)
                                                                                                                                                                                                                                     aC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGGTCCAGATGGGCATAATCAACGCCAGAACCAATTGGCATGTGAATCCCGACCAGCTGCAGCGGATTACGTTGTATAACGCCCAGGGAGGGGGGCGGagagggggagc      >  2:96286/1‑229 (MQ=255)
                                                                                                                                                                                                                                        cATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCAAACACCTGTTCCCATTGGGCATAATCAACGCACGAACCAATTTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCCGATCGGGGTACCaa    >  2:1060460/1‑238 (MQ=255)
                                                                                                                                                                                                                                          taataaTGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCGGTGCGGCGGGAAGATCGGTGTACCaccg  >  1:683485/1‑236 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
TGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGT‑GGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAGAT  >  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: ‑