Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 817,654 +G coding (377/513 nt) moaB → molybdopterin biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913817,6511.G97.5% 365.4 / 2.1 120coding (374/513 nt)moaBmolybdopterin biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (69/48);  total (70/50)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738
                                                                                              |                                                                                        
tttCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGc                                                                                           <  1:2718625/92‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgg                                                                                       <  2:1498362/93‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgg                                                                                       <  1:3343223/93‑1 (MQ=255)
     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggc                                                                                    <  2:4150906/93‑3 (MQ=255)
     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggc                                                                                    <  1:1730336/93‑3 (MQ=255)
       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCCGGGCGt                                                                                   <  2:3901354/93‑1 (MQ=255)
       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggcgt                                                                                  >  2:2695629/1‑89 (MQ=255)
       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggcgt                                                                                  >  1:1706960/1‑89 (MQ=255)
       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggcgt                                                                                  <  2:4131024/93‑5 (MQ=255)
          gAAGTTGAAGGTTTTGGGGAAGTGTTCCGTATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGGGCGGAAGCGGGGCGTCGc                                                                                >  1:3579014/1‑93 (MQ=255)
           aaGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGcc                                                                               >  2:3041244/1‑93 (MQ=255)
           aaGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGcc                                                                               <  2:462449/93‑1 (MQ=255)
            aGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCca                                                                              <  1:526784/93‑1 (MQ=255)
             gTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaa                                                                             >  2:2649132/1‑93 (MQ=255)
               tGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaata                                                                           >  1:3926104/1‑91 (MQ=255)
               tGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaaca                                                                           <  2:379558/93‑1 (MQ=255)
               tGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaaca                                                                           <  1:4533995/93‑1 (MQ=255)
                gAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaacaa                                                                          <  1:1662282/93‑1 (MQ=255)
                gAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaacaa                                                                          >  2:2308083/1‑93 (MQ=255)
                gAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaacaa                                                                          >  2:3646755/1‑93 (MQ=255)
                gAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaacaa                                                                          <  1:1070886/93‑1 (MQ=255)
                 agggTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGGTGCAATGTCGTGCGGTAGCGGGGCGTCGCCAACaaa                                                                         >  2:1301289/3‑93 (MQ=255)
                 aaGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACaaa                                                                         >  1:1829835/1‑93 (MQ=255)
                  aGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACaaaa                                                                        <  2:514761/93‑1 (MQ=255)
                    gTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACg                                                                      >  1:1364374/1‑93 (MQ=255)
                    gTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACg                                                                      >  2:3575800/1‑93 (MQ=255)
                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGc                                                                     >  2:2250158/1‑93 (MQ=255)
                        tGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGa                                                                  <  1:3961139/93‑1 (MQ=255)
                        tGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGa                                                                  <  2:4241587/93‑1 (MQ=255)
                         ggTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGGCAACAAAACGCTGAt                                                                 >  1:90631/1‑93 (MQ=255)
                         ggTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAt                                                                 <  2:2464140/93‑1 (MQ=255)
                          gTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCGACAAAACGCTGAtt                                                                <  1:2408447/93‑1 (MQ=255)
                          gTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtt                                                                <  2:3750266/93‑1 (MQ=255)
                           tGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCGACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAttt                                                               >  1:3305851/1‑93 (MQ=255)
                            gAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtttt                                                              <  2:886606/93‑1 (MQ=255)
                              aGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCg                                                            >  1:375848/1‑93 (MQ=255)
                               gtgtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGc                                                           <  1:195279/93‑1 (MQ=255)
                               gtgtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGc                                                           >  1:883220/1‑93 (MQ=255)
                                tgtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGcc                                                          <  1:101331/93‑1 (MQ=255)
                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGCGGCGTCGCCAACAAAACGCTGATTTTCGcct                                                         >  1:2734594/1‑92 (MQ=255)
                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        <  1:2916224/93‑1 (MQ=255)
                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        <  1:2546361/93‑1 (MQ=255)
                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        >  2:953218/1‑93 (MQ=255)
                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    >  1:37241/1‑93 (MQ=255)
                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    >  2:2469177/1‑93 (MQ=255)
                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    >  2:1045088/1‑93 (MQ=255)
                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   <  2:3866252/93‑1 (MQ=255)
                                        aTGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCggg                                                  <  1:1910849/93‑1 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  2:2883243/1‑93 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 <  1:987291/93‑1 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  1:2261821/1‑93 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGATGATTTTCGCGATGCCGGGt                                                 >  1:2051805/1‑93 (MQ=255)
                                          tttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGtt                                                <  2:540203/92‑1 (MQ=255)
                                          gttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGtt                                                >  1:4478335/1‑93 (MQ=255)
                                          gttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAACGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGtt                                                <  2:4566522/93‑1 (MQ=255)
                                            tgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGGAGCGGGGCGGCGACAACAAAACGCTGATTTTTGCCATGCCGGGCTCg                                              >  2:2661047/1‑93 (MQ=255)
                                                gTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCa                                          >  2:3028806/1‑93 (MQ=255)
                                                 tttGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCaa                                         >  2:3619484/1‑93 (MQ=255)
                                                  ttGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGGAGCGGGGCGGCGCCAACAAAATGCTGATTTTCGCCATGCCGGGTTCGAccgag                                        >  1:2102766/1‑90 (MQ=255)
                                                  tagAAGAGATTGGCTCTTCCATGTTGCAATCTCGTGCGGTAGCTGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCaaa                                        <  2:2561337/91‑1 (MQ=255)
                                                   tGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAg                                       <  1:3062868/93‑1 (MQ=255)
                                                    gAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGc                                      <  2:921803/93‑1 (MQ=255)
                                                    gAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGc                                      >  2:121340/1‑93 (MQ=255)
                                                    gAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGc                                      >  2:2770680/1‑93 (MQ=255)
                                                    gAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGGCGCGAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGc                                      >  1:1777108/1‑93 (MQ=255)
                                                        ggaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGc                                  <  2:1376267/92‑1 (MQ=255)
                                                        agaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGc                                  >  2:1474575/1‑93 (MQ=255)
                                                          aTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCg                                <  2:4409925/93‑1 (MQ=255)
                                                               cACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTAccc                           >  1:1455196/1‑92 (MQ=255)
                                                                ccTTCCACGCTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTCCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGc                          <  1:3195299/92‑1 (MQ=255)
                                                                 cTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCa                         >  1:2488053/1‑93 (MQ=255)
                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        >  1:279926/1‑93 (MQ=255)
                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        >  2:2460689/1‑93 (MQ=255)
                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATg                       <  2:4622633/93‑1 (MQ=255)
                                                                     cACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATggg                     >  1:2750776/1‑93 (MQ=255)
                                                                      aCGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGa                    >  1:460347/1‑93 (MQ=255)
                                                                      aCGTTGCAATCTCGTGCGGAAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGa                    >  2:4579510/1‑93 (MQ=255)
                                                                        gTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaa                  >  1:2704577/1‑93 (MQ=255)
                                                                        gTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaa                  >  1:4333808/1‑93 (MQ=255)
                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGCACCGCATGGGaaaa                 >  2:30670/1‑93 (MQ=255)
                                                                          tGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAat                >  1:3411700/1‑93 (MQ=255)
                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               <  1:1972290/93‑1 (MQ=255)
                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               <  1:1030992/93‑1 (MQ=255)
                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               <  1:4076499/93‑1 (MQ=255)
                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               >  2:2525549/1‑93 (MQ=255)
                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               <  2:1825104/93‑1 (MQ=255)
                                                                            cAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAatat              >  1:1682796/1‑93 (MQ=255)
                                                                            cAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAatat              <  1:1555401/93‑1 (MQ=255)
                                                                             aaTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAatata             >  2:2085886/1‑92 (MQ=255)
                                                                             aaTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatc             <  2:17847/93‑1 (MQ=255)
                                                                             aaTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatc             >  1:16060/1‑93 (MQ=255)
                                                                              aTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            >  2:4568077/1‑93 (MQ=255)
                                                                              aTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            <  1:136068/93‑1 (MQ=255)
                                                                              aTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            >  1:3184037/1‑93 (MQ=255)
                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           >  1:3489820/1‑93 (MQ=255)
                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           <  2:4010143/93‑1 (MQ=255)
                                                                                ctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcatc          <  2:1618673/93‑1 (MQ=255)
                                                                                 tcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcg         >  2:2103435/1‑93 (MQ=255)
                                                                                 tcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcg         >  2:1850659/1‑93 (MQ=255)
                                                                                   gTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgcg       >  1:292906/1‑93 (MQ=255)
                                                                                              |                                                                                        
GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738

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