Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 817,441 A→G Y55C (TAC→TGC)  moaB → molybdopterin biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913817,4410AG96.0% 671.8 / 11.3 202Y55C (TAC→TGC) moaBmolybdopterin biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (91/103);  total (94/108)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 210 total aligned reads displayed.

AGACGATACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGAAGGTG  >  NC_000913/817352‑817530
                                                                                         |                                                                                         
agaCGATACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACcgctgcgc                                                                                        <  2:3803469/93‑1 (MQ=255)
     atgccTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTc                                                                                   <  1:5454418/90‑1 (MQ=255)
     aTACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTc                                                                                   <  1:5901518/93‑1 (MQ=255)
     aTACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTc                                                                                   <  1:4790472/93‑1 (MQ=255)
     aTACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTc                                                                                   >  1:4257278/1‑93 (MQ=255)
      tACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTcg                                                                                  >  1:3308742/1‑93 (MQ=255)
       aCCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTcgc                                                                                 >  1:1602555/1‑93 (MQ=255)
       aCCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTcgc                                                                                 <  1:3744743/93‑1 (MQ=255)
         cTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCgcgc                                                                               >  2:5155762/1‑93 (MQ=255)
          tCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCt                                                                              >  2:4240708/1‑93 (MQ=255)
            cGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCa                                                                            <  2:4968123/93‑1 (MQ=255)
              gTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAgg                                                                          <  1:3424139/93‑1 (MQ=255)
                cACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTa                                                                        >  2:198675/1‑93 (MQ=255)
                cACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTa                                                                        <  1:5832204/93‑1 (MQ=255)
                 aCTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAGAACCGCTGCGCTATTCGCGCTCAGGTAt                                                                       <  2:4805215/93‑1 (MQ=255)
                  cTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATc                                                                      >  2:5606191/1‑93 (MQ=255)
                   tATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCt                                                                     <  1:1697763/93‑1 (MQ=255)
                   tATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCt                                                                     <  2:2802680/93‑1 (MQ=255)
                    aTCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTg                                                                    <  1:5543134/93‑1 (MQ=255)
                     tCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGc                                                                   <  1:3999682/93‑1 (MQ=255)
                      cTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCAGGTATCTGCg                                                                  <  1:3891921/93‑1 (MQ=255)
                      cTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCAATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCCGGTATCTGCg                                                                  >  2:2385947/1‑93 (MQ=255)
                          gcgATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGa                                                              >  1:1667674/1‑93 (MQ=255)
                           cgATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAACCCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGAt                                                             <  2:413373/93‑1 (MQ=255)
                           cgATTCGGCGCAAGAAGCGGGCCAGCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGAt                                                             <  2:3374045/93‑1 (MQ=255)
                            gATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATc                                                            >  2:3040424/1‑93 (MQ=255)
                              ttCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGc                                                          <  1:3464171/93‑1 (MQ=255)
                               tCGGCCCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCTCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGcc                                                         <  1:357031/93‑1 (MQ=255)
                                  gcgcAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGc                                                      <  2:4769256/93‑1 (MQ=255)
                                  gcgcAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGc                                                      >  1:2903875/1‑93 (MQ=255)
                                    gcAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGcga                                                    <  1:5769105/93‑1 (MQ=255)
                                      aagaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGcgacg                                                  >  2:2339431/1‑93 (MQ=255)
                                      aagaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGcgacg                                                  >  1:4275466/1‑93 (MQ=255)
                                       agaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCTCTCAGGTATCTGCGTGGATCGCCAGcgacga                                                 <  1:1975564/93‑1 (MQ=255)
                                       agaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGcgacga                                                 >  2:6047193/1‑93 (MQ=255)
                                       agaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGcgacga                                                 <  1:4910729/93‑1 (MQ=255)
                                       agaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGcgacga                                                 <  2:3684906/93‑1 (MQ=255)
                                        gaagCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGAt                                                <  2:2272698/93‑1 (MQ=255)
                                          agCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGt                                              <  2:5319336/93‑1 (MQ=255)
                                          agCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGt                                              <  1:362662/93‑1 (MQ=255)
                                           gCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTa                                             >  1:2180205/1‑93 (MQ=255)
                                              ggCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACaa                                          <  1:2827619/93‑1 (MQ=255)
                                              ggCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACaa                                          <  1:742787/93‑1 (MQ=255)
                                               gCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAg                                         <  1:1629181/93‑1 (MQ=255)
                                                ccATCACGTTGTCGATAAAGCCATTGTGAAAGATAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGt                                        <  2:5941737/93‑1 (MQ=255)
                                                ccATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGt                                        >  2:3232068/1‑93 (MQ=255)
                                                 cATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTg                                       <  1:1462093/93‑1 (MQ=255)
                                                  aTCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGTTTGCCATCTATGATGTACAAGTgg                                      >  2:5495162/1‑93 (MQ=255)
                                                  aTCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTgg                                      <  1:1516579/93‑1 (MQ=255)
                                                  aTCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTgg                                      <  2:458937/93‑1 (MQ=255)
                                                  aTCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTgg                                      <  2:4080189/93‑1 (MQ=255)
                                                  aTCACGTTGTCGATAAAGCCATTGTGAAAGAAAACAGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTgg                                      >  1:6083367/1‑93 (MQ=255)
                                                    cACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGtt                                    >  2:1856379/1‑92 (MQ=255)
                                                     aCGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTAt                                   >  1:1347025/1‑93 (MQ=255)
                                                        ttGTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGa                                >  1:5049379/1‑93 (MQ=255)
                                                          gTCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGAtt                              <  2:110292/93‑1 (MQ=255)
                                                          gTCGATAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGAtt                              >  2:364281/1‑93 (MQ=255)
                                                           tCGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCATGGATCGCCAGCGACGATGTACAAGTGGTATTGATTa                             <  2:4686402/93‑1 (MQ=255)
                                                            cGATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTAc                            <  1:209596/93‑1 (MQ=255)
                                                             gATAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACg                           >  2:1306015/1‑93 (MQ=255)
                                                              aTAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACgg                          >  1:1066030/1‑93 (MQ=255)
                                                              aTAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACgg                          <  2:5357601/93‑1 (MQ=255)
                                                              aTAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACgg                          <  1:72230/93‑1 (MQ=255)
                                                               tAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACggg                         >  1:2709118/1‑93 (MQ=255)
                                                               tAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACggg                         >  1:2572462/1‑93 (MQ=255)
                                                               tAAAGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACggg                         >  1:2077249/1‑93 (MQ=255)
                                                               tAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACggg                         >  1:2326002/1‑93 (MQ=255)
                                                                aaaGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggt                        <  2:5105238/93‑1 (MQ=255)
                                                                 aaGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggtg                       >  1:3102938/1‑93 (MQ=255)
                                                                  aGCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggtgg                      >  2:1315320/1‑93 (MQ=255)
                                                                  aGCCATTGCGAAAGAAAACCTCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggtgg                      <  1:4174835/93‑1 (MQ=255)
                                                                   gCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggttgt                     >  1:1986397/1‑93 (MQ=255)
                                                                   gCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggtggt                     >  1:4879514/1‑93 (MQ=255)
                                                                   gCCATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGggtggt                     >  1:1009109/1‑93 (MQ=255)
                                                                    ccATTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTa                    >  2:5670519/1‑93 (MQ=255)
                                                                      aTTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACt                  >  1:1768063/1‑93 (MQ=255)
                                                                      aTTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACt                  >  1:5331609/1‑93 (MQ=255)
                                                                      aTTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACt                  <  2:3893450/93‑1 (MQ=255)
                                                                      aTTGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGAGACGATGTACAAGTGGTATTGATTACGGGTGGTACt                  >  2:2230451/1‑93 (MQ=255)
                                                                       ttGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTg                 <  2:500075/93‑1 (MQ=255)
                                                                       ttGTGAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTg                 <  1:4470981/93‑1 (MQ=255)
                                                                        tgtgAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTgg                <  2:847565/93‑1 (MQ=255)
                                                                          tgAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGcc              >  1:3009706/1‑93 (MQ=255)
                                                                          tgAAAGAAAACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGcc              >  2:3970294/1‑93 (MQ=255)
                                                                           gaaagaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCt             >  1:1300490/1‑93 (MQ=255)
                                                                           gaaagaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCt             <  1:5178831/93‑1 (MQ=255)
                                                                            aaagaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTg            >  2:742364/1‑93 (MQ=255)
                                                                            aaagaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTg            >  2:674367/1‑93 (MQ=255)
                                                                             aagaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGa           <  2:3496574/93‑1 (MQ=255)
                                                                             aagaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGGGGTACTGGCCTgc           >  1:3986002/1‑92 (MQ=255)
                                                                              agaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGAc          >  2:5725266/1‑93 (MQ=255)
                                                                               gaaaACCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACg         <  1:5229071/93‑1 (MQ=255)
                                                                               gaaaACCGCTACGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGGGGTAATGGCCTGACg         >  1:4438198/1‑93 (MQ=255)
                                                                                 aaaCCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGa       <  1:4610049/93‑1 (MQ=255)
                                                                                  aaCCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGaa      <  2:2075542/93‑1 (MQ=255)
                                                                                  aaCCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGaa      <  1:2170821/93‑1 (MQ=255)
                                                                                  aaCCGCTGCGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGaa      <  2:4818759/93‑1 (MQ=255)
                                                                                     cgctgcgctaTTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGAAGGt   <  1:5112896/93‑1 (MQ=255)
                                                                                      gctgcgctaTTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGTTGGTACTGGCCTGACGGAAGGTg  >  1:3526639/1‑93 (MQ=255)
                                                                                      gctgcgctaTTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGAAGGTg  <  2:1640871/93‑1 (MQ=255)
                                                                                         |                                                                                         
AGACGATACCTCCGGTCACTATCTGCGCGATTCGGCGCAAGAAGCGGGCCATCACGTTGTCGATAAAGCCATTGTGAAAGAAAACCGCTACGCTATTCGCGCTCAGGTATCTGCGTGGATCGCCAGCGACGATGTACAAGTGGTATTGATTACGGGTGGTACTGGCCTGACGGAAGGTG  >  NC_000913/817352‑817530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: