Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913233,1180TC42.1% 86.7 / 15.6 57I280V (ATT→GTT) mltDputative membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (31/2);  new base C (0/24);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.32e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTCAATGCCAGCATTTT  >  NC_000913/232969‑233259
                                                                                                                                                         |                                                                                                                                             
cGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑t                                                                                                                                               <  1:807184/150‑1 (MQ=255)
     gTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCg                                                                                                                                          <  2:754302/150‑1 (MQ=255)
                aTTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCt                                                                                                                               >  1:586365/1‑150 (MQ=255)
                      ccTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCg                                                                                                                        >  1:304206/1‑151 (MQ=255)
                      ccTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCg                                                                                                                        >  1:989568/1‑151 (MQ=255)
                           aGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCAtt                                                                                                                   >  2:417880/1‑151 (MQ=255)
                             ccaggcatccACGCATTTTATCTGCATGTTTTTTTGGCCCCCTCCCGCCGTGCGGCCCTCTGGCGCCGAGGGGGGACCT‑TTTT‑CCCCCAGGTTTGAT‑TTTTTTCAGTTTGCGGAGGAAAT‑TCCCCGCCATATCTGCTACCTTCGCCa                                                                                                                     <  1:711162/138‑1 (MQ=255)
                             caggggtttCCCGCGGTTGATCTGCTTCTTTTTTGGGCCCCTCCACGTATGGGGGGCCTCTTGCCCCCGGCGTGGCGCT‑TTTT‑CCGCCAGCTTTGAA‑TTTCTTCAGTTTGGTGAGGAAAA‑CCCCCGCCATATCTGCTACCTTCGcc                                                                                                                      <  1:494021/139‑1 (MQ=255)
                              cAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTc                                                                                                                 >  2:759892/1‑150 (MQ=255)
                                   aTTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCgg                                                                                                           >  2:886073/1‑151 (MQ=255)
                                      cccGCTGTTGTTCTGCTTTCTTTTTTGCCCCCCCCACGTGCGGCGGCCCCCTCGCCCCCGGGGTGGACCT‑TTTC‑CCGCCAGTGTGGAT‑TTTCTTCGGTTTGCGGAGGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGgct                                                                                                        <  2:90303/149‑1 (MQ=255)
                                      cACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGgct                                                                                                        >  2:1072479/1‑151 (MQ=255)
                                        cGCAGTTGCTCCGCTTGCTTCTTGGGCACCATCACGTTCGGCGGCCCTCTCGCCCCGGGGGTGG‑GCC‑TTTTCACGCACGCGGGAAT‑TTTCTACGGCTTGCTGCCGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGgctg                                                                                                       <  2:705852/150‑1 (MQ=255)
                                           aGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTc                                                                                                    >  1:406593/1‑150 (MQ=255)
                                           aGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTc                                                                                                    >  1:718051/1‑150 (MQ=255)
                                             ttGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAgg                                                                                                 >  1:570192/1‑151 (MQ=255)
                                               gATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTg                                                                                               >  1:970100/1‑151 (MQ=255)
                                                  cTGCTTGTTTTTTTGCCACCACCCCGTCCGGGGGCCCACTGCCGCCCGGCGGGGCGCT‑TTTT‑CCGCCAGGTTTTAT‑TTTTTTCAGCTTGCTGACGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGc                                                                                              <  2:923215/149‑1 (MQ=255)
                                                    gcttgtttcttGGGCACCACCGCGTCCGGCGGCCCCCGGCCCCCCGGCGGGGAGCT‑TTTT‑CCGCCGGGTTTGAA‑TTTCTTCAGCTTGCGAAGGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGca                                                                                             <  1:915506/142‑1 (MQ=255)
                                                    gCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCacac                                                                                          >  2:646011/1‑151 (MQ=255)
                                                          ttcttTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGcc                                                                                     >  1:885206/1‑150 (MQ=255)
                                                          ttcttTGCCACCACCACGTTCGGCGGGCCCTTCGCCCCCGGGGTGGAGTT‑TTTC‑ACCCGGGCGTTGAT‑TGTCTTCCGCTTGCTGGGGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACgcgc                                                                                      <  2:934805/149‑1 (MQ=255)
                                                           tcttTGCCCCCACCACGTACGGGGGCCCCCTGGCCCCCGGGGGGGAGCT‑TTTT‑CCGCCAGCGTTGAT‑TGTCTCCAGCTTGCTGACGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCa                                                                                    <  1:100720/150‑1 (MQ=255)
                                                             tttGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCag                                                                                   >  1:905281/1‑149 (MQ=255)
                                                                  cACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTAAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACgg                                                                            >  1:629682/1‑151 (MQ=255)
                                                                     cATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCtt                                                                         >  2:353748/1‑151 (MQ=255)
                                                                      caccgggtcctggggCCCCTTGCCCCCCGGCGGGGGCCTTTTT‑ACGCCGGCTTTGAT‑TTTCTTC‑GTTTTCTGGAGGAAACCCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCtt                                                                         <  1:496453/139‑1 (MQ=255)
                                                                                gCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGt                                                              >  1:712423/1‑151 (MQ=255)
                                                                                gCGGCCCCCTCGCCCCCAGGGGGGAGCT‑TTTT‑CCGCCGGGGTTGAA‑TGTCTCCAGCTTGCGGAAGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTc                                                                <  1:309089/149‑1 (MQ=255)
                                                                                gCGGCCCACGTGCCCCGGGGGGGGACCT‑TTTT‑CCGCCAGGGTTAAT‑TCTCTCCAGCTTGCTGAGGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCttt                                                                        <  2:716658/141‑1 (MQ=255)
                                                                                   gccctctgccgcCCACGGGGGGGCC‑TTTC‑CCGCCAGTGTTAAT‑TTTCTTCAGTTGGGGGACTGAAG‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTgg                                                           <  1:771513/142‑1 (MQ=255)
                                                                                     cccctgccccccGGCGGGGCGCT‑TTTT‑CCGCCGGCTTTGAT‑TGTCTTCAGTTTGCGGACGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTgg                                                           <  1:240917/140‑1 (MQ=255)
                                                                                     ccccgtGCCCCCGGGGTGGACCT‑TTTT‑CCGCCAGCTTTTTT‑TTTTTCCAGTTTGCGGGCGGAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCa                                                         <  2:777738/146‑1 (MQ=255)
                                                                                      cACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGc                                                          >  1:818650/1‑149 (MQ=255)
                                                                                         ttGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACg                                                     >  2:428503/1‑151 (MQ=255)
                                                                                            cgcCCAGCGTGGTGCC‑TTTT‑CCCCGACGGTTGAA‑TGTCCTCAGCTTCCTGACAGAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGt                                                    <  1:594640/149‑1 (MQ=255)
                                                                                              cccGGCGTGGCGCT‑TTTT‑CCCCCGGTGTTGAT‑TGTCTTCGGTTTGCTGACGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTAc                                                  <  1:753434/149‑1 (MQ=255)
                                                                                               ccAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGcc                                               >  2:1066141/1‑151 (MQ=255)
                                                                                                   cGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATa                                            >  1:659891/1‑150 (MQ=255)
                                                                                                     tGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACg                                         >  2:1075526/1‑151 (MQ=255)
                                                                                                             tttt‑CCGCCAGTGTTAAA‑TGTCTTCAGCTTGCTGACGAAAA‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGc                                    <  2:780527/149‑1 (MQ=255)
                                                                                                              ttttccgcagtcGTTTAA‑TTTCTTCAGTTTCCAGACGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGc                                    <  2:548146/138‑1 (MQ=255)
                                                                                                                    cgcaggtttGAATTGTTTTCGGTTTGCTGAGGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGtttt                              <  2:353483/143‑1 (MQ=255)
                                                                                                                      cAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTg                            >  2:301261/1‑149 (MQ=255)
                                                                                                                         cGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAat                       >  1:462765/1‑151 (MQ=255)
                                                                                                                         cGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAat                       >  2:669006/1‑151 (MQ=255)
                                                                                                                          gTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAat                       >  1:483976/1‑150 (MQ=255)
                                                                                                                             gAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATc                    >  1:898851/1‑150 (MQ=255)
                                                                                                                              aa‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCAc                  >  1:152855/1‑151 (MQ=255)
                                                                                                                                 tttttccggctggctgCCGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCa                   <  1:747765/140‑1 (MQ=255)
                                                                                                                                 tttctcCAGTTTGCTGACGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTgaga                         <  2:138150/140‑1 (MQ=255)
                                                                                                                                  gTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTc                >  2:663368/1‑150 (MQ=255)
                                                                                                                                   tCTTCATCTTGCTGACGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCa                   <  1:370433/146‑1 (MQ=255)
                                                                                                                                     tccAGCGTGCGGACGAAAT‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAatat                     <  2:712400/140‑1 (MQ=255)
                                                                                                                                              cTGACGAAAT‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTCAATGCCAGCAt     <  1:607093/149‑1 (MQ=255)
                                                                                                                                              cTGACGAAAC‑CCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTCAATGCCAGCa      <  1:267621/148‑1 (MQ=255)
                                                                                                                                                gACGGAAA‑TCCCCGCCATATCTGCTCCATTCGCCATTTCAACCGGGCTGCTCAGGTGAACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTCAATGCCAGCAtttt  >  2:990137/1‑150 (MQ=255)
                                                                                                                                                 aCGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTCAATGCCAGCAtttt  >  1:928992/1‑149 (MQ=255)
                                                                                                                                                         |                                                                                                                                             
CGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCC‑TTTT‑ACGCCAGCGTTGAA‑TGTCTTCAGCTTGCTGACGGAAA‑TCCCCGCCATATCTGCTACCTTCGCCATTTCAACCGGGCTGCTCAGGTGCACACGCGCCAGAGCACGGCTTTCATCGGTCGTTGGCAGACGTACGCCATAACGCTTGCTGTTTTTGAGAATATCACTCAATGCCAGCATTTT  >  NC_000913/232969‑233259

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