Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 321,039 G→A D70N (GAC→AAC)  ykgE → putative dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913321,0390GA99.7% 1309.6 / ‑5.6 356D70N (GAC→AAC) ykgEputative dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (169/186);  total (170/186)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 378 total aligned reads displayed.

GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCA  >  NC_000913/320949‑321128
                                                                                          |                                                                                         
 gATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAacaacg                                                                                        >  2:5336503/1‑93 (MQ=255)
  aTGCTGCGGTCAGCGTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAAcaacga                                                                                       <  1:1843711/93‑1 (MQ=255)
  aTGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAAcaacga                                                                                       >  1:4572498/1‑93 (MQ=255)
   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                      >  1:8673323/1‑93 (MQ=255)
    gctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATc                                                                                     <  2:1940714/93‑1 (MQ=255)
     ctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATcc                                                                                    >  1:6111853/1‑93 (MQ=255)
      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCTCACTGGAGGATAACAACGATccc                                                                                   >  2:5563667/1‑93 (MQ=255)
      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATccc                                                                                   <  1:4788864/93‑1 (MQ=255)
        cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCat                                                                                 <  2:3867928/93‑1 (MQ=255)
        cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAAGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCACat                                                                                 >  2:5070375/1‑93 (MQ=255)
          gTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatta                                                                               >  2:474334/1‑93 (MQ=255)
          gTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatta                                                                               >  2:2516957/1‑93 (MQ=255)
             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            >  2:6297986/1‑93 (MQ=255)
             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            >  1:6240958/1‑93 (MQ=255)
             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            <  2:9379301/93‑1 (MQ=255)
             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            <  1:7158460/93‑1 (MQ=255)
             aGCCTGCGATCAATAACGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGACGCACTGGAGGATAACAACTATACCATTAttt                                                                            >  2:105753/1‑93 (MQ=255)
                cTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAc                                                                         >  2:8877139/1‑93 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        <  2:443809/93‑1 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        >  1:3499663/1‑93 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        <  2:6710947/93‑1 (MQ=255)
                  gCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCg                                                                       <  2:9881396/93‑1 (MQ=255)
                  gCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCg                                                                       >  2:9759971/1‑93 (MQ=255)
                   cGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCgg                                                                      <  1:7920196/93‑1 (MQ=255)
                    gATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggc                                                                     >  2:7500414/1‑93 (MQ=255)
                    gATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggc                                                                     <  1:6812391/93‑1 (MQ=255)
                     aTCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggct                                                                    <  1:9762243/93‑1 (MQ=255)
                      tCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctg                                                                   >  1:5114380/1‑93 (MQ=255)
                      tCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctg                                                                   >  1:8523638/1‑93 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                 >  2:63386/1‑93 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATAAAAAATCCGATCGCCGCCCAGGAGGATAACAACCACCCCATTATTTCACCgtctggc                                                                 <  2:8315992/93‑1 (MQ=255)
                          tAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctc                                                               <  1:3367008/93‑1 (MQ=255)
                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              <  2:1797415/93‑1 (MQ=255)
                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              <  2:1153139/93‑1 (MQ=255)
                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              <  1:8296635/93‑1 (MQ=255)
                            gCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCtt                                                             <  1:9021943/93‑1 (MQ=255)
                            gCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCtt                                                             <  2:1001799/93‑1 (MQ=255)
                            gCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCtt                                                             >  2:8491061/1‑93 (MQ=255)
                             cGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTg                                                            <  2:7663333/93‑1 (MQ=255)
                             cGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTg                                                            <  1:7411303/93‑1 (MQ=255)
                              ggTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGc                                                           <  2:5164542/93‑1 (MQ=255)
                              ggTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTAGc                                                           >  1:1676695/1‑93 (MQ=255)
                              ggTTATACCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGc                                                           >  1:1321365/1‑93 (MQ=255)
                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctgg                                                                  >  1:1493721/1‑85 (MQ=255)
                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTggc                                                          >  1:5444362/1‑91 (MQ=255)
                                 tataTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCAcc                                                        >  2:3630227/1‑93 (MQ=255)
                                  atatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCAcca                                                       >  1:2778397/1‑92 (MQ=255)
                                  atatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCt                                                       <  1:9311378/93‑1 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     <  1:2450284/93‑1 (MQ=255)
                                     tCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATg                                                    >  1:4122475/1‑93 (MQ=255)
                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                  <  2:9467528/93‑1 (MQ=255)
                                           aagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaa                                              <  2:2352120/93‑1 (MQ=255)
                                            agCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaa                                             <  1:94169/93‑1 (MQ=255)
                                             gCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaa                                            >  1:3923721/1‑93 (MQ=255)
                                               gATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAg                                          >  2:3025363/1‑93 (MQ=255)
                                                 ttCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGtt                                        <  2:1455762/93‑1 (MQ=255)
                                                   ccAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAc                                      <  2:4380704/93‑1 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  1:5583130/93‑1 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  1:9538470/93‑1 (MQ=255)
                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                    <  2:7028292/93‑1 (MQ=255)
                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                    <  2:7767785/93‑1 (MQ=255)
                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                    <  1:816105/93‑1 (MQ=255)
                                                       ggATGAAAAATCTGATCACCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGa                                  >  1:1252352/1‑93 (MQ=255)
                                                       ggATGAAAAAGCTGATAGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGa                                  <  1:1015168/93‑1 (MQ=255)
                                                        gATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGAc                                 <  1:7737064/93‑1 (MQ=255)
                                                         aTGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACg                                <  2:3302361/93‑1 (MQ=255)
                                                         aTGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAACCGACg                                >  1:2010006/1‑93 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTa                              >  2:4990673/1‑93 (MQ=255)
                                                            aaaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGCCGTAt                             >  2:2005201/1‑93 (MQ=255)
                                                            aaaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTAt                             <  1:1805839/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                            <  2:3221559/93‑1 (MQ=255)
                                                               aaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTg                          <  1:6632248/93‑1 (MQ=255)
                                                                aTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTgg                         <  1:9676691/93‑1 (MQ=255)
                                                                 tCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGc                        <  2:9547187/93‑1 (MQ=255)
                                                                   tGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCgg                      <  1:2289718/93‑1 (MQ=255)
                                                                    gATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGa                     <  1:8602176/93‑1 (MQ=255)
                                                                     aTCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGAt                    <  2:8577143/93‑1 (MQ=255)
                                                                     aTCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGAt                    >  1:6002273/1‑93 (MQ=255)
                                                                       cgccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGa                  <  2:583823/93‑1 (MQ=255)
                                                                        gccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGaa                 <  2:8252454/93‑1 (MQ=255)
                                                                         gcgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAc                <  1:3186695/92‑1 (MQ=255)
                                                                           gcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAGAGTTACCCGACGTATCTGGCGGATGAACCt              <  2:7880352/93‑1 (MQ=255)
                                                                             aCTGGAGGATAACATCGATCCCTGTATTACACCGGCTGGCTCTTGCACCTATGCTGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGa            <  2:5813348/93‑1 (MQ=255)
                                                                              cTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGaa           >  2:7128630/1‑93 (MQ=255)
                                                                              cTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGaa           <  2:7354310/93‑1 (MQ=255)
                                                                              cTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGaa           >  1:3088722/1‑93 (MQ=255)
                                                                               tGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAAt          <  1:2348146/93‑1 (MQ=255)
                                                                                ggaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGTATg         >  1:3733225/1‑93 (MQ=255)
                                                                                ggaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATg         >  2:9168617/1‑93 (MQ=255)
                                                                                  aggaTACCAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTACAAAGTTACCCGACGTATCTGGCGGATGAAACTGAATggg       >  1:7009251/1‑93 (MQ=255)
                                                                                   ggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGc      >  1:4335057/1‑93 (MQ=255)
                                                                                   ggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAACTGAATGGGc      >  2:4056097/1‑93 (MQ=255)
                                                                                    gaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCa     <  2:2162901/93‑1 (MQ=255)
                                                                                    gaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCa     >  1:2166361/1‑93 (MQ=255)
                                                                                    gaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCa     <  2:3793913/93‑1 (MQ=255)
                                                                                    gaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCa     >  1:5810934/1‑93 (MQ=255)
                                                                                     aTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCAt    >  2:610100/1‑93 (MQ=255)
                                                                                      tAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATc   >  1:4467967/1‑93 (MQ=255)
                                                                                      tAACAACGATCCCATTATTTCACCGGCTGGCATGTGCCGCTATGCCGTAAAAAGTTCCCCGACGTACCTGGCGAAGGACCCTGAATGGGCATc   <  2:6809465/93‑1 (MQ=255)
                                                                                       aacaacGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCa  <  2:8740312/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCA  >  NC_000913/320949‑321128

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