Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,329,743 G→A 29.6% 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,7430GA29.6% 206.0 / 100.3 135V138M (GTG→ATG) sohBinner membrane protein, S49 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (53/42);  new base A (25/15);  total (78/57)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

GTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGT  >  NC_000913/1329644‑1329840
                                                                                                   |                                                                                                 
gTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTtg                                                                                                   <  1:1053316/100‑1 (MQ=255)
 tGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTtgt                                                                                                  <  2:2120875/100‑1 (MQ=255)
   cTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGc                                                                                                >  1:266400/1‑100 (MQ=255)
    tGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGCTCCGGTTGTGc                                                                                                >  1:1345939/1‑99 (MQ=255)
    tGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCt                                                                                               >  2:4405737/1‑100 (MQ=255)
     ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTc                                                                                              >  2:3716387/1‑100 (MQ=255)
     ggATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTc                                                                                              >  2:1429207/1‑100 (MQ=255)
     ggATTTGAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTc                                                                                              <  1:2642102/100‑1 (MQ=255)
      gATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTcc                                                                                             <  2:2083236/100‑1 (MQ=255)
      gATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTcc                                                                                             >  1:1531770/1‑100 (MQ=255)
      gATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTcc                                                                                             >  2:2897113/1‑100 (MQ=255)
           tAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTg                                                                                        <  2:3294878/100‑1 (MQ=255)
             aaGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGTATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGa                                                                                      <  1:1568542/100‑1 (MQ=255)
              aGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGaa                                                                                     >  1:1288023/1‑100 (MQ=255)
              aGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGaa                                                                                     >  2:3421686/1‑100 (MQ=255)
              aGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGa                                                                                      >  2:682975/1‑99 (MQ=255)
               ggCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGaaa                                                                                    >  2:1870842/1‑100 (MQ=255)
                gcagcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAg                                                                                   >  2:4008720/1‑100 (MQ=255)
                 cagcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGc                                                                                  >  1:956740/1‑100 (MQ=255)
                 cagcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGc                                                                                  >  2:2276350/1‑100 (MQ=255)
                 cagcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGc                                                                                  >  2:1098115/1‑100 (MQ=255)
                   gcaTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGccc                                                                                <  2:2998869/100‑1 (MQ=255)
                    caTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCTCAGCATTCAAACCTCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCt                                                                               >  1:596494/1‑100 (MQ=255)
                    caTGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCt                                                                               <  2:1126798/100‑1 (MQ=255)
                          cGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgc                                                                         >  1:2479834/1‑100 (MQ=255)
                           gCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgca                                                                        <  2:1087162/100‑1 (MQ=255)
                              cATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatgg                                                                     <  1:3332964/100‑1 (MQ=255)
                                tGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtg                                                                   <  2:231641/100‑1 (MQ=255)
                                 gAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgc                                                                  >  2:2399747/1‑100 (MQ=255)
                                  aaGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAAACGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgca                                                                 >  1:4150146/1‑100 (MQ=255)
                                    gTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatg                                                               <  2:4230222/100‑1 (MQ=255)
                                    gTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatg                                                               <  1:3421665/100‑1 (MQ=255)
                                     tGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatgg                                                              <  1:3124079/100‑1 (MQ=255)
                                      gAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatggt                                                             >  1:2294456/1‑100 (MQ=255)
                                      gAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatggt                                                             >  2:2425537/1‑100 (MQ=255)
                                      gAACTCGCTACGCGAAGAGATAACGGCTGTACTCGCAGCCTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGgcatggtgcatggt                                                             >  1:2020615/1‑100 (MQ=255)
                                        aCTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTa                                                           <  1:3565388/100‑1 (MQ=255)
                                        aCTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTa                                                           <  2:3136452/100‑1 (MQ=255)
                                           cGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACgg                                                        >  2:351519/1‑100 (MQ=255)
                                           cGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACgg                                                        >  2:2553581/1‑100 (MQ=255)
                                           cGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACgg                                                        <  2:2947649/100‑1 (MQ=255)
                                           cGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACgg                                                        <  2:3131738/100‑1 (MQ=255)
                                           cGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACg                                                         <  2:4045863/99‑1 (MQ=255)
                                             cTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGc                                                      >  2:2495849/1‑100 (MQ=255)
                                             cTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGc                                                      <  2:3043321/100‑1 (MQ=255)
                                             cTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGc                                                      <  1:2615100/100‑1 (MQ=255)
                                               aCGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTg                                                    <  1:2534415/100‑1 (MQ=255)
                                                cGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTgg                                                   <  1:1429201/100‑1 (MQ=255)
                                                cGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTgg                                                   >  2:2826867/1‑100 (MQ=255)
                                                 gTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggc                                                  >  2:1397783/1‑100 (MQ=255)
                                                   gAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcgg                                                >  1:1564296/1‑100 (MQ=255)
                                                   gAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcgg                                                >  2:3676862/1‑100 (MQ=255)
                                                   gAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTggcgg                                                >  1:105547/1‑100 (MQ=255)
                                                      gagaTAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCtt                                             >  2:1765959/1‑100 (MQ=255)
                                                        gaTAACGGCTGTCCTCGCAGCATTCAAACCGCAGTCTCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCg                                           <  2:3699566/100‑1 (MQ=255)
                                                        gaTAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCg                                           >  2:3402961/1‑100 (MQ=255)
                                                        gaTAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCg                                           >  1:479593/1‑100 (MQ=255)
                                                         aTAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGc                                          <  2:4139305/100‑1 (MQ=255)
                                                          tAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCa                                         >  2:3228007/1‑100 (MQ=255)
                                                          tAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCa                                         >  2:1539128/1‑100 (MQ=255)
                                                             cGGCTGTACTCGCAGCATTCGAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACt                                      <  1:3196803/100‑1 (MQ=255)
                                                               gCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGc                                    >  1:1929327/1‑100 (MQ=255)
                                                               gCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGc                                    <  1:2849733/100‑1 (MQ=255)
                                                                cTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCa                                   <  2:1504547/100‑1 (MQ=255)
                                                                   tACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCg                                >  2:1720284/1‑100 (MQ=255)
                                                                     cTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTc                              <  2:4328587/100‑1 (MQ=255)
                                                                     cTCGCAGCATTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTc                              >  2:3946416/1‑100 (MQ=255)
                                                                       cGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTg                            >  2:1214229/1‑100 (MQ=255)
                                                                       cGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTg                            >  2:634521/1‑100 (MQ=255)
                                                                          agcaTTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTgcgc                         <  2:3484469/100‑1 (MQ=255)
                                                                           gcaTTCAAACCGTAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGcgcg                        >  2:2898586/1‑100 (MQ=255)
                                                                           gcaTTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTGTGcgcg                        >  1:1960513/1‑100 (MQ=255)
                                                                           gcaTTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGcgcg                        >  1:1687288/1‑100 (MQ=255)
                                                                            caTTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGa                       >  1:3743545/1‑100 (MQ=255)
                                                                             aTTCAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGAt                      <  2:4378523/100‑1 (MQ=255)
                                                                              ttCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATa                     >  1:628547/1‑100 (MQ=255)
                                                                               tCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATaa                    >  1:3300221/1‑100 (MQ=255)
                                                                               tCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATaa                    >  1:235612/1‑100 (MQ=255)
                                                                                cAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCCACTGCAGCGTCTGCGCGATaaa                   >  1:699436/1‑100 (MQ=255)
                                                                                cAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATaaa                   >  1:3116625/1‑100 (MQ=255)
                                                                                cAAACCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCACTGGGGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGcgcg                        >  1:2089168/1‑95 (MQ=255)
                                                                                 aaaCCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATaaaa                  >  2:3421322/1‑100 (MQ=255)
                                                                                  aaCCGCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATaaaaa                 >  1:4012484/1‑100 (MQ=255)
                                                                                   aCCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCACTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAAc                <  1:2313351/100‑1 (MQ=255)
                                                                                      gCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACAtt             <  1:589130/100‑1 (MQ=255)
                                                                                      gCAGGATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAAc                <  1:1179980/97‑1 (MQ=255)
                                                                                         ggATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCt          >  2:546133/1‑100 (MQ=255)
                                                                                         ggATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCt          >  1:3546836/1‑100 (MQ=255)
                                                                                         ggATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCt          >  2:2656166/1‑100 (MQ=255)
                                                                                         ggATCAGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTcc           <  2:863152/99‑1 (MQ=255)
                                                                                          gATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCtt         >  2:431963/1‑100 (MQ=255)
                                                                                           aTCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCttt        <  1:1126061/100‑1 (MQ=255)
                                                                                            tCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTa       <  2:779498/100‑1 (MQ=255)
                                                                                             cAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTaa      <  2:941749/100‑1 (MQ=255)
                                                                                              aGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAAc     <  2:2233083/100‑1 (MQ=255)
                                                                                              aGGTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAAc     >  1:2830582/1‑100 (MQ=255)
                                                                                               ggTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACt    >  2:1036823/1‑100 (MQ=255)
                                                                                                gTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTg   >  1:4274551/1‑100 (MQ=255)
                                                                                                gTTATGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTg   >  1:4189453/1‑100 (MQ=255)
                                                                                                 ttGTGCGCCGTCTGGAAAGCCCTGGTGGCAGGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGt  <  1:2075056/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GTGCTGGATTTTAAAGGCAGCATGGACGCCCATGAAGTGAACTCGCTACGTGAAGAGATAACGGCTGTACTCGCAGCATTCAAACCGCAGGATCAGGTTGTGCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGT  >  NC_000913/1329644‑1329840

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