Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,908 G→T 17.4% intergenic (+15/+20) mog → / ← satP molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,9080GT17.4% 432.2 / 71.7 174intergenic (+15/+20)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base G (80/62);  new base T (19/12);  total (100/74)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 176 total aligned reads displayed.

GCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCG  >  NC_000913/9809‑10006
                                                                                                   |                                                                                                  
gCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAg                                                                                                    >  2:2401691/1‑100 (MQ=255)
 cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGc                                                                                                   >  2:1549098/1‑100 (MQ=255)
 cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATACCGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCGCCCGAg                                                                                                    >  1:2107921/1‑99 (MQ=255)
  tGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCg                                                                                                  >  2:1339118/1‑100 (MQ=255)
  tGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCg                                                                                                  >  2:3519114/1‑100 (MQ=255)
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCGCGAGCgg                                                                                                 >  1:2388362/1‑100 (MQ=255)
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                 >  2:1093272/1‑100 (MQ=255)
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                 >  2:2612655/1‑100 (MQ=255)
     aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                               >  2:743873/1‑100 (MQ=255)
     aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                               >  2:3642587/1‑100 (MQ=255)
       ggcccATACGTAGAAACGGCACCGGAAGTGGTTGCAGCAGTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggggg                                                                                             <  1:956289/97‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                           >  2:3288033/1‑100 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                          <  2:3123650/100‑1 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                          >  2:2170183/1‑100 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttt                                                                                         <  2:1454763/100‑3 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                         >  2:2967287/1‑100 (MQ=255)
            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtctc                                                                                        <  2:1030557/100‑1 (MQ=255)
             tACCTTGAAACGGCACCGGAAGCGGTTGCAGCATTCAGACCGAAAAGTTCAATACGCGGCGTTAGCGAATAAAAAAATCCCCCCGATCGGGGGGAttttt                                                                                       <  1:532382/100‑5 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGTCGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                      <  2:1878648/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                      <  2:2983796/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaa                                                                                      <  2:1292238/100‑1 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                       <  2:2106332/99‑1 (MQ=255)
              aCGTTGAAACGGCACCGGAAGGGGTTGCAGCATTAAGACTGAGGAGTGCCAGACGCGACGGTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                      <  1:3410357/100‑6 (MQ=255)
                gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAACTCCCCCCGCTCGGGGGGAttttttta                                                                                    <  2:2321185/100‑8 (MQ=38)
                gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                    <  1:1758295/100‑1 (MQ=255)
                gTTGAAAAGCCACCGGAAGTCGTTTCATCATTCAGACCAAAGAGTTCAATACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCTGGGGGATCTCaaaa                                                                                    <  1:15589/100‑1 (MQ=255)
                  tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttttatt                                                                                  <  2:2007279/100‑10 (MQ=38)
                   gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttttatt                                                                                  <  1:2198359/99‑10 (MQ=38)
                    aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAt                                                                                >  1:1429887/1‑100 (MQ=255)
                     aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCATAACAAtt                                                                               <  2:1400125/100‑1 (MQ=255)
                         gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                           <  1:2459702/100‑1 (MQ=255)
                            ccGGAAGTGGTTGCAGCCTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                        >  2:1782880/1‑100 (MQ=255)
                             cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                       >  2:1558334/1‑100 (MQ=255)
                                 aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                   <  2:3794912/100‑1 (MQ=255)
                                 aGTGGTTGCAGCATCTAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                   <  1:2597718/100‑1 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  >  1:1405972/1‑100 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  <  2:2502225/100‑1 (MQ=255)
                                    ggTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCaa                                                                >  1:382812/1‑100 (MQ=255)
                                     gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                               <  2:3465460/100‑1 (MQ=255)
                                     gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                               >  2:824448/1‑100 (MQ=255)
                                      ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                              >  2:320140/1‑100 (MQ=255)
                                       tGCAGCATTGAGAACGCAGAGTGCAAGACGCCAGGTTCGCAAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGGttaccccaccg                                                             >  1:1683282/1‑91 (MQ=255)
                                        gcagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                            >  1:2854528/1‑100 (MQ=255)
                                        gcagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGAAGCTCAAAACAATTAGTGGGATACACCAAt                                                               >  1:410400/1‑97 (MQ=255)
                                          agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                          <  2:1920387/100‑1 (MQ=255)
                                            caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAACCAATTAGTGGGCTTCACCAATCGGCAGa                                                        >  1:3158968/1‑100 (MQ=255)
                                                cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                    >  1:1644303/1‑100 (MQ=255)
                                                cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACg                                                     >  2:1703750/1‑99 (MQ=255)
                                                 ngACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGt                                                   >  2:1232740/2‑100 (MQ=255)
                                                  gagcGAGGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                  <  1:3136552/97‑1 (MQ=255)
                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                  >  2:3440318/1‑100 (MQ=255)
                                                   aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                  >  1:1333793/1‑99 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                <  2:109877/100‑1 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                >  2:2173016/1‑100 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                >  1:3499002/1‑100 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                              >  2:462270/1‑100 (MQ=255)
                                                        agagTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCa                                            <  1:2149792/100‑1 (MQ=255)
                                                          agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCaaa                                          >  1:2739049/1‑100 (MQ=255)
                                                          agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCaaa                                          >  2:3067041/1‑100 (MQ=255)
                                                              cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                      <  2:1429895/100‑1 (MQ=255)
                                                              cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                      >  2:1779188/1‑100 (MQ=255)
                                                              cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                      >  2:849179/1‑100 (MQ=255)
                                                              aaaGACGCGACGTTAGCGAATAAAAAAATACCCCCGAGCGGGGGGATCTCAAAACAATGAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                         >  2:934317/2‑97 (MQ=255)
                                                                 gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                   <  2:3353044/100‑1 (MQ=255)
                                                                   cgcgACGTTAGAGAAGAAAAAAATCCCCCCGCGCGGGGGGATGTTAAAAGAATTAGTGGGTTTCACCAATCGGCCGAACGGTGCGACCAAAGTGCTCg                                   >  1:913729/1‑98 (MQ=255)
                                                                     cgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                               >  2:2480681/1‑100 (MQ=255)
                                                                       aCGTTAGCGAATAAAAAAATCCCCCCGCGCGGGGGGAGCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                             <  1:2666863/100‑1 (MQ=255)
                                                                         gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                           >  1:1173865/1‑100 (MQ=255)
                                                                         gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                           <  1:1726563/100‑1 (MQ=255)
                                                                           tAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                         <  1:1416431/100‑1 (MQ=255)
                                                                           tAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                         >  1:769528/1‑100 (MQ=255)
                                                                            ngCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                        >  1:3007375/2‑100 (MQ=255)
                                                                               gAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                     <  1:1519959/100‑1 (MQ=255)
                                                                                 gttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                   >  2:2770657/9‑100 (MQ=38)
                                                                                  ttttgagaTCCCCCCGCTCGGGGGGCTCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                  >  1:651492/8‑100 (MQ=255)
                                                                                  ttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                  >  2:2882581/8‑100 (MQ=38)
                                                                                  tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                  >  1:2530483/1‑100 (MQ=255)
                                                                                   tttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                    >  1:2057058/7‑97 (MQ=39)
                                                                                   aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  2:3435792/1‑100 (MQ=255)
                                                                                    ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  2:2630498/6‑99 (MQ=39)
                                                                                    ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCAGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  2:3193853/6‑100 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  1:2651081/1‑100 (MQ=255)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGCACAATGCGACCAAACTGCTCGTTCAGTACTTCACCCCTCGc               >  1:3641990/5‑100 (MQ=255)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc               >  1:192966/5‑100 (MQ=255)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTACGACCAAACTGCTCGTTCAGTACTTCACCCAGCg                >  1:3098296/5‑99 (MQ=255)
                                                                                     aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc               <  1:3776738/100‑1 (MQ=255)
                                                                                     aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc               <  2:3450699/100‑1 (MQ=255)
                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc              >  2:586478/4‑100 (MQ=255)
                                                                                      aataTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGATCGTTCAGTACTTCACCCATCGcc              >  1:1631755/4‑100 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  2:1854346/2‑100 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  2:2369380/2‑100 (MQ=255)
                                                                                          tCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCaga           >  2:2450180/1‑99 (MQ=255)
                                                                                          tCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagat          >  2:1915158/1‑100 (MQ=255)
                                                                                           ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata         <  1:2289945/100‑1 (MQ=255)
                                                                                           ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            <  1:1090846/97‑1 (MQ=255)
                                                                                            cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatag        >  2:2400988/1‑100 (MQ=255)
                                                                                              cccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatagat      >  2:1401847/1‑100 (MQ=255)
                                                                                               ccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt     >  1:1766450/1‑100 (MQ=255)
                                                                                                 gAGCGGGGGGGTCTAAAGCCAATGAGAAGGTTTCCCCACTCTGCAGCACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                <  1:2308448/87‑1 (MQ=255)
                                                                                                  aGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACAAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcg  >  2:1603715/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCG  >  NC_000913/9809‑10006

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