Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,388,496 T→A 5.2% Q292L (CAG→CTG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,4960TA 5.2% 310.4 / 3.9 93Q292L (CAG→CTG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base A (3/2);  total (48/45)
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.81e-01

GCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTA  >  NC_000913/4388349‑4388642
                                                                                                                                                   |                                                                                                                                                  
gcCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTctg                                                                                                                                                   <  6:325233/149‑1 (MQ=255)
   aTTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTctggct                                                                                                                                                >  2:445361/1‑149 (MQ=255)
        ggTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCa                                                                                                                                            <  5:364378/148‑1 (MQ=255)
         gTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGc                                                                                                                                          >  1:334045/1‑149 (MQ=255)
         gTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGACACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGc                                                                                                                                          <  7:57259/149‑1 (MQ=255)
         ctATCAAGCTGTTGTGGTTTCGGCATTTCCGTCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTTGCAGc                                                                                                                                          <  6:13827/148‑1 (MQ=255)
           aTCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCt                                                                                                                                        >  4:192251/1‑149 (MQ=255)
            tCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCt                                                                                                                                        <  2:340499/148‑1 (MQ=255)
              acgCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGAc                                                                                                                                     <  1:224553/147‑1 (MQ=255)
              aaGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGATCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                      <  7:350599/148‑1 (MQ=255)
               aGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACg                                                                                                                                    >  5:81662/1‑149 (MQ=255)
                    ttgtGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGtt                                                                                                                               <  1:204008/149‑1 (MQ=255)
                      gtgGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTg                                                                                                                              <  5:32194/148‑1 (MQ=255)
                      gtgGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                                                                             >  5:414263/1‑149 (MQ=255)
                        ggTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGa                                                                                                                           <  1:445361/149‑1 (MQ=255)
                         gTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAg                                                                                                                          >  7:433718/1‑149 (MQ=255)
                         gTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAg                                                                                                                          <  2:249679/149‑1 (MQ=255)
                          tttCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAgg                                                                                                                         >  4:454922/1‑149 (MQ=255)
                              ggCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                                                                                      <  1:69294/148‑1 (MQ=255)
                              ggCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                                                                                     >  8:169486/1‑149 (MQ=255)
                               gCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGACgt                                                                                                                                   >  2:96240/1‑133 (MQ=255)
                               gCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGACgt                                                                                                                                   <  1:96240/134‑2 (MQ=255)
                                cATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                                                                                     >  7:58048/1‑147 (MQ=255)
                                 aTTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGGGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                                                                                     >  1:155156/1‑146 (MQ=255)
                                            ccGGGCCCCCTGCGCCCGCAGCGCGTCGCCGAGCAGATTGGACGCGCCCCGCCATTGCGACTGTTCACGCAGCGTCTGCCACGCCTGCCGGACCTGTAACGTCTTGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                                                                                       <  7:21238/149‑1 (MQ=255)
                                              gTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                                                                                     >  3:206456/1‑149 (MQ=255)
                                                 ccACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCg                                                                                                         >  3:86695/1‑142 (MQ=255)
                                                 ccACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCg                                                                                                         <  4:86695/142‑1 (MQ=255)
                                                 ccACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGcc                                                                                                  >  1:116225/1‑149 (MQ=255)
                                                   aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                                                                                               >  8:232760/1‑134 (MQ=255)
                                                   aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                                                                                               <  7:232760/134‑1 (MQ=255)
                                                   aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTg                                                                                                >  6:388050/1‑149 (MQ=255)
                                                   aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTg                                                                                                >  5:109157/1‑149 (MQ=255)
                                                    ccTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGt                                                                                               <  5:367205/149‑1 (MQ=255)
                                                    ccTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGt                                                                                               <  6:81662/149‑1 (MQ=255)
                                                     cTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGtt                                                                                              <  5:388050/149‑1 (MQ=255)
                                                        cgcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGa                                                                                            <  8:433718/148‑1 (MQ=255)
                                                        cgcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAt                                                                                           <  8:318716/149‑1 (MQ=255)
                                                         gcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAtt                                                                                          <  2:152758/149‑1 (MQ=255)
                                                           ccGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTc                                                                                         >  2:379415/1‑148 (MQ=255)
                                                             gCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa                                                                                      <  4:206212/149‑1 (MQ=255)
                                                                  gcTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                                                                                 <  2:116225/149‑1 (MQ=255)
                                                                                aTTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGt                                                                   <  1:39600/149‑1 (MQ=255)
                                                                                    gACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGtt                                                                                                                               <  7:332815/85‑1 (MQ=255)
                                                                                    gACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGtt                                                                                                                               >  8:332815/1‑85 (MQ=255)
                                                                                     aCGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAtt                                                              <  6:414263/149‑1 (MQ=255)
                                                                                           ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGaa                                                         <  2:336146/148‑1 (MQ=255)
                                                                                           ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGaa                                                         <  2:334045/148‑1 (MQ=255)
                                                                                           ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAt                                                        >  1:5685/1‑149 (MQ=255)
                                                                                           ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAt                                                        >  1:207676/1‑149 (MQ=255)
                                                                                            cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa                                                                                      <  4:152096/118‑1 (MQ=255)
                                                                                            cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa                                                                                      >  3:152096/1‑118 (MQ=255)
                                                                                            cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAtt                                                       >  3:69864/1‑149 (MQ=255)
                                                                                            cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAtt                                                       >  7:505999/1‑149 (MQ=255)
                                                                                             aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                      >  7:18789/1‑149 (MQ=255)
                                                                                              gCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgc                                                     <  2:207676/149‑1 (MQ=255)
                                                                                               ccATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgcg                                                    >  7:356142/1‑149 (MQ=255)
                                                                                                  ttGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                 <  8:58048/149‑1 (MQ=255)
                                                                                                  ttGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                 >  2:183177/1‑149 (MQ=255)
                                                                                                    gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                          >  5:189894/1‑122 (MQ=255)
                                                                                                    gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                          <  6:189894/122‑1 (MQ=255)
                                                                                                     cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGAc                                                                                                                                     <  1:166054/62‑1 (MQ=255)
                                                                                                     cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGAc                                                                                                                                     >  2:166054/1‑62 (MQ=255)
                                                                                                      gACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATa                                                                            <  2:497724/118‑1 (MQ=255)
                                                                                                      gACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATa                                                                            >  1:497724/1‑118 (MQ=255)
                                                                                                      gACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcgacga                                              >  2:275655/1‑148 (MQ=255)
                                                                                                             cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                       <  4:69864/148‑1 (MQ=255)
                                                                                                             cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                       <  7:169486/148‑1 (MQ=255)
                                                                                                              aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                                                                                 <  8:348814/105‑1 (MQ=255)
                                                                                                              aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                                                                                 >  7:348814/1‑105 (MQ=255)
                                                                                                               cGCAGCGTATTCAACGCCTGACGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCgg                                    <  6:109157/149‑1 (MQ=255)
                                                                                                                   gCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCaaa                                >  4:404015/1‑149 (MQ=255)
                                                                                                                       aTTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCgg                            >  5:167042/1‑149 (MQ=255)
                                                                                                                        ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGGCGATAGTTCGCGGTTAATTATGAATTGCGCGACGATATCTTTCGGCGAATCGgc                           <  2:155156/149‑1 (MQ=255)
                                                                                                                           aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGgcgc                         <  4:206456/148‑1 (MQ=255)
                                                                                                                           aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                        >  3:441936/1‑149 (MQ=255)
                                                                                                                           aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                        >  7:457431/1‑149 (MQ=255)
                                                                                                                           aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                        >  4:421597/1‑149 (MQ=255)
                                                                                                                             cGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGt                      <  2:287463/149‑1 (MQ=255)
                                                                                                                               ccggCCGGACCTGTAACGTCTGGCTGGCCGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATGCAAAGCCGCCGATAGTTAGCGCTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGttt                    >  8:352889/4‑149 (MQ=255)
                                                                                                                                   ccGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCt                 <  4:441936/148‑1 (MQ=255)
                                                                                                                                   ccGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCt                 <  1:379415/148‑1 (MQ=255)
                                                                                                                                     ggACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCt                                         >  7:455888/1‑122 (MQ=255)
                                                                                                                                          aacagcGTCAGGCTGCCAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGAGAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCaa         >  1:440264/6‑149 (MQ=255)
                                                                                                                                           gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt        >  4:159436/1‑149 (MQ=255)
                                                                                                                                           gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt        >  3:35240/1‑149 (MQ=255)
                                                                                                                                           gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt        >  8:189081/1‑149 (MQ=255)
                                                                                                                                           acagcGTCAGGCTGCCAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                 >  3:10500/5‑108 (MQ=255)
                                                                                                                                           acagcGTCAGGCTGCCAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                 <  4:10500/104‑1 (MQ=255)
                                                                                                                                            tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCACCGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAGTTTGCATTGCGCGACGATATCGTTCGGCAAATCGGCGGTGGTTTCTGCCAGCACtt       >  6:283711/1‑149 (MQ=255)
                                                                                                                                            tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCGGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGGAAGCAAta       >  4:259185/1‑148 (MQ=255)
                                                                                                                                            tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAtt       >  6:82959/1‑149 (MQ=255)
                                                                                                                                            tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAGTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCTCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAtt       >  6:276794/1‑149 (MQ=255)
                                                                                                                                              aCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCg     <  3:454922/149‑1 (MQ=255)
                                                                                                                                                 tCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTa  >  5:338803/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTA  >  NC_000913/4388349‑4388642

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