Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,999,680 (G)7→8 coding (212/756 nt) ygeR ← LysM domain‑containing M23 family putative peptidase; septation lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,999,6731.G93.7% 260.9 / 12.4 80coding (219/756 nt)ygeRLysM domain‑containing M23 family putative peptidase; septation lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (35/39);  total (37/43)
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.73e-01

ATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTG  >  NC_000913/2999576‑2999766
                                                                                                   |                                                                                               
aTGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑g‑g                                                                                             >  2:2220656/1‑100 (MQ=255)
 tGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑g‑gg                                                                                            >  1:454667/1‑100 (MQ=255)
     gTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑g‑gggggg                                                                                        <  2:133423/100‑1 (MQ=255)
         tGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑G‑GGGGGGgaaa                                                                                    <  1:2548111/100‑2 (MQ=255)
          ggTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGAGTA‑A‑‑G‑GGGGGGgaaat                                                                                   >  1:214235/1‑98 (MQ=255)
           gTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATg                                                                                   >  2:667693/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGcc                                                                                 >  2:1043632/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGcc                                                                                 >  2:130073/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGcc                                                                                 >  2:1353372/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAAGACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGcc                                                                                 <  1:2412663/100‑1 (MQ=255)
              gTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCg                                                                                <  1:308549/100‑1 (MQ=255)
              gTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCg                                                                                <  1:2143135/100‑1 (MQ=255)
               ttGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGt                                                                               >  2:1383949/1‑100 (MQ=255)
               ttGATTTGGCGGTTGATTTACGTGTAATACGGCTACTTTTCGCCCCACCCAGTTTTACTTTCTGACCAACTTTAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGt                                                                               <  2:1929773/100‑1 (MQ=255)
                 gATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTc                                                                             <  1:718144/100‑1 (MQ=255)
                  aTTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCa                                                                            <  2:711849/100‑1 (MQ=255)
                    ttGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGt                                                                          >  2:688792/1‑100 (MQ=255)
                      ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AGGG‑GGGGGGAAATGCCGTTCAGTc                                                                         <  1:1063961/100‑1 (MQ=255)
                      ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTcg                                                                        >  2:327327/1‑100 (MQ=255)
                      ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTcg                                                                        >  2:327424/1‑100 (MQ=255)
                        cGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTcgcg                                                                      <  2:2127867/100‑1 (MQ=255)
                          gTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGcc                                                                    <  2:1705984/100‑1 (MQ=255)
                            tGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAg                                                                  <  1:1859441/100‑1 (MQ=255)
                            tGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTGTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAg                                                                  <  1:2377526/100‑1 (MQ=255)
                             gATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGt                                                                 >  2:2287353/1‑100 (MQ=255)
                               tttACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTATA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTc                                                               >  2:2334094/1‑100 (MQ=255)
                                ttACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCt                                                              <  2:1097112/100‑1 (MQ=255)
                                 tACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCtt                                                             >  2:1284828/1‑100 (MQ=255)
                                  aCGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCttt                                                            <  1:1387081/100‑1 (MQ=255)
                                    gtgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                          >  2:96393/1‑100 (MQ=255)
                                    gtgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                          <  1:986884/100‑1 (MQ=255)
                                    gtgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                          >  1:2020073/1‑100 (MQ=255)
                                     tgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAc                                                         >  2:886611/1‑100 (MQ=255)
                                     tgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAc                                                         <  1:2644590/100‑1 (MQ=255)
                                     tgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGCTCAGTCGCGCCAGTTCTTTTAc                                                         >  1:834880/1‑100 (MQ=255)
                                      gtAATACTGCTCCTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGCAATGCCGTTCAGTCGCGCCAGTTCTTTTACg                                                        >  1:1644319/1‑100 (MQ=255)
                                      gtAATACTGCTACTTTTCGCCCCACCCCGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACg                                                        <  1:2242151/100‑1 (MQ=255)
                                       tAATATTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGc                                                       <  1:595812/100‑1 (MQ=255)
                                          tACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTgg                                                    <  1:1039305/100‑1 (MQ=255)
                                          tACTGCGACTTTGCTCCCCACGCAGTTTTCGTTTCTGGCGAACTTCCATCGTGTA‑AT‑G‑GGCGGGGAGTGGCTTTAAGTCGCGCCAGTTCTTTTACGCggg                                                    <  2:2236439/100‑1 (MQ=255)
                                           aCTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGt                                                   >  2:188550/1‑100 (MQ=255)
                                             tGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTc                                                 >  2:1420262/1‑100 (MQ=255)
                                               cTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTccc                                               <  1:1414515/100‑1 (MQ=255)
                                               cTACTTTTCGCCCCACCCAGTTTTAGTTTCCGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTccc                                               <  1:289100/100‑1 (MQ=255)
                                                tACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTGAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                              >  2:1994227/1‑100 (MQ=255)
                                                tACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                              >  1:553951/1‑100 (MQ=255)
                                                tACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑GTGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                              <  2:2006808/100‑1 (MQ=255)
                                                 aCTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgt                                             <  1:1667320/100‑1 (MQ=255)
                                                   ttttCGCCCCACCCAGTTTTCGTCTCAGACCAACTTCAAAGGTGTC‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCTGTTCTTTTACGCCGGTTCCCgtgg                                           <  1:2469934/100‑1 (MQ=255)
                                                   ttttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtgg                                           <  1:2448355/100‑1 (MQ=255)
                                                   ttttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtgg                                           <  2:2042704/100‑1 (MQ=255)
                                                    tttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtggt                                          >  2:1396999/1‑100 (MQ=255)
                                                    tttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCCCCAGTTCTTTTACGCTGGTTCCCgtggt                                          >  2:2121644/1‑100 (MQ=255)
                                                        gCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTgcgc                                      >  1:358990/1‑100 (MQ=255)
                                                          cccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGa                                    <  1:199287/100‑1 (MQ=255)
                                                           ccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGaa                                   <  1:1546924/100‑1 (MQ=255)
                                                            caccaaGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCCGTTCTTTTCCGCTGGTTCCCCTGGTGGGCGaaa                                  >  2:2205496/1‑100 (MQ=255)
                                                                cggTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACg                              <  2:2366171/98‑1 (MQ=255)
                                                                caGTTTTAGTTTCTGACAAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACg                              <  2:1456092/100‑1 (MQ=255)
                                                                 aGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGa                             <  2:2182077/100‑1 (MQ=255)
                                                                 aGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGat                            <  2:329286/100‑1 (MQ=255)
                                                                  gTTTTAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGat                            >  1:1094317/1‑100 (MQ=255)
                                                                   ttttAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGata                           >  2:1637273/1‑100 (MQ=255)
                                                                      tAGTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAg                        >  1:598992/1‑100 (MQ=255)
                                                                        gTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                      >  2:238940/1‑100 (MQ=255)
                                                                        gTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                      <  1:1318746/100‑1 (MQ=255)
                                                                        gTTTCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                      <  1:267792/100‑1 (MQ=255)
                                                                         tttCTGACCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGt                     >  2:328186/1‑100 (MQ=255)
                                                                               aCCAACTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccc               >  2:1548686/1‑100 (MQ=255)
                                                                                  aaCTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCg            >  2:219803/1‑100 (MQ=255)
                                                                                    cTTCAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGtt          >  2:2211574/1‑100 (MQ=255)
                                                                                       cAATGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCa       >  2:1372977/1‑100 (MQ=255)
                                                                                         aTGGTGTATT‑‑G‑GGGGGGATATGCCGTTCATTTGTTTTAGTTCTTTTACGCTTGTTCGCGTTGTGCGCTAAGTACGATATAGCGTATCCCGCCGTTTCTCg     >  2:119263/1‑100 (MQ=255)
                                                                                         aTGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACg     <  1:2331957/100‑1 (MQ=255)
                                                                                         aTGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACg     <  1:46790/100‑1 (MQ=255)
                                                                                          ttgtgtA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACgg    <  1:1579473/98‑1 (MQ=255)
                                                                                          tGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACgg    <  1:1571202/100‑1 (MQ=255)
                                                                                          tGGTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACgg    <  1:433488/100‑1 (MQ=255)
                                                                                           ggTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgt   <  1:17344/100‑1 (MQ=255)
                                                                                           ggTGTA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgt   <  1:2253283/100‑1 (MQ=255)
                                                                                            gtgtA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtg  >  2:1088578/1‑100 (MQ=255)
                                                                                            gtgtA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtg  <  1:2448812/100‑1 (MQ=255)
                                                                                            gtgtA‑AG‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtg  >  2:1907260/1‑100 (MQ=255)
                                                                                                   |                                                                                               
ATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTA‑A‑‑G‑GGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTG  >  NC_000913/2999576‑2999766

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