Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,907 A→C 17.5% intergenic (+14/+21) 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,9070AC17.5% 395.1 / 82.6 175intergenic (+14/+21)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base A (81/63);  new base C (20/11);  total (101/74)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

Only 100 of 180 total aligned reads displayed.

TGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCG  >  NC_000913/9808‑10006
                                                                                                   |                                                                                                   
tGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGa                                                                                                     >  2:3661887/1‑100 (MQ=255)
 gCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAg                                                                                                    >  2:2401691/1‑100 (MQ=255)
  cTGGAAGGGCCATACGTTGAAACGGCACCTGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGc                                                                                                   <  1:2299671/100‑1 (MQ=255)
  cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGc                                                                                                   <  1:2173870/100‑1 (MQ=255)
  cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATACCGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCGCCCGAg                                                                                                    >  1:2107921/1‑99 (MQ=255)
   tGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCg                                                                                                  >  2:3519114/1‑100 (MQ=255)
   tGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCg                                                                                                  >  2:1416439/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCGCGAGCgg                                                                                                 >  1:2388362/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                 >  2:2612655/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                 >  2:1093272/1‑100 (MQ=255)
      aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                               >  2:743873/1‑100 (MQ=255)
      aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCAAAGAGCGAAAGTCGCGACGTTAGCGAAGAAAAAAATCCCCCAGAGCggg                                                                                                >  1:552982/1‑99 (MQ=255)
         ggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGa                                                                                            <  2:1907084/100‑1 (MQ=255)
           ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                          <  2:3123650/100‑1 (MQ=255)
            cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                         >  2:2967287/1‑100 (MQ=255)
               aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                       <  2:2106332/99‑1 (MQ=255)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAACTCCCCCCGCTCGGGGGGAttttttta                                                                                    <  2:2321185/100‑8 (MQ=38)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                    <  1:1758295/100‑1 (MQ=255)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                    >  1:3098332/1‑100 (MQ=255)
                 gTTGAAAAGCCACCGGAAGTCGTTTCATCATTCAGACCAAAGAGTTCAATACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCTGGGGGATCTCaaaa                                                                                    <  1:15589/100‑1 (MQ=255)
                   tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttttatt                                                                                  <  2:2007279/100‑10 (MQ=38)
                      aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCATAACAAtt                                                                               <  2:1400125/100‑1 (MQ=255)
                         ggCACCGGAAGTGGTTGCAGCATTCAAACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATCAGt                                                                            <  1:2260611/100‑1 (MQ=255)
                          gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                           <  1:2459702/100‑1 (MQ=255)
                           cACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                          >  1:3377723/1‑100 (MQ=255)
                            aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAATATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                         <  1:1694731/100‑1 (MQ=255)
                             ccGGAAGTGGTTGCAGCCTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                        >  2:1782880/1‑100 (MQ=255)
                             ccGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                        >  2:2149802/1‑100 (MQ=255)
                              cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                       >  2:1558334/1‑100 (MQ=255)
                              cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                       <  2:3447583/100‑1 (MQ=255)
                                  aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                   <  2:300098/100‑1 (MQ=255)
                                  aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                   <  2:3794912/100‑1 (MQ=255)
                                  aGTGGTTGCAGCATCTAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                   <  1:2597718/100‑1 (MQ=255)
                                   gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  <  2:2502225/100‑1 (MQ=255)
                                   gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  <  2:1436714/100‑1 (MQ=255)
                                   gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  >  1:1405972/1‑100 (MQ=255)
                                     ggTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCaa                                                                >  1:382812/1‑100 (MQ=255)
                                      gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                               >  2:824448/1‑100 (MQ=255)
                                       ntGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                              >  2:1367680/2‑100 (MQ=255)
                                        tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                             <  1:2091713/100‑1 (MQ=255)
                                         gcagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                            >  1:513303/1‑100 (MQ=255)
                                           agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                          <  2:1920387/100‑1 (MQ=255)
                                             caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAACAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGa                                                        <  1:1652609/100‑1 (MQ=255)
                                               ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                      >  2:3247706/1‑100 (MQ=255)
                                                    aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                  >  1:1333793/1‑99 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                >  2:2173016/1‑100 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                >  1:3499002/1‑100 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                <  2:1137462/100‑1 (MQ=255)
                                                       gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                              >  2:2381192/1‑100 (MQ=255)
                                                       gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                              >  2:462270/1‑100 (MQ=255)
                                                        aaGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAAAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                              >  1:1439643/1‑99 (MQ=255)
                                                           agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCaaa                                          >  2:3067041/1‑100 (MQ=255)
                                                             tGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACt                                        <  1:994686/100‑1 (MQ=255)
                                                             tGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACt                                        >  2:1247361/1‑100 (MQ=255)
                                                               cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                      >  1:1847547/1‑100 (MQ=255)
                                                               cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                      >  2:849179/1‑100 (MQ=255)
                                                               aaaGACGCGACGTTAGCGAATAAAAAAATACCCCCGAGCGGGGGGATCTCAAAACAATGAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                         >  2:934317/2‑97 (MQ=255)
                                                                 aGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                    >  2:1726572/1‑100 (MQ=255)
                                                                  gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGTGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                   <  1:934310/100‑1 (MQ=255)
                                                                  gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                   <  2:3353044/100‑1 (MQ=255)
                                                                  gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                   >  1:2721634/1‑100 (MQ=255)
                                                                   aCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                  >  2:1327007/1‑100 (MQ=255)
                                                                     gtcaCGTTAGCGAATAAAAAGATCCCCCCGATCGAGGGGATCTCAAAACAATGAGTGGGATCCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                <  1:2891423/97‑1 (MQ=255)
                                                                      cgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                               >  2:2480681/1‑100 (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)
                                                                             aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCATAACAATTAGTGGGATTCACCAAGCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                        >  1:2047888/1‑100 (MQ=255)
                                                                              gCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCa                       >  2:3488185/1‑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)
                                                                                   ttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                  >  2:1032787/8‑100 (MQ=38)
                                                                                   tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                  >  1:2530483/1‑100 (MQ=255)
                                                                                    aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  2:3435792/1‑100 (MQ=255)
                                                                                     ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  2:2630498/6‑99 (MQ=39)
                                                                                     ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCAGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  2:3193853/6‑100 (MQ=255)
                                                                                     aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  2:3349278/1‑100 (MQ=255)
                                                                                     aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  1:2651081/1‑100 (MQ=255)
                                                                                     aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATAAGTGGGATTCACCAATGGGAAGAACGGTGCGACCAAACTGGTCGTTGAGTACTTCACCCATCg                >  1:3142643/1‑100 (MQ=255)
                                                                                      tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGCACAATGCGACCAAACTGCTCGTTCAGTACTTCACCCCTCGc               >  1:3641990/5‑100 (MQ=255)
                                                                                      tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTACGACCAAACTGCTCGTTCAGTACTTCACCCAGCg                >  1:3098296/5‑99 (MQ=255)
                                                                                      aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc               <  1:3776738/100‑1 (MQ=255)
                                                                                       gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc              >  2:586478/4‑100 (MQ=255)
                                                                                       aataTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGATCGTTCAGTACTTCACCCATCGcc              >  1:1631755/4‑100 (MQ=255)
                                                                                       aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc              <  2:3460959/100‑1 (MQ=255)
                                                                                        aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa             <  1:2646814/100‑1 (MQ=255)
                                                                                        aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa             <  2:1527660/100‑1 (MQ=255)
                                                                                         gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  2:2369380/2‑100 (MQ=255)
                                                                                         gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  2:1854346/2‑100 (MQ=255)
                                                                                         gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  2:1561401/2‑100 (MQ=255)
                                                                                           tCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCaga           >  2:2450180/1‑99 (MQ=255)
                                                                                            ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagat          <  1:3519093/99‑1 (MQ=255)
                                                                                            ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            <  1:1090846/97‑1 (MQ=255)
                                                                                             cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagcgca        >  1:2435810/1‑96 (MQ=255)
                                                                                             cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatag        >  1:2745086/1‑100 (MQ=255)
                                                                                              ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagataga       >  2:1870989/1‑100 (MQ=255)
                                                                                                 cGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTg    <  2:1682084/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCG  >  NC_000913/9808‑10006

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