Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,907 A→C 15.8% 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,9070AC15.8% 289.8 / 47.8 123intergenic (+14/+21)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base A (40/63);  new base C (10/10);  total (50/73)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 125 total aligned reads displayed.

TGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTG  >  NC_000913/9808‑10004
                                                                                                   |                                                                                                 
tGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGa                                                                                                   >  2:1251232/1‑100 (MQ=255)
  cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGGGACGTTAGCGAGTAAAAAAATCCCCCCGAGc                                                                                                 >  2:131522/1‑100 (MQ=255)
  cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGc                                                                                                 >  1:40649/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                               >  1:2045295/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                               >  1:2189408/1‑100 (MQ=255)
      aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                              >  1:1019744/1‑99 (MQ=255)
       aGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggggg                                                                                            >  2:1845603/1‑100 (MQ=255)
       aGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggggg                                                                                            <  2:1163/100‑1 (MQ=255)
       aGGGCCATACGTTGAAACGGCACCGGAAGTGGCTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggggg                                                                                            <  2:176328/100‑1 (MQ=255)
        gggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCgggggg                                                                                           <  2:952333/100‑1 (MQ=255)
        gggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCgggggg                                                                                           <  1:1445220/100‑1 (MQ=255)
         ggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCTAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggggg                                                                                           <  1:1796743/99‑1 (MQ=255)
         ggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGa                                                                                          <  2:1431678/100‑1 (MQ=255)
          gcgagaCGGTGACACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCACGACGCGACGCAATAAAAAAAAAAAATCCCCCCGAGCGGGGGGa                                                                                          <  1:1906235/94‑1 (MQ=255)
          gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                         >  1:1496886/1‑100 (MQ=255)
           ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGGGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                        <  1:2771911/100‑1 (MQ=255)
            cagaCGTTGAAACGGCACCGGAAGTCGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                        <  1:827748/96‑1 (MQ=255)
            cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttt                                                                                       <  2:2205375/100‑3 (MQ=255)
              tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                     <  2:2358546/100‑1 (MQ=255)
               aggTTGAAACGGCACCGGAAGTGGTTGCATCATTCAGACCTAAGAGCGCAAGACGAGACGTTATCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                    <  1:370017/98‑6 (MQ=255)
               aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaa                                                                                    <  2:1386631/100‑1 (MQ=255)
                cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAGGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttttttt                                                                                   <  1:2389478/100‑7 (MQ=38)
                cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttttttt                                                                                   <  1:1917544/100‑7 (MQ=39)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                  <  2:2678957/100‑1 (MQ=255)
                  ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                 <  1:973742/100‑1 (MQ=255)
                  tcgAAACTTCACTGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCTAGACGCGACCTTAGCTAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaacc                                                                                 <  1:1657953/98‑3 (MQ=255)
                   tGAAACGGCACCGGAAGTGGTTTCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCCAATAAAAATTTCCCCCCGCGCGGGGGGCTCTCAAAACa                                                                                <  1:650603/100‑1 (MQ=255)
                   tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                 <  1:261030/99‑1 (MQ=255)
                    gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACaa                                                                               >  1:1383157/1‑100 (MQ=255)
                     aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACaa                                                                               <  1:280776/99‑1 (MQ=255)
                     aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAt                                                                              <  1:813293/100‑1 (MQ=255)
                     aaaCGGCACCGGAAGTGGTTGCAGCAATCAGACCGAAGAGTGCAAGACGCGACGGTAGCGAATAAAAAAATCCCCCCGCGCGGGGGGTTCTCAAAACAAt                                                                              <  1:2740475/100‑1 (MQ=255)
                       aCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTa                                                                            <  2:326053/100‑1 (MQ=255)
                        cGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                           >  2:112292/1‑100 (MQ=255)
                         ggCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGt                                                                          >  1:1721006/1‑100 (MQ=255)
                          gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                         <  2:2669209/100‑1 (MQ=255)
                             ccGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                       <  1:183771/99‑1 (MQ=255)
                              cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                     <  1:2094843/100‑1 (MQ=255)
                                gAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTc                                                                   >  1:2650873/1‑100 (MQ=255)
                                  aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                 <  2:1507564/100‑1 (MQ=255)
                                   gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAACAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                <  1:1615139/100‑1 (MQ=255)
                                   gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                >  1:1847015/1‑100 (MQ=255)
                                     ggTTCCAGCATTCAGACAGAAGAGCGCAAGACGCGACGTTAGCGCATAAAAAAAACCCCACGCGCGGGGGGAGCGCAAAACAATTAGTGGCATTCACCa                                                               >  1:1594513/1‑99 (MQ=255)
                                       ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                            <  1:1689046/100‑1 (MQ=255)
                                        tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                           <  2:1350162/100‑1 (MQ=255)
                                        tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                           <  2:1350445/100‑1 (MQ=255)
                                        tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                           <  1:2219793/100‑1 (MQ=255)
                                          catcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                         <  1:2566591/97‑1 (MQ=255)
                                          cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATATCCAAACAATTAGTGGGATTCACCAACCGGc                                                         >  1:1026988/1‑100 (MQ=255)
                                            gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCaa                                                       >  2:2729291/1‑99 (MQ=255)
                                             caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGATCGCGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGa                                                      <  1:897233/100‑1 (MQ=255)
                                             caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGa                                                      <  2:2645328/100‑1 (MQ=255)
                                              aTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGaa                                                     >  1:2249271/1‑100 (MQ=255)
                                               ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                    >  2:457342/1‑100 (MQ=255)
                                                 cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                  >  2:1660365/1‑100 (MQ=255)
                                                 cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                  <  2:2535429/100‑1 (MQ=255)
                                                   gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                <  2:2497288/100‑1 (MQ=255)
                                                    aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                               <  1:954510/100‑1 (MQ=255)
                                                    aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                               <  1:1188152/100‑1 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGTGATTCACCAATCGGCAGAACGGTGCg                                              >  1:1583005/1‑100 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGTATTCACCAATCGGCAGAACGGTGCg                                              <  2:1242414/100‑1 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAACTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                              >  1:2096001/1‑100 (MQ=255)
                                                       gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                            <  2:221666/100‑1 (MQ=255)
                                                       gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                            >  1:1383472/1‑100 (MQ=255)
                                                               cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGATCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                    <  1:1954616/100‑1 (MQ=255)
                                                               cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                    >  1:2783764/1‑100 (MQ=255)
                                                                aaGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCt                                   <  1:2367627/100‑1 (MQ=255)
                                                                aaGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCt                                   <  2:1651875/100‑1 (MQ=255)
                                                                 aGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                  >  2:614620/1‑100 (MQ=255)
                                                                 aGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                  >  1:2604456/1‑100 (MQ=255)
                                                                   aCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                >  2:1951691/1‑100 (MQ=255)
                                                                        aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                           <  1:244416/100‑1 (MQ=255)
                                                                        aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                           >  1:587084/1‑100 (MQ=255)
                                                                        aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                           <  2:20485/100‑1 (MQ=255)
                                                                         cGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                          <  1:979682/100‑1 (MQ=255)
                                                                          gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCGCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                         <  1:1022739/100‑1 (MQ=255)
                                                                            tAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                       >  1:520546/1‑100 (MQ=255)
                                                                              gCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCa                     >  2:2769846/1‑100 (MQ=255)
                                                                               cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                    <  2:2237107/100‑1 (MQ=255)
                                                                               cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                    <  2:2763767/100‑1 (MQ=255)
                                                                                 tgttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTAAcc                   >  1:2511248/10‑99 (MQ=37)
                                                                                 aaTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                  <  2:387264/100‑1 (MQ=255)
                                                                                 aaTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                  <  2:2571190/100‑1 (MQ=255)
                                                                                  gttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                 >  2:2573495/9‑100 (MQ=38)
                                                                                  aTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                 >  2:399032/1‑100 (MQ=255)
                                                                                   ttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGGGGGGTTGACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                >  2:2037894/8‑100 (MQ=37)
                                                                                   tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                >  1:1446667/1‑100 (MQ=255)
                                                                                      aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc             <  2:2010752/100‑1 (MQ=255)
                                                                                       gagaTCCCCCCGCTCGGGGGGACCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACAGAGCCACCACACTGCTCGTTCAGTACTTCACCCATCGcc            >  1:1286613/4‑100 (MQ=255)
                                                                                       aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc            <  2:1611266/100‑1 (MQ=255)
                                                                                        agaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa           >  2:821100/3‑100 (MQ=255)
                                                                                        aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa           <  2:469307/100‑1 (MQ=255)
                                                                                         gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag          >  2:725764/2‑100 (MQ=255)
                                                                                         aaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag          <  1:1308617/100‑1 (MQ=255)
                                                                                            ccccccGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata       >  2:102901/1‑100 (MQ=255)
                                                                                            acccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCTCCCATCGCAagat        <  1:1121640/98‑1 (MQ=255)
                                                                                             cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatag      >  2:147339/1‑100 (MQ=255)
                                                                                              ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagataga     >  2:335712/1‑100 (MQ=255)
                                                                                               cccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatagat    >  1:1227947/1‑100 (MQ=255)
                                                                                                 cGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTg  <  1:2229740/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
TGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTG  >  NC_000913/9808‑10004

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: