Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,388,960 C→T S137S (TCG→TCA mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,9600CT99.2% 414.9 / ‑4.7 123S137S (TCG→TCAmscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (57/65);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 123 total aligned reads displayed.

GCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGT  >  NC_000913/4388815‑4389105
                                                                                                                                                 |                                                                                                                                                 
gcgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAAt                                                                                                                                               >  1:395709/1‑150 (MQ=255)
gcgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATc                                                                                                                                              <  2:1558117/151‑1 (MQ=255)
 tgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTGAGTGaa                                                                                                                                                >  2:786430/2‑148 (MQ=255)
  gcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCgg                                                                                                                                            <  4:1403639/151‑1 (MQ=255)
      tAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCaa                                                                                                                                         >  4:79911/1‑150 (MQ=255)
      tAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCa                                                                                                                                          <  2:690740/149‑1 (MQ=255)
       aaGACGTGCAGAGCCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGGATCGGCAAt                                                                                                                                        >  4:1364330/1‑150 (MQ=255)
        aGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAAtct                                                                                                                                      <  4:1317361/151‑1 (MQ=255)
         gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAAtct                                                                                                                                      >  2:989577/1‑150 (MQ=255)
           cGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTcgc                                                                                                                                   <  4:723384/151‑1 (MQ=255)
           cGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTcgc                                                                                                                                   <  4:737748/151‑1 (MQ=255)
            gTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGGATCGGCAATCTcgc                                                                                                                                   >  4:1588577/1‑150 (MQ=255)
            gTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTcgcg                                                                                                                                  >  3:220354/1‑151 (MQ=255)
            gTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTcgcg                                                                                                                                  <  4:1286820/151‑1 (MQ=255)
              gCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCgg                                                                                                                                 <  4:94638/150‑1 (MQ=255)
              gCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCgg                                                                                                                                 <  4:874177/150‑1 (MQ=255)
                agagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGgcg                                                                                                                              >  4:1535549/1‑151 (MQ=255)
                agagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGgcg                                                                                                                              <  4:1238927/151‑1 (MQ=255)
                 gagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGgcgc                                                                                                                             <  1:259962/151‑1 (MQ=255)
                 gagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGgcg                                                                                                                              >  4:335975/1‑150 (MQ=255)
                     cAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGGTTCAGTGAATCGGCAATCTCGCGGgggcg                                                                                                                            >  2:1113273/1‑148 (MQ=255)
                     cAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCt                                                                                                                          <  4:1524305/150‑1 (MQ=255)
                     cAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCt                                                                                                                          <  4:32693/150‑1 (MQ=255)
                      aGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGGTTGCTGTTGCGGCAGCTGATTCAGTGAATCGGCGATCTCGCGGGCGCGCTcc                                                                                                                        >  2:82221/1‑151 (MQ=255)
                          tGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCCCGGGCGCGCTCCt                                                                                                                       >  4:1034212/1‑148 (MQ=255)
                           gCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGctc                                                                                                                   >  4:471269/1‑151 (MQ=255)
                            cAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGctct                                                                                                                  <  2:1679637/151‑1 (MQ=255)
                            cAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGctct                                                                                                                  <  4:1512829/151‑1 (MQ=255)
                            cAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGctct                                                                                                                  <  2:650279/151‑1 (MQ=255)
                              aCGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCtt                                                                                                                 >  1:28736/1‑150 (MQ=255)
                              aCGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTg                                                                                                                <  2:307943/151‑1 (MQ=255)
                              aCGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTg                                                                                                                <  1:142008/151‑1 (MQ=255)
                               cgcgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGc                                                                                                               >  4:850877/1‑151 (MQ=255)
                                 cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCtt                                                                                                                 >  1:172603/1‑147 (MQ=255)
                                    aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGc                                                                                                          <  4:423995/151‑1 (MQ=255)
                                       tttGTGCCTGATTGAGCGGATTTGTGCCGGTAAGCGTACCCAGGCGGCGCTCGATCTCATTTAGCGGGCGGCGGGCGTCGCTTTGCTGTTGCGGCTGTTGATTCAGTCAATCGTCAATCTCGGGGGCGCGCTCCTCCTctggctggtcgt                                                                                                        <  2:1036871/150‑5 (MQ=255)
                                       tttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGcc                                                                                                         >  1:106181/1‑149 (MQ=255)
                                         tgtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTTCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGa                                                                                                      <  2:214418/150‑1 (MQ=255)
                                         tgtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGa                                                                                                      <  2:233474/150‑1 (MQ=255)
                                         tgtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                                                     >  4:604098/1‑151 (MQ=255)
                                          gtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAca                                                                                                    <  2:7847/151‑2 (MQ=255)
                                          gtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                                                     >  4:1667561/1‑150 (MQ=255)
                                             ccTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACgg                                                                                                   >  1:391258/1‑149 (MQ=255)
                                              cTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtt                                                                                                >  4:588754/1‑151 (MQ=255)
                                                gATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                                              <  2:1381294/151‑1 (MQ=255)
                                                gATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                                              <  2:171759/151‑1 (MQ=255)
                                                gATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCttt                                                                                               <  4:619195/150‑1 (MQ=255)
                                                    gAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTTCTGGGCCTGACGGCTTTTATcc                                                                                          >  2:105351/1‑151 (MQ=255)
                                                      gCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATcc                                                                                          >  1:198098/1‑149 (MQ=255)
                                                        ggAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTTAGTGAATCGGCAATCTCGCGGGGGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAg                                                                                        >  1:477187/1‑149 (MQ=255)
                                                        ggAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGc                                                                                       <  1:383921/150‑1 (MQ=255)
                                                          aGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCaa                                                                                     >  2:776049/1‑150 (MQ=255)
                                                              ttGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgc                                                                                >  4:1038940/1‑151 (MQ=255)
                                                               tGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgc                                                                                >  4:268384/1‑150 (MQ=255)
                                                                gCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgc                                                                                >  2:772885/1‑149 (MQ=255)
                                                                  cGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgct                                                                            <  3:444031/151‑1 (MQ=255)
                                                                   ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                           <  2:670789/151‑1 (MQ=255)
                                                                   ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                           >  4:540676/1‑151 (MQ=255)
                                                                   ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                           <  2:1246021/151‑1 (MQ=255)
                                                                    gTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                           >  3:129300/1‑150 (MQ=255)
                                                                      aaGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGCTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGAc                                                                         <  2:192249/150‑1 (MQ=255)
                                                                       aGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCt                                                                       <  2:599677/151‑1 (MQ=255)
                                                                       aGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCt                                                                       <  4:709440/151‑1 (MQ=255)
                                                                        gCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTg                                                                      <  3:207193/151‑1 (MQ=255)
                                                                           ttCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                     >  4:1116529/1‑149 (MQ=255)
                                                                             cccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGgag                                                                   >  2:806585/1‑149 (MQ=255)
                                                                              ccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAt                                                                >  3:449010/1‑151 (MQ=255)
                                                                               cAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAtt                                                               <  2:399352/151‑1 (MQ=255)
                                                                                aGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAttt                                                              <  4:840882/151‑1 (MQ=255)
                                                                                    ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                           >  2:1393201/1‑150 (MQ=255)
                                                                                     gcgcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGa                                                         <  1:14047/151‑1 (MQ=255)
                                                                                          ccaTCTAATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGTCAGTTGATGCAGTGAATCGGCATTCTCTCTGGCTTGTTCCTGCTCGTGCTGGCCCTGACGCCTTTTATCAAGCAACCCGCTTCTCCCCTACAGAACTTTCTGATTCac                                                    <  2:1679028/149‑2 (MQ=255)
                                                                                          cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACTGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTc                                                      <  4:738918/149‑1 (MQ=255)
                                                                                          cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                                                       >  1:40564/1‑148 (MQ=255)
                                                                                           gATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                    >  4:1624492/1‑150 (MQ=255)
                                                                                                cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcgc                                                 <  4:787464/148‑1 (MQ=255)
                                                                                                cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcgc                                                 <  2:1569989/148‑1 (MQ=255)
                                                                                                cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgc                                                   <  4:1508852/146‑1 (MQ=255)
                                                                                                cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                    <  2:1362632/145‑1 (MQ=255)
                                                                                                cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                    <  4:80786/145‑1 (MQ=255)
                                                                                                 aTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcgc                                                 >  2:1200806/1‑147 (MQ=255)
                                                                                                     aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCgg                                            >  1:465590/1‑148 (MQ=255)
                                                                                                     aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTa                                          <  3:55189/150‑1 (MQ=255)
                                                                                                     aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAg                                         <  3:10723/151‑1 (MQ=255)
                                                                                                         gACGGCGGGCGTCGGGTTGCTGTTGCGGCAGGTGATTCAATGAATCGGCAATCTCGCGGGCGCGCTCCTGCCCTTGCTGGGGCTGACGGGTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTTCTGATTCAGCGCGTCGGTAGACa                                      >  1:284456/1‑150 (MQ=255)
                                                                                                              cgggcgTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTccc                                 >  2:1538911/1‑150 (MQ=255)
                                                                                                                  cgTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGa                            <  4:496773/151‑1 (MQ=255)
                                                                                                                   gTCGGTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGac                           >  3:373681/1‑151 (MQ=255)
                                                                                                                       gTTTGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGgaga                                                                  <  1:68764/108‑1 (MQ=255)
                                                                                                                         ttGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCt                      <  2:315223/150‑1 (MQ=255)
                                                                                                                         ttGCTGTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCt                      <  1:306415/150‑1 (MQ=255)
                                                                                                                              gTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTgcgc                  <  4:323315/149‑1 (MQ=255)
                                                                                                                              gTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGcgcg                 <  4:328613/150‑1 (MQ=255)
                                                                                                                              gTTGCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGcgcg                 <  4:812256/150‑1 (MQ=255)
                                                                                                                                 gCGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCt              <  4:1114681/150‑1 (MQ=255)
                                                                                                                                  cGGCAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCg            <  4:664957/151‑1 (MQ=255)
                                                                                                                                     cAGTTGATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa         <  2:79763/151‑1 (MQ=255)
                                                                                                                                          gATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAt    >  4:436851/1‑151 (MQ=255)
                                                                                                                                          gATTCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAt    >  4:1593576/1‑151 (MQ=255)
                                                                                                                                             tCAGTGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAtgt  <  2:849761/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGT  >  NC_000913/4388815‑4389105

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