Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,347,285 C→G 5.3% M14I (ATG→ATC 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,2850CG 5.3% 269.9 / 5.0 75M14I (ATG→ATCfumBanaerobic 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 = 9.94e-01

CATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGG  >  NC_000913/4347144‑4347389
                                                                                                                                             |                                                                                                        
cATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCggg                                                                                                    >  4:473529/1‑148 (MQ=255)
           cgtcgtGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                         >  3:688845/1‑148 (MQ=255)
                tGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                         <  1:577091/143‑1 (MQ=255)
                tGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                         >  2:577091/1‑143 (MQ=255)
                 gAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaa                                                                                                  <  2:186323/133‑1 (MQ=255)
                 gAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGaa                                                                                                  >  1:186323/1‑133 (MQ=255)
                       cTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAggg                                                                             <  3:473529/148‑1 (MQ=255)
                         tgctgcGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCAt                                                                                                        >  1:147735/1‑119 (MQ=255)
                         tgctgcGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCAt                                                                                                        <  2:147735/119‑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                                                                                    <  3:1076957/129‑1 (MQ=255)
                                   gcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGa                                                                                    >  4:1076957/1‑129 (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                                                                                               <  4:1097274/107‑1 (MQ=255)
                                              gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAgg                                                                                               >  3:1097274/1‑107 (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                                                                        >  3:182780/1‑115 (MQ=255)
                                                             ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgt                                                                        <  4:182780/115‑1 (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                                    <  3:670739/140‑1 (MQ=255)
                                                                        ttCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGtt                                    >  4:670739/1‑140 (MQ=255)
                                                                            ccGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTa                                          >  3:765234/1‑130 (MQ=255)
                                                                            ccGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTa                                          <  4:765234/130‑1 (MQ=255)
                                                                                            cGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATa                                                                <  4:955765/92‑1 (MQ=255)
                                                                                            cGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATa                                                                >  3:955765/1‑92 (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                                                                         >  1:1120190/1‑79 (MQ=255)
                                                                                                aaCGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttg                                                                         <  2:1120190/79‑1 (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                                                                     >  4:1095871/1‑77 (MQ=255)
                                                                                                      aTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGtttgtttg                                                                     <  3:1095871/77‑1 (MQ=255)
                                                                                                       tCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcgg  >  4:358959/1‑143 (MQ=255)
                                                                                                       tCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcgg  <  3:358959/143‑1 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCtt                                                            <  1:707952/84‑1 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCtt                                                            >  2:707952/1‑84 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTa                                                  <  4:416240/94‑1 (MQ=255)
                                                                                                        cGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTa                                                  >  3:416240/1‑94 (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                          <  2:998866/94‑1 (MQ=255)
                                                                                                                                tATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAAt                          >  1:998866/1‑94 (MQ=255)
                                                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcg                 <  3:121483/102‑1 (MQ=255)
                                                                                                                                 aTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcg                 >  4:121483/1‑102 (MQ=255)
                                                                                                                                             |                                                                                                        
CATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGG  >  NC_000913/4347144‑4347389

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: