Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,329,743 G→A 100% V138M (GTG→ATG)  sohB → inner membrane protein, S49 peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,329,7430GA86.1% 205.8 / 19.6 85V138M (GTG→ATG) sohBinner membrane protein, S49 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base A (45/28);  total (51/34)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACT  >  NC_000913/1329646‑1329838
                                                                                                 |                                                                                               
gCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTtatg                                                                                               >  1:1134083/1‑100 (MQ=255)
gCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTtatg                                                                                               <  2:1513369/100‑1 (MQ=255)
   ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTc                                                                                            <  2:639010/100‑1 (MQ=255)
   ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTc                                                                                            <  2:413099/100‑1 (MQ=255)
   ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTc                                                                                            >  1:1290902/1‑100 (MQ=255)
   ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTc                                                                                            >  1:2501094/1‑100 (MQ=255)
   ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTc                                                                                            >  1:2407755/1‑100 (MQ=255)
      ttttAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGt                                                                                         <  2:865611/100‑1 (MQ=255)
      ttttAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGt                                                                                         >  1:179053/1‑100 (MQ=255)
       tttAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCCTTCAAACCGCAGGATCAGGTTATGCTCCGTc                                                                                        >  1:357279/1‑100 (MQ=255)
         tAAGGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTTAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCt                                                                                       <  1:45657/99‑1 (MQ=255)
            aGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGaa                                                                                   <  1:1361802/100‑1 (MQ=255)
            aGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGaa                                                                                   >  2:1018473/1‑100 (MQ=255)
            aGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGaa                                                                                   <  1:2407801/100‑1 (MQ=255)
             ggCAGCATGGACGCCCATGAAGTGCACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATCCAAACCGCAGGCTCAGGTTATGCTCCGTCTGGaaa                                                                                  >  2:402622/1‑100 (MQ=255)
               cagcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGc                                                                                >  1:402229/1‑100 (MQ=255)
                 gcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGccc                                                                              <  2:1426667/100‑1 (MQ=255)
                 gcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGccc                                                                              >  2:1274356/1‑100 (MQ=255)
                 gcaGGGAGGCCCATTAAGTGAACTCGCTACGTGAAGCGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGccc                                                                              <  1:1100268/100‑1 (MQ=255)
                   agggACGCCCATGAAGTGAGCTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCtg                                                                            <  1:877844/98‑1 (MQ=255)
                    tGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCtgg                                                                           >  2:1520923/1‑100 (MQ=255)
                    tGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCtgg                                                                           <  2:148500/100‑1 (MQ=255)
                      gACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCtggtg                                                                         >  2:395896/1‑100 (MQ=255)
                       aCGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCtggtgg                                                                        <  2:686975/100‑1 (MQ=255)
                       aCGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCtggtgg                                                                        <  2:1843219/100‑1 (MQ=255)
                          cccATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCAGATGgcat                                                                     >  2:212397/1‑100 (MQ=255)
                           ccATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatg                                                                    >  1:2285100/1‑100 (MQ=255)
                             aTGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggt                                                                  <  2:1164563/100‑1 (MQ=255)
                               gAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggtgc                                                                <  2:799276/100‑1 (MQ=255)
                               gAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggtgc                                                                <  2:1268588/100‑1 (MQ=255)
                                 aGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACTGCAGGATCAGGTTGTGCTCCGTCTAGAAAGCCCTGGTGgcatggtgcat                                                              >  1:662981/1‑100 (MQ=255)
                                 aGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcat                                                              >  1:1351117/1‑100 (MQ=255)
                                 aGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAAACGCAGGATCAGGTGATGCTCCGTCTGGAAAGCCCTGGTGgcatagtgcat                                                              >  1:799025/1‑100 (MQ=255)
                                  gTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatg                                                             >  2:1467683/1‑100 (MQ=255)
                                  gTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatg                                                             <  1:585626/100‑1 (MQ=255)
                                   tGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatgg                                                            >  2:1434836/1‑100 (MQ=255)
                                    gAACTCGCTACGTGAAGAGATAACGGCTGTATACGCAGCATTCAAACCGCAGGATCAGATTAGGCGCCGTATGGAAACCACTGGGGgcatggtgcatggt                                                           >  1:726642/1‑100 (MQ=255)
                                      aCTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTAGAAAGCCCTGGTGGCATGGTGCATGGTTa                                                         >  2:1021893/1‑100 (MQ=255)
                                       cTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTAc                                                        >  2:1822640/1‑100 (MQ=255)
                                        tCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACg                                                       <  1:56023/100‑1 (MQ=255)
                                          gCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACggg                                                     <  1:748786/100‑1 (MQ=255)
                                          gCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACggg                                                     <  2:1557197/100‑1 (MQ=255)
                                           cTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGc                                                    <  2:521920/100‑1 (MQ=255)
                                            tACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCGGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGc                                                    <  1:426883/99‑1 (MQ=255)
                                               gTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggc                                                >  1:303561/1‑100 (MQ=255)
                                               gTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggc                                                >  1:1202799/1‑100 (MQ=255)
                                                tGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcg                                               <  2:1809025/100‑1 (MQ=255)
                                                tGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcg                                               >  1:2052362/1‑100 (MQ=255)
                                                tGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCAGTCTGGAAAGCCCTGGTGGCATGGTGCATGGGTACTGGCTggcg                                               >  2:2255304/1‑100 (MQ=255)
                                                tGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGGTATGCTCCGTCTGGAAAGCCCTGGTGGAATGGTGCATGGTTACGGGCTggcg                                               >  1:589761/1‑100 (MQ=255)
                                                  aaGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcggc                                             >  2:2391318/1‑100 (MQ=255)
                                                   agagATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCt                                            >  1:1742718/1‑100 (MQ=255)
                                                    gagaTAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCtt                                           >  2:682781/1‑100 (MQ=255)
                                                     agaTAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTc                                          >  2:29176/1‑100 (MQ=255)
                                                        tAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCa                                       <  1:2392556/100‑1 (MQ=255)
                                                        tAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCa                                       <  2:402231/100‑1 (MQ=255)
                                                        tAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCa                                       >  1:2276330/1‑100 (MQ=255)
                                                           cGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACt                                    >  1:777906/1‑100 (MQ=255)
                                                            ggCTGTACGCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCATGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTg                                   >  1:1266913/1‑100 (MQ=255)
                                                                gTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGc                               >  1:1267313/1‑100 (MQ=255)
                                                                 tACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCg                              <  2:4553/100‑1 (MQ=255)
                                                                 tACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTAGCGGCTTCGCAACTGCAGc                               >  1:823906/1‑99 (MQ=255)
                                                                   cTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTc                            <  2:9782/100‑1 (MQ=255)
                                                                   cTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTc                            >  2:2500132/1‑100 (MQ=255)
                                                                   cTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTAGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTc                            >  1:831599/1‑100 (MQ=255)
                                                                     cGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGCAAGCCCTGGTGGAATGGTGCATGGTTACGGGCTGGCCGCTTCGCAACTGCAGCGTCt                           >  1:357161/1‑99 (MQ=255)
                                                                     cGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTg                          >  2:2426812/1‑100 (MQ=255)
                                                                     cGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTAGCATGGTGCTTGGTTACGGGCTGGGGGCTTCGCAACTGCAGCGTCTg                          >  1:1750854/1‑100 (MQ=255)
                                                                      gcagcaTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTgc                         >  2:1320928/1‑100 (MQ=255)
                                                                       cagcaTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTgcg                        <  1:1274351/100‑1 (MQ=255)
                                                                        agcaTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTgcgc                       >  1:30842/1‑100 (MQ=255)
                                                                        agcaTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTgcgc                       <  1:1493547/100‑1 (MQ=255)
                                                                           aTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGTCTTCGCAACTGCAGCGTCTGCGCGAt                    <  2:1579091/100‑1 (MQ=255)
                                                                            ttCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATa                   >  2:1666113/1‑100 (MQ=255)
                                                                                 aCCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAAc              >  1:1437198/1‑100 (MQ=255)
                                                                                 aCCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAAc              >  2:1957821/1‑100 (MQ=255)
                                                                                  ccGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACa             >  2:35940/1‑100 (MQ=255)
                                                                                    gcgggATCAGGTTGTGCTCCGTGGGGAACGCCCCGGTGGCATGGTGCATGGTTACGGCCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAAc              <  1:576133/94‑1 (MQ=255)
                                                                                      aGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTcc         <  2:272238/100‑1 (MQ=255)
                                                                                      aGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTcc         <  2:1316814/100‑1 (MQ=255)
                                                                                      aGGATCAGGTTATGCACCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTcc         <  2:2549663/100‑1 (MQ=255)
                                                                                         aTCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGTGCGATAAAAACATTCCttt      >  2:806280/1‑100 (MQ=255)
                                                                                         aTCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCttt      >  1:1736818/1‑100 (MQ=255)
                                                                                          tCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTa     <  2:1893547/100‑1 (MQ=255)
                                                                                            aGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAAc   >  2:2580657/1‑100 (MQ=255)
                                                                                            aGGTTATGCTCCGTCTGGAAAGCACTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAAc   >  1:1176972/1‑100 (MQ=255)
                                                                                             ggTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCCTCTCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACt  >  1:1818864/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACT  >  NC_000913/1329646‑1329838

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: