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

ATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTG  >  NC_000913/2999576‑2999766
                                                                                                   |                                                                                               
aTGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑gg                                                                                             >  2:2220656/1‑100 (MQ=255)
 tGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑ggg                                                                                            >  1:454667/1‑100 (MQ=255)
     gTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑ggggggg                                                                                        <  2:133423/100‑1 (MQ=255)
         tGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGgaaa                                                                                    <  1:2548111/100‑2 (MQ=255)
          ggTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGAGTAA‑‑GGGGGGGgaaat                                                                                   >  1:214235/1‑98 (MQ=255)
           gTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATg                                                                                   >  2:667693/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGcc                                                                                 >  2:1043632/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGcc                                                                                 >  2:130073/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGcc                                                                                 >  2:1353372/1‑100 (MQ=255)
             cGTTGATTTGGCGGTTGATTTACGTGTAAGACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGcc                                                                                 <  1:2412663/100‑1 (MQ=255)
              gTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCg                                                                                <  1:308549/100‑1 (MQ=255)
              gTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCg                                                                                <  1:2143135/100‑1 (MQ=255)
               ttGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGt                                                                               >  2:1383949/1‑100 (MQ=255)
               ttGATTTGGCGGTTGATTTACGTGTAATACGGCTACTTTTCGCCCCACCCAGTTTTACTTTCTGACCAACTTTAATGGTGTAAGGGGGGGGAAATGCCGt                                                                               <  2:1929773/100‑1 (MQ=255)
                 gATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTc                                                                             <  1:718144/100‑1 (MQ=255)
                  aTTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCa                                                                            <  2:711849/100‑1 (MQ=255)
                    ttGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGt                                                                          >  2:688792/1‑100 (MQ=255)
                      ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGGAAATGCCGTTCAGTc                                                                         <  1:1063961/100‑1 (MQ=255)
                      ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcg                                                                        >  2:327327/1‑100 (MQ=255)
                      ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcg                                                                        >  2:327424/1‑100 (MQ=255)
                        cGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcgcg                                                                      <  2:2127867/100‑1 (MQ=255)
                          gTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGcc                                                                    <  2:1705984/100‑1 (MQ=255)
                            tGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAg                                                                  <  1:1859441/100‑1 (MQ=255)
                            tGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTGTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAg                                                                  <  1:2377526/100‑1 (MQ=255)
                             gATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGt                                                                 >  2:2287353/1‑100 (MQ=255)
                               tttACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTATAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTc                                                               >  2:2334094/1‑100 (MQ=255)
                                ttACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCt                                                              <  2:1097112/100‑1 (MQ=255)
                                 tACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCtt                                                             >  2:1284828/1‑100 (MQ=255)
                                  aCGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCttt                                                            <  1:1387081/100‑1 (MQ=255)
                                    gtgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                          >  2:96393/1‑100 (MQ=255)
                                    gtgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                          <  1:986884/100‑1 (MQ=255)
                                    gtgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                          >  1:2020073/1‑100 (MQ=255)
                                     tgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAc                                                         >  2:886611/1‑100 (MQ=255)
                                     tgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAc                                                         <  1:2644590/100‑1 (MQ=255)
                                     tgtAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGCTCAGTCGCGCCAGTTCTTTTAc                                                         >  1:834880/1‑100 (MQ=255)
                                      gtAATACTGCTCCTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGCAATGCCGTTCAGTCGCGCCAGTTCTTTTACg                                                        >  1:1644319/1‑100 (MQ=255)
                                      gtAATACTGCTACTTTTCGCCCCACCCCGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACg                                                        <  1:2242151/100‑1 (MQ=255)
                                       tAATATTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGc                                                       <  1:595812/100‑1 (MQ=255)
                                          tACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTgg                                                    <  1:1039305/100‑1 (MQ=255)
                                          tACTGCGACTTTGCTCCCCACGCAGTTTTCGTTTCTGGCGAACTTCCATCGTGTAATGGGCGGGGAGTGGCTTTAAGTCGCGCCAGTTCTTTTACGCggg                                                    <  2:2236439/100‑1 (MQ=255)
                                           aCTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGt                                                   >  2:188550/1‑100 (MQ=255)
                                             tGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTc                                                 >  2:1420262/1‑100 (MQ=255)
                                               cTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTccc                                               <  1:1414515/100‑1 (MQ=255)
                                               cTACTTTTCGCCCCACCCAGTTTTAGTTTCCGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTccc                                               <  1:289100/100‑1 (MQ=255)
                                                tACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTGAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                              >  2:1994227/1‑100 (MQ=255)
                                                tACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                              >  1:553951/1‑100 (MQ=255)
                                                tACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GTGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                              <  2:2006808/100‑1 (MQ=255)
                                                 aCTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgt                                             <  1:1667320/100‑1 (MQ=255)
                                                   ttttCGCCCCACCCAGTTTTCGTCTCAGACCAACTTCAAAGGTGTCAGGGGGGGGAAATGCCGTTCAGTCGCGCCTGTTCTTTTACGCCGGTTCCCgtgg                                           <  1:2469934/100‑1 (MQ=255)
                                                   ttttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtgg                                           <  1:2448355/100‑1 (MQ=255)
                                                   ttttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtgg                                           <  2:2042704/100‑1 (MQ=255)
                                                    tttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtggt                                          >  2:1396999/1‑100 (MQ=255)
                                                    tttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCCCCAGTTCTTTTACGCTGGTTCCCgtggt                                          >  2:2121644/1‑100 (MQ=255)
                                                        gCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTgcgc                                      >  1:358990/1‑100 (MQ=255)
                                                          cccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGa                                    <  1:199287/100‑1 (MQ=255)
                                                           ccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGaa                                   <  1:1546924/100‑1 (MQ=255)
                                                            caccaaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCCGTTCTTTTCCGCTGGTTCCCCTGGTGGGCGaaa                                  >  2:2205496/1‑100 (MQ=255)
                                                                cggTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACg                              <  2:2366171/98‑1 (MQ=255)
                                                                caGTTTTAGTTTCTGACAAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACg                              <  2:1456092/100‑1 (MQ=255)
                                                                 aGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGa                             <  2:2182077/100‑1 (MQ=255)
                                                                 aGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGat                            <  2:329286/100‑1 (MQ=255)
                                                                  gTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGat                            >  1:1094317/1‑100 (MQ=255)
                                                                   ttttAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGata                           >  2:1637273/1‑100 (MQ=255)
                                                                      tAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAg                        >  1:598992/1‑100 (MQ=255)
                                                                        gTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                      >  2:238940/1‑100 (MQ=255)
                                                                        gTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                      <  1:1318746/100‑1 (MQ=255)
                                                                        gTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                      <  1:267792/100‑1 (MQ=255)
                                                                         tttCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGt                     >  2:328186/1‑100 (MQ=255)
                                                                               aCCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccc               >  2:1548686/1‑100 (MQ=255)
                                                                                  aaCTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCg            >  2:219803/1‑100 (MQ=255)
                                                                                    cTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGtt          >  2:2211574/1‑100 (MQ=255)
                                                                                       cAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCa       >  2:1372977/1‑100 (MQ=255)
                                                                                         aTGGTGTATT‑‑GGGGGGGATATGCCGTTCATTTGTTTTAGTTCTTTTACGCTTGTTCGCGTTGTGCGCTAAGTACGATATAGCGTATCCCGCCGTTTCTCg     >  2:119263/1‑100 (MQ=255)
                                                                                         aTGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACg     <  1:2331957/100‑1 (MQ=255)
                                                                                         aTGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACg     <  1:46790/100‑1 (MQ=255)
                                                                                          ttgtgtAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACgg    <  1:1579473/98‑1 (MQ=255)
                                                                                          tGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACgg    <  1:1571202/100‑1 (MQ=255)
                                                                                          tGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACgg    <  1:433488/100‑1 (MQ=255)
                                                                                           ggTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgt   <  1:17344/100‑1 (MQ=255)
                                                                                           ggTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgt   <  1:2253283/100‑1 (MQ=255)
                                                                                            gtgtAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtg  >  2:1088578/1‑100 (MQ=255)
                                                                                            gtgtAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtg  <  1:2448812/100‑1 (MQ=255)
                                                                                            gtgtAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtg  >  2:1907260/1‑100 (MQ=255)
                                                                                                   |                                                                                               
ATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTG  >  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: