Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,347,288 C→G 5.3% P13P (CCG→CCC fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,347,2880CG 5.3% 273.5 / 5.3 75P13P (CCG→CCCfumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base G (2/2);  total (37/38)
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.45e-01

CATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAG  >  NC_000913/4347144‑4347431
                                                                                                                                                |                                                                                                                                               
cATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCggg                                                                                                                                              >  4:473529/1‑148 (MQ=255)
           cgtcgtGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                   >  3:688845/1‑148 (MQ=255)
                tGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                   >  2:577091/1‑143 (MQ=255)
                tGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                   <  1:577091/143‑1 (MQ=255)
                 gAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaa                                                                                                                                            >  1:186323/1‑133 (MQ=255)
                 gAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaa                                                                                                                                            <  2:186323/133‑1 (MQ=255)
                       cTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAggg                                                                                                                       <  3:473529/148‑1 (MQ=255)
                          gctgcGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgttt                                                                                                                    <  4:688845/148‑1 (MQ=255)
                                caagcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGa                                                                                                              <  2:452607/146‑1 (MQ=255)
                                  agcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCGATGGGGAAAg                                                                                                                                          >  3:986138/1‑118 (MQ=255)
                                  agcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCGATGGGGAAAg                                                                                                                                          <  4:986138/118‑1 (MQ=255)
                                   gcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGa                                                                                                                              >  4:1076957/1‑129 (MQ=255)
                                   gcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGa                                                                                                                              <  3:1076957/129‑1 (MQ=255)
                                          cAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACAt                                                                                                           <  3:400163/141‑1 (MQ=255)
                                          cAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACAt                                                                                                           >  4:400163/1‑141 (MQ=255)
                                            gggCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAgg                                                                                                                        <  4:916064/126‑1 (MQ=255)
                                            gggCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAgg                                                                                                                        >  3:916064/1‑126 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaa                                                                                                                                            >  1:779467/1‑104 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaa                                                                                                                                            <  2:779467/104‑1 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAgg                                                                                                                                         >  3:1097274/1‑107 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAgg                                                                                                                                         <  4:1097274/107‑1 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                   <  4:991927/113‑1 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                   >  3:991927/1‑113 (MQ=255)
                                               cTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgt                                                                                                                      >  4:388331/1‑125 (MQ=255)
                                               cTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgt                                                                                                                      <  3:388331/125‑1 (MQ=255)
                                                           cTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTAt                                                                                   >  3:930081/1‑148 (MQ=255)
                                                             ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgt                                                                                                                  <  4:182780/115‑1 (MQ=255)
                                                             ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgt                                                                                                                  >  3:182780/1‑115 (MQ=255)
                                                             ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTAttt                                                                                 >  2:180016/1‑148 (MQ=255)
                                                              tCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTg                                                                                <  1:180016/148‑1 (MQ=255)
                                                                  gATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgt                                                                                                                  <  4:770796/110‑1 (MQ=255)
                                                                  gATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgt                                                                                                                  >  3:770796/1‑110 (MQ=255)
                                                                  gATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTAt                                                                            <  4:930081/148‑1 (MQ=255)
                                                                   aTGGTTTCGCTGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaaa                                                                                                                                           <  3:978984/84‑1 (MQ=255)
                                                                   aTGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaaa                                                                                                                                           >  4:978984/1‑84 (MQ=255)
                                                                   aTGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgtt                                                                                                                 >  3:1040693/1‑110 (MQ=255)
                                                                   aTGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgtt                                                                                                                 <  4:1040693/110‑1 (MQ=255)
                                                                        ttCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGtt                                                                              >  4:670739/1‑140 (MQ=255)
                                                                        ttCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGtt                                                                              <  3:670739/140‑1 (MQ=255)
                                                                            ccGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTa                                                                                    >  3:765234/1‑130 (MQ=255)
                                                                            ccGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTa                                                                                    <  4:765234/130‑1 (MQ=255)
                                                                                            cGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATa                                                                                                          >  3:955765/1‑92 (MQ=255)
                                                                                            cGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATa                                                                                                          <  4:955765/92‑1 (MQ=255)
                                                                                              cTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTAt                                                                            >  3:796354/1‑120 (MQ=255)
                                                                                              cTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTAt                                                                            <  4:796354/120‑1 (MQ=255)
                                                                                               tAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGt                                                        <  2:836185/139‑1 (MQ=255)
                                                                                               tAACGTAATCAGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGt                                                        >  1:836185/1‑139 (MQ=255)
                                                                                                aaCGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttg                                                                                                                   <  2:1120190/79‑1 (MQ=255)
                                                                                                aaCGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttg                                                                                                                   >  1:1120190/1‑79 (MQ=255)
                                                                                                aaCGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACgg                                                >  3:484093/1‑146 (MQ=255)
                                                                                                aaCGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACgg                                                <  4:484093/146‑1 (MQ=255)
                                                                                                    tAATCGGAAGTGAGTAGATAGTATTCGGTATTGCCTTTCCCGATGGGGAAAga                                                                                                                                         >  3:463580/1‑52 (MQ=25)
                                                                                                    tAATCGGAAGTGAGTAGATAGTATTCGGTATTGCCTTTCCCGATGGGGAAAga                                                                                                                                         <  4:463580/53‑2 (MQ=25)
                                                                                                      aTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGtttgtttg                                                                                                               <  3:1095871/77‑1 (MQ=255)
                                                                                                      aTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGtttgtttg                                                                                                               >  4:1095871/1‑77 (MQ=255)
                                                                                                       tCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcgg                                            <  3:358959/143‑1 (MQ=255)
                                                                                                       tCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcgg                                            >  4:358959/1‑143 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCtt                                                                                                      <  1:707952/84‑1 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCtt                                                                                                      >  2:707952/1‑84 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTa                                                                                            >  3:416240/1‑94 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTa                                                                                            <  4:416240/94‑1 (MQ=255)
                                                                                                           aaGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAggg                                                                                                                       >  1:420836/1‑64 (MQ=255)
                                                                                                           aaGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAggg                                                                                                                       <  2:420836/64‑1 (MQ=255)
                                                                                                            aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTAtt                                                                                  >  1:1079239/1‑100 (MQ=255)
                                                                                                            aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTAtt                                                                                  <  2:1079239/100‑1 (MQ=255)
                                                                                                            aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGt                                                        <  2:47906/126‑1 (MQ=255)
                                                                                                            aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGt                                                        >  1:47906/1‑126 (MQ=255)
                                                                                                                     aTAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCg                                                                >  4:650724/1‑109 (MQ=255)
                                                                                                                     aTAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCg                                                                <  3:650724/109‑1 (MQ=255)
                                                                                                                                tATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAAt                                                                    >  1:998866/1‑94 (MQ=255)
                                                                                                                                tATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAAt                                                                    <  2:998866/94‑1 (MQ=255)
                                                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcg                                                           >  4:121483/1‑102 (MQ=255)
                                                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcg                                                           <  3:121483/102‑1 (MQ=255)
                                                                                                                                               tCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAg  >  4:866058/1‑145 (MQ=255)
                                                                                                                                               tCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAg  <  3:866058/145‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAG  >  NC_000913/4347144‑4347431

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