Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,435,589 G→A H1074Y (CAT→TAT)  ydbK ← putative oxidoreductase Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,5890GA95.0% 379.1 / 3.9 121H1074Y (CAT→TAT) ydbKputative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (64/50);  total (66/55)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 121 total aligned reads displayed.

GCACGACGCGGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGCTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGCTGCG  >  NC_000913/1435498‑1435678
                                                                                           |                                                                                         
gCACGACGCGGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGc                                                                                          <  2:3452201/93‑1 (MQ=255)
  acgacgCGGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATActc                                                                                        >  1:3087012/1‑93 (MQ=255)
    gacgCGGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCtt                                                                                      >  2:1544253/1‑93 (MQ=255)
     acgCGGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTc                                                                                     <  1:1749151/93‑1 (MQ=255)
        cgGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTcaca                                                                                  >  1:3969992/1‑93 (MQ=255)
        cgGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTcaca                                                                                  >  1:4361017/1‑93 (MQ=255)
          gctCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACg                                                                                <  1:172935/91‑1 (MQ=255)
          gATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACg                                                                                >  1:4204740/1‑93 (MQ=255)
          gATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACg                                                                                >  1:1831323/1‑93 (MQ=255)
          gATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACg                                                                                <  2:3838123/93‑1 (MQ=255)
           aTCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACgg                                                                               >  2:2843710/1‑93 (MQ=255)
           aTCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACgg                                                                               <  1:3485257/93‑1 (MQ=255)
           aTCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACgg                                                                               >  1:434072/1‑93 (MQ=255)
            tCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACggg                                                                              >  1:4176098/1‑93 (MQ=255)
            tCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACggg                                                                              >  1:1354933/1‑93 (MQ=255)
             cAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGc                                                                             <  1:1637275/93‑1 (MQ=255)
              aaaGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCt                                                                            >  2:880044/1‑93 (MQ=255)
               aaGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCta                                                                           >  2:2116636/1‑93 (MQ=255)
                 gCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGCTCTTCACACGGGCtata                                                                         <  2:3389625/93‑1 (MQ=255)
                  cGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATaa                                                                        <  1:3954410/93‑1 (MQ=255)
                   gATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAg                                                                       >  2:3966528/1‑93 (MQ=255)
                   gATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAg                                                                       <  2:4646449/93‑1 (MQ=255)
                    atatAGCGGCCCGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGc                                                                      <  1:3229907/93‑1 (MQ=255)
                    atatAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAga                                                                      >  1:3671635/1‑92 (MQ=255)
                    atatAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGc                                                                      <  1:617587/93‑1 (MQ=255)
                     tataGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCa                                                                     >  2:3206238/1‑93 (MQ=255)
                     tataGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCa                                                                     >  2:2885893/1‑93 (MQ=255)
                      ataGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCaa                                                                    <  1:4398054/93‑1 (MQ=255)
                       taGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGCTCTTCACACGGGCTATAAGCAAt                                                                   >  2:4318280/1‑93 (MQ=255)
                       taGCGACCAGAAGCCGGTAGCTGTGAGTTGGCGCAGCTGGTCGTGGCTGAGTGCCAGATCGTACCCATACTCTTCACACGGGCTATAAGCAAt                                                                   <  2:1124397/93‑1 (MQ=255)
                         gCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGa                                                                 >  2:793533/1‑93 (MQ=255)
                          cGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGAt                                                                <  1:3287355/93‑1 (MQ=255)
                          cGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGAt                                                                <  1:631612/93‑1 (MQ=255)
                           ggCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATc                                                               >  2:1989523/1‑93 (MQ=255)
                           ggCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATc                                                               <  2:1817773/93‑1 (MQ=255)
                           ggCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATc                                                               <  2:179652/93‑1 (MQ=255)
                            gCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCGAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCa                                                              <  1:3741220/93‑1 (MQ=255)
                            gCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCa                                                              <  1:2481255/93‑1 (MQ=255)
                            gCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCa                                                              <  1:4490299/93‑1 (MQ=255)
                             ccAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAg                                                             <  1:4313632/93‑1 (MQ=255)
                              cAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGc                                                            >  1:3851487/1‑93 (MQ=255)
                              cAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGc                                                            <  1:3947310/93‑1 (MQ=255)
                               aGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCg                                                           >  1:197235/1‑93 (MQ=255)
                                gAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGa                                                          <  2:3739004/93‑1 (MQ=255)
                                gAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGa                                                          >  2:1060364/1‑93 (MQ=255)
                                 aaGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGCTCTTCACACGGGCTATAAGCAATGATCAGCGAt                                                         >  1:2110389/1‑93 (MQ=255)
                                 aaGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGAt                                                         <  2:3506410/93‑1 (MQ=255)
                                  aGCCGGTAGCTGTGAGTTGGCGCATCTGTTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATg                                                        <  1:1226209/93‑1 (MQ=255)
                                  aGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGCTCTTCACACGGGCTATAAGCAATGATCAGCGATg                                                        <  2:3943332/93‑1 (MQ=255)
                                    ccGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGc                                                      <  1:4095354/93‑1 (MQ=255)
                                       gTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGcccc                                                   >  2:4422697/1‑93 (MQ=255)
                                           cTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATGGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGaa                                               >  1:3849437/1‑92 (MQ=255)
                                            tgtgAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAGCCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATa                                              <  2:766458/93‑1 (MQ=255)
                                            tgtgAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATa                                              <  1:734137/93‑1 (MQ=255)
                                            tgtgAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATa                                              <  1:1848806/93‑1 (MQ=255)
                                            tgtgAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATa                                              >  1:1437262/1‑93 (MQ=255)
                                             gtgAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAc                                             >  2:3614207/1‑93 (MQ=255)
                                               gAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATActc                                           <  2:2959409/93‑3 (MQ=255)
                                               gAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgc                                           <  1:433549/93‑1 (MQ=255)
                                               gAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgc                                           <  2:3052366/93‑1 (MQ=255)
                                               gAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgc                                           <  2:3618249/93‑1 (MQ=255)
                                                 gTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgctt                                         <  1:3352546/93‑1 (MQ=255)
                                                  gtGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgcttc                                        <  2:1038736/92‑1 (MQ=255)
                                                     gcgcATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgcttccgc                                     <  2:2946486/93‑1 (MQ=255)
                                                        cATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAcgcttccgcttc                                  >  2:2480406/1‑93 (MQ=255)
                                                         aTCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACgcttccgcttcc                                 >  2:177654/1‑93 (MQ=255)
                                                         aTCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACgcttccgcttcc                                 >  2:2445837/1‑93 (MQ=255)
                                                         aTCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACgcttccgcttcc                                 <  1:3148117/93‑1 (MQ=255)
                                                         aTCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACgcttccgcttcc                                 <  2:389188/93‑1 (MQ=255)
                                                          tCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCt                                <  2:391013/93‑1 (MQ=255)
                                                          tCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCt                                >  1:4101141/1‑93 (MQ=255)
                                                          tCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCt                                >  1:2776727/1‑93 (MQ=255)
                                                          tCTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCt                                >  2:629412/1‑93 (MQ=255)
                                                           cTGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTg                               >  1:2749641/1‑93 (MQ=255)
                                                            tGGTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGa                              >  2:2744523/1‑93 (MQ=255)
                                                             ggTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGaa                             <  1:361644/93‑1 (MQ=255)
                                                              gTCGTGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAAt                            >  2:3439091/1‑93 (MQ=255)
                                                                  tGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCt                        >  1:3076303/1‑93 (MQ=255)
                                                                  tGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCt                        <  2:2818616/93‑1 (MQ=255)
                                                                  tGGCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCt                        <  2:1080681/93‑1 (MQ=255)
                                                                   ggCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCtt                       <  2:3560271/93‑1 (MQ=255)
                                                                   ggCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCtt                       >  2:1214801/1‑93 (MQ=255)
                                                                   ggCTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCtt                       <  1:4516725/93‑1 (MQ=255)
                                                                     cTGAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTc                     >  1:2578643/1‑93 (MQ=255)
                                                                       gAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCAc                   >  1:2521273/1‑93 (MQ=255)
                                                                       gAGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCAc                   >  2:4540337/1‑93 (MQ=255)
                                                                       gAGCGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGACCCCGAATACGCTTCCGCTTCCTGAATCGCTTTCAc                   >  2:3006311/1‑93 (MQ=255)
                                                                        aGTGCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCAcc                  >  1:3753145/1‑93 (MQ=255)
                                                                           gCCAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTc               >  2:205532/1‑93 (MQ=255)
                                                                            ccAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCt              <  1:3969180/93‑1 (MQ=255)
                                                                            ccAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCt              >  1:1007522/1‑93 (MQ=255)
                                                                             cAGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTg             >  2:3768202/1‑93 (MQ=255)
                                                                              aGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTgg            >  1:253688/1‑93 (MQ=255)
                                                                              aGATCGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTgg            >  1:3912524/1‑93 (MQ=255)
                                                                                  cGTAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCa        <  1:3937327/93‑1 (MQ=255)
                                                                                    tAACCATCCTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGc      <  2:2729150/93‑1 (MQ=255)
                                                                                    tAACCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGc      <  1:2789946/93‑1 (MQ=255)
                                                                                      aCCATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGCTg    >  2:849850/1‑93 (MQ=255)
                                                                                        cATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATAGGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGCTgcg  >  1:2685840/1‑93 (MQ=255)
                                                                                        cATACTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGCTgcg  >  1:1116530/1‑93 (MQ=255)
                                                                                           |                                                                                         
GCACGACGCGGATCAAAGCGATATAGCGGCCAGAAGCCGGTAGCTGTGAGTTGGCGCATCTGGTCGTGGCTGAGTGCCAGATCGTAACCATGCTCTTCACACGGGCTATAAGCAATGATCAGCGATGGCCCCGGATACGCTTCCGCTTCCTGAATCGCTTTCACCGTCTGGTTCAGCTGCG  >  NC_000913/1435498‑1435678

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