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.G98.9% 596.2 / ‑1.6 185coding (374/513 nt)moaBmolybdopterin biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (95/88);  total (97/88)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738
                                                                                             |                                                                                        
gttCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGc                                                                                           >  1:3487659/1‑93 (MQ=255)
gttCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGc                                                                                           <  1:2685743/93‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgg                                                                                        <  1:1157590/93‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgg                                                                                        <  1:6015655/93‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgg                                                                                        <  2:2703694/93‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgg                                                                                        <  1:5944357/93‑1 (MQ=255)
  tCGACCGTGAAGTTGAAGGTTTTGGGGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTGGTGCGGTAGCgg                                                                                        >  2:1924889/1‑93 (MQ=255)
   cGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggg                                                                                      >  2:5181600/1‑93 (MQ=255)
    gACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCgggg                                                                                     <  2:6552959/93‑2 (MQ=255)
     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggc                                                                                    <  1:1746312/93‑3 (MQ=255)
     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggc                                                                                    >  2:6076156/1‑91 (MQ=255)
      ccGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggcg                                                                                   <  1:1610679/93‑4 (MQ=255)
       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggcgt                                                                                  >  2:635123/1‑89 (MQ=255)
       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCggggcgt                                                                                  >  2:665146/1‑89 (MQ=255)
         tGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCg                                                                                 <  2:4502907/93‑1 (MQ=255)
          gAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGc                                                                                >  2:5726297/1‑93 (MQ=255)
            aGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCca                                                                              >  1:2838442/1‑93 (MQ=255)
            aGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCca                                                                              <  1:1175369/93‑1 (MQ=255)
              ttGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaac                                                                            <  2:3740546/93‑1 (MQ=255)
              ttGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaac                                                                            >  2:3027376/1‑93 (MQ=255)
               tGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaaca                                                                           <  1:369596/93‑1 (MQ=255)
                 aaGGTTTTGGTGAAGTGTTCCGGATGGTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCCTGCGGTAGCGGCGCGTCGCccacata                                                                         >  1:4751531/1‑91 (MQ=255)
                  aGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACaaaa                                                                        >  2:3077383/1‑93 (MQ=255)
                    gTNTTGGTNANGTGTTCCGGATGTTGTCNNTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACg                                                                      <  1:4380922/93‑1 (MQ=255)
                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGc                                                                     <  1:4603999/93‑1 (MQ=255)
                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGc                                                                     <  1:1867911/93‑1 (MQ=255)
                        tGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGa                                                                  <  2:1799409/93‑1 (MQ=255)
                        tGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGa                                                                  <  2:4264645/93‑1 (MQ=255)
                         ggTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAt                                                                 >  2:2383226/1‑93 (MQ=255)
                         ggTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTCAt                                                                 <  2:4734289/93‑1 (MQ=255)
                          gTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtt                                                                <  1:5344897/93‑1 (MQ=255)
                          gTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtt                                                                <  1:3078867/93‑1 (MQ=255)
                           tGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAttt                                                               <  1:4902143/93‑1 (MQ=255)
                             aaGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTc                                                             >  2:5052340/1‑93 (MQ=255)
                               gtgtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGc                                                           >  2:3758439/1‑93 (MQ=255)
                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCGACAAAACGCTGATTTTCGCCa                                                         >  1:3368004/1‑93 (MQ=255)
                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCa                                                         <  2:4139927/93‑1 (MQ=255)
                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        <  2:2068096/93‑1 (MQ=255)
                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        >  1:6706511/1‑93 (MQ=255)
                                    ccGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGc                                                      <  2:5664176/93‑1 (MQ=255)
                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    >  2:2195645/1‑93 (MQ=255)
                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    <  1:6402817/93‑1 (MQ=255)
                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:1450127/1‑93 (MQ=255)
                                        aTGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCggg                                                  >  1:106266/1‑93 (MQ=255)
                                        aTGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGATGACTTTCGCCATGCCggg                                                  >  1:4523215/1‑93 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCCACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  2:3507525/1‑93 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 <  2:1017802/93‑1 (MQ=255)
                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 <  2:295891/93‑1 (MQ=255)
                                          gttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGtt                                                <  1:5612272/93‑1 (MQ=255)
                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               >  1:2976296/1‑93 (MQ=255)
                                            tgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCg                                              <  2:6176605/93‑1 (MQ=255)
                                            tgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCg                                              >  1:730769/1‑93 (MQ=255)
                                             gtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGa                                             >  1:2240460/1‑93 (MQ=255)
                                              tCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGAc                                            <  2:6285651/93‑1 (MQ=255)
                                                  ttGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCaaa                                        <  1:3507320/93‑1 (MQ=255)
                                                   tGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAg                                       <  2:5128295/93‑1 (MQ=255)
                                                     aaGAGATTGGCACTTCGACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCg                                     <  1:5212940/93‑1 (MQ=255)
                                                      agagATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCACCAAAACGCTGATTTTCGCCATCCCGGGTTCGACCAAAGCgg                                    >  1:1863131/1‑92 (MQ=255)
                                                      agagATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGt                                    >  2:2634802/1‑93 (MQ=255)
                                                      agagATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGt                                    >  2:338242/1‑93 (MQ=255)
                                                       gagaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTg                                   >  2:6700257/1‑93 (MQ=255)
                                                        agaTTGGCACTTCCCCGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGc                                  <  1:1311375/93‑1 (MQ=255)
                                                        agaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGc                                  <  2:46141/93‑1 (MQ=255)
                                                        agaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGc                                  >  1:4947030/1‑93 (MQ=255)
                                                        agaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAACGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGc                                  >  2:6218968/1‑93 (MQ=255)
                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 >  2:4148700/1‑93 (MQ=255)
                                                           ttGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGt                               >  1:2026275/1‑93 (MQ=255)
                                                           ttGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGt                               <  1:2385221/93‑1 (MQ=255)
                                                              gCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTAcc                            <  1:753512/93‑1 (MQ=255)
                                                               cACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCg                           <  1:3649364/93‑1 (MQ=255)
                                                               cACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCg                           >  1:3929698/1‑93 (MQ=255)
                                                                aCTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGc                          >  1:5504685/1‑93 (MQ=255)
                                                                aCTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGc                          >  1:6498507/1‑93 (MQ=255)
                                                                aCTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGc                          >  2:5271897/1‑93 (MQ=255)
                                                                aCTTCCACGTTCCAATCTCGTGCCGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGc                          <  2:5587273/93‑1 (MQ=255)
                                                                 cTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCa                         >  1:2050204/1‑93 (MQ=255)
                                                                 cTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCa                         >  1:4126263/1‑93 (MQ=255)
                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        <  2:2824823/93‑1 (MQ=255)
                                                                   tCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATg                       <  1:5875650/93‑1 (MQ=255)
                                                                    ccACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATgg                      >  2:1490816/1‑93 (MQ=255)
                                                                     cACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATggg                     <  1:2188804/93‑1 (MQ=255)
                                                                      aCGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGa                    >  1:3957367/1‑93 (MQ=255)
                                                                      aCGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGa                    >  2:1145752/1‑93 (MQ=255)
                                                                        gTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaa                  <  2:1382672/93‑1 (MQ=255)
                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATgggcaca                 >  2:1453654/1‑89 (MQ=255)
                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  2:3258285/1‑93 (MQ=255)
                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  2:388546/1‑93 (MQ=255)
                                                                          tGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGACAat                >  1:626597/1‑93 (MQ=255)
                                                                          tGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAat                <  1:576171/93‑1 (MQ=255)
                                                                            cAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAatat              <  2:4079650/93‑1 (MQ=255)
                                                                             aaTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatc             >  1:2143193/1‑93 (MQ=255)
                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           >  1:5220789/1‑93 (MQ=255)
                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           <  1:6193074/93‑1 (MQ=255)
                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           <  2:60262/93‑1 (MQ=255)
                                                                                   gTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgcg       <  1:4175478/93‑1 (MQ=255)
                                                                                     gCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgcc     <  1:4864968/93‑1 (MQ=255)
                                                                                      cGGTAGCGGGGCGTCGCCAACACAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccg    >  1:1656298/1‑93 (MQ=255)
                                                                                       ggTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccgc   >  1:3334540/1‑93 (MQ=255)
                                                                                       ggTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccgc   >  2:2541879/1‑93 (MQ=255)
                                                                                        gTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCa  >  2:358946/1‑93 (MQ=255)
                                                                                             |                                                                                        
GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738

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