Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,356,737 T→C 15.4% intergenic (+26/‑167) glpC → / → yfaD anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster/transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,356,7370TC15.4% 312.5 / 43.1 134intergenic (+26/‑167)glpC/yfaDanaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base T (63/50);  new base C (4/17);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 3.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 142 total aligned reads displayed.

ATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAAC  >  NC_000913/2356639‑2356835
                                                                                                  |                                                                                                  
aTGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtg                                                                                                   <  2:2755350/100‑1 (MQ=255)
 tGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                  <  2:6211/100‑1 (MQ=255)
 tGGCAGATTGAGATGTCCACAAGTCTGCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                  <  2:3350826/100‑1 (MQ=255)
  ggCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                 >  1:1334955/1‑100 (MQ=255)
  ggCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                 <  1:3123399/100‑1 (MQ=255)
  ggCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                 >  1:2598741/1‑100 (MQ=255)
  ggCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                 >  1:2301237/1‑100 (MQ=255)
   gCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTgg                                                                                                <  1:3000823/100‑1 (MQ=255)
   gCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTgg                                                                                                >  1:2837800/1‑100 (MQ=255)
      gATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTa                                                                                             >  2:2793606/1‑100 (MQ=255)
         tGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCg                                                                                          <  1:280121/100‑1 (MQ=255)
          gagaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCg                                                                                          >  2:1972781/1‑99 (MQ=255)
          gagaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCgg                                                                                         <  2:1775050/100‑1 (MQ=39)
          gagaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTAATGCCCGGTAACCACATGCTTACCgg                                                                                         <  2:271884/100‑1 (MQ=38)
           agaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCggg                                                                                        <  2:477304/100‑1 (MQ=38)
           agaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCggg                                                                                        <  1:400618/100‑1 (MQ=38)
           agaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCggg                                                                                        <  2:885907/100‑1 (MQ=38)
           agaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCggg                                                                                        <  1:1261288/100‑1 (MQ=38)
            gaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCGGGc                                                                                       <  1:1703795/100‑1 (MQ=38)
             agggccACTAGTCCTCGCTGCGACCATCCGACTACGCGACTGGCCCCGGCGCTGGCTTAAACTCCGTTCTGATGCCCGGTAACCACATGCTTACCggg                                                                                        <  2:3053363/94‑1 (MQ=19)
             aTGTCCACATGTCTTCGCTGCGAACATCCGATTACGCTACTGTCCCACGCCCTGGCTTAAACTCCTTCCTGATGCCCGGTAACCACATGCTTACCggg                                                                                        <  2:1342228/98‑1 (MQ=37)
             aTGTCCACAGGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCGTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACTggg                                                                                        <  2:1967394/98‑1 (MQ=255)
              tGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAt                                                                                     >  2:1200741/1‑100 (MQ=255)
               gTCCACAAGTCTTCGCTGCGAACATCAGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCa                                                                                      <  2:245287/98‑1 (MQ=255)
                 ccACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAtttt                                                                                  >  1:2411283/1‑100 (MQ=255)
                 ccACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGGTACCGGGCAgttt                                                                                  >  1:1982899/1‑100 (MQ=255)
                 ccACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCGGGCAtcag                                                                                  <  2:1119315/100‑4 (MQ=37)
                  cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                 >  1:281528/1‑100 (MQ=255)
                   acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTg                                                                                <  1:1591473/100‑1 (MQ=255)
                   acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTg                                                                                <  2:2207843/100‑1 (MQ=255)
                   acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTg                                                                                <  2:2771213/100‑1 (MQ=255)
                   aaaaGTCTTCGCTGCGAACATCAGATTACGCTACTGGCCCAGGAGCTGGCTTAACCTTCTTTCTGAGGCCCGGTAAGCATGTGGTTACCGGTCATTTTTg                                                                                <  2:2119781/98‑1 (MQ=255)
                    caAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGc                                                                               >  1:691508/1‑100 (MQ=255)
                    caAGTCGTCGCTGCGAACTTCCGATTATGCTACTGGCCCAGGCTATGGCTTAAACGCGTTTCTGATGACCGGTAAGCATGTGGTTACCGGGCATTTTTga                                                                               <  2:1338654/100‑2 (MQ=255)
                     aaGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCg                                                                              <  2:1463053/100‑1 (MQ=255)
                     aaGTCTTCGCTGAGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTATACTTCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCg                                                                              <  2:1204870/100‑1 (MQ=255)
                      aGTCTTCGCTGCGAACATCCGATTACGCTACTGGCACAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCg                                                                              <  2:2868246/99‑1 (MQ=255)
                       gTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTa                                                                            >  2:764802/1‑100 (MQ=255)
                         cTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTaca                                                                          >  1:413567/1‑100 (MQ=255)
                               tGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTc                                                                    >  1:2268064/1‑100 (MQ=255)
                                     cATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGcc                                                              >  1:1003419/1‑100 (MQ=255)
                                     cATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGcc                                                              <  1:1056425/100‑1 (MQ=255)
                                      aTCCGATTACGCTACTGGCCCATGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTg                                                                <  2:1062522/97‑1 (MQ=255)
                                      aTCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTCCACGATTCCGTGccc                                                             >  1:3582598/1‑100 (MQ=255)
                                       tCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCa                                                            >  1:2973138/1‑100 (MQ=255)
                                       tCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTACTGGCCCa                                                            >  2:3262686/1‑100 (MQ=255)
                                         cGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAAt                                                          >  1:1822475/1‑100 (MQ=255)
                                          gATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATg                                                         >  1:2029177/1‑100 (MQ=255)
                                              aCGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATg                                                     >  1:1229042/1‑100 (MQ=255)
                                               cGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGc                                                    >  1:868016/1‑100 (MQ=255)
                                               cGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGc                                                    >  1:1702221/1‑100 (MQ=255)
                                               cGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGc                                                    <  1:1476586/100‑1 (MQ=255)
                                                gCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCg                                                   <  2:3287549/100‑1 (MQ=255)
                                                 cTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGt                                                  >  1:2508973/1‑100 (MQ=255)
                                                 cTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGt                                                  >  2:938368/1‑100 (MQ=255)
                                                  tACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGtt                                                 >  1:2448075/1‑100 (MQ=255)
                                                    cTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGc                                               >  1:243063/1‑100 (MQ=255)
                                                     tGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCa                                              <  1:2638944/100‑1 (MQ=255)
                                                      ggCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCaa                                             >  1:3110571/1‑100 (MQ=255)
                                                        cccAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACg                                           >  1:3176882/1‑100 (MQ=255)
                                                         ccAGACGCTGGCTTAACCTCCTTTCTGATGCCCGGTAAGCATGTGGTTCCCGGGCATTTTTGCGTCCACGATCCCGTGCCCAATGTATGCGTTGCAACGc                                          <  2:1442562/100‑1 (MQ=255)
                                                              cgcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAAAGCACTg                                      >  2:3024877/1‑99 (MQ=255)
                                                               gcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaa                                    >  2:1705931/1‑100 (MQ=255)
                                                               gcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaa                                    <  2:2950073/100‑1 (MQ=255)
                                                                 tGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaaaa                                  <  1:180316/100‑1 (MQ=255)
                                                                  ggCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCatacattgg                                 <  2:3437041/100‑9 (MQ=255)
                                                                    cTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                               >  1:3487012/1‑100 (MQ=255)
                                                                      tAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                             >  2:301005/1‑100 (MQ=255)
                                                                       aaaCTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCctc                            <  2:379246/100‑1 (MQ=255)
                                                                       aaaCTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCctc                            <  1:301007/100‑1 (MQ=255)
                                                                          cTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGa                         >  1:3161918/1‑100 (MQ=255)
                                                                           tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGaa                        <  1:2541608/100‑1 (MQ=255)
                                                                           tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGaa                        <  1:571924/100‑1 (MQ=255)
                                                                            ccTTTCTTATGCCCAGTAAGAATGTAGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGaaa                       >  2:480317/1‑100 (MQ=255)
                                                                            ccTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGaaa                       >  2:2835896/1‑100 (MQ=255)
                                                                             cTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGaaaa                      <  2:3110548/100‑1 (MQ=255)
                                                                               ttCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACg                    <  1:2346420/100‑1 (MQ=255)
                                                                               caaaaaTGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCAGatcgga                    >  1:1055726/6‑95 (MQ=25)
                                                                                  tGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtct                 >  1:2203075/1‑100 (MQ=255)
                                                                                  tGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtct                 >  2:3082688/1‑100 (MQ=255)
                                                                                   gATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtctc                <  1:315591/100‑1 (MQ=255)
                                                                                   gATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGAGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtct                 >  1:940/1‑99 (MQ=255)
                                                                                    aTGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCg               <  2:832303/100‑1 (MQ=255)
                                                                                    aTGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCg               >  1:334233/1‑100 (MQ=38)
                                                                                      gCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCa             <  1:3350831/100‑1 (MQ=255)
                                                                                      gCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCa             >  2:1137155/1‑100 (MQ=38)
                                                                                       cccGGTAATCATGTGGTTACCGGGCAGTTTTGCTTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCa             <  2:1425043/99‑1 (MQ=255)
                                                                                       cccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaa            >  1:2201516/1‑100 (MQ=255)
                                                                                       cccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaa            >  1:3606644/1‑100 (MQ=255)
                                                                                        ccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATTCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaaa           >  2:2563094/1‑100 (MQ=255)
                                                                                        ccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaaa           >  1:521763/1‑100 (MQ=255)
                                                                                        ccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaaa           <  1:2854756/100‑1 (MQ=255)
                                                                                         cGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAg          >  1:669142/1‑100 (MQ=255)
                                                                                          ggTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAgg         <  2:3403266/100‑1 (MQ=255)
                                                                                            tAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGc        <  1:3295312/99‑1 (MQ=255)
                                                                                               gCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaa    >  1:2080033/1‑100 (MQ=255)
                                                                                               gCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaa    <  1:657028/100‑1 (MQ=255)
                                                                                               gCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaa    >  2:1915219/1‑100 (MQ=255)
                                                                                                cATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaaa   >  1:3291815/1‑100 (MQ=255)
                                                                                                 aTGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAAc  <  2:3616812/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
ATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAAC  >  NC_000913/2356639‑2356835

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: