Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,907 A→C 15.7% 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.7% 504.5 / 84.4 207intergenic (+14/+21)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base A (91/83);  new base C (17/16);  total (108/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 208 total aligned reads displayed.

TGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGC  >  NC_000913/9808‑10005
                                                                                                   |                                                                                                  
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                >  1:589374/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                >  1:3199989/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                >  1:1602795/1‑100 (MQ=255)
    ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                >  2:2846264/1‑100 (MQ=255)
     gAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                               <  2:237235/100‑1 (MQ=255)
         ggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGa                                                                                           >  1:630795/1‑100 (MQ=255)
          gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                          >  1:3660927/1‑100 (MQ=255)
          gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                          >  2:1411515/1‑100 (MQ=255)
           ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                         <  1:3998418/100‑1 (MQ=255)
           ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                         >  2:4221922/1‑100 (MQ=255)
           ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                         >  1:2174720/1‑100 (MQ=255)
            cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                        <  1:2832864/100‑1 (MQ=255)
            cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                        >  2:3483670/1‑100 (MQ=255)
            cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                         >  1:3499613/1‑99 (MQ=255)
             aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                       <  1:2044363/100‑4 (MQ=255)
             aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                       <  2:364717/100‑4 (MQ=255)
             aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtctc                                                                                       <  2:793729/100‑1 (MQ=255)
             aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTAAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                        <  1:3918236/99‑1 (MQ=255)
              tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                      >  2:2228560/1‑100 (MQ=255)
              tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGAAGCGACGTTAACGAATAAAAAAACACCCCCGAGCGGGAGGCGCTCa                                                                                      >  1:3718144/1‑100 (MQ=255)
               aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                     <  1:1075876/100‑6 (MQ=255)
               aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                     <  2:1825215/100‑6 (MQ=255)
               aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGCGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttcttt                                                                                     <  1:325928/100‑6 (MQ=255)
                cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttttttt                                                                                    <  2:4157512/100‑7 (MQ=39)
                cGTTGAAACGGCACCGGAAGTGGTTGCAGCATCCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaa                                                                                    <  1:1336779/100‑1 (MQ=255)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATATAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                   >  1:2590481/1‑100 (MQ=255)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                   >  2:3699248/1‑100 (MQ=255)
                 gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                   >  1:2951247/1‑100 (MQ=255)
                  ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttttat                                                                                  <  2:3844974/100‑9 (MQ=38)
                  ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                  <  2:2445619/100‑1 (MQ=255)
                  ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                  <  2:16282/100‑1 (MQ=255)
                  ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                  <  2:1542262/100‑1 (MQ=255)
                   tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAATAAAATCCCCCCGCTCGGGGGGAtttttttatt                                                                                 <  2:4184729/100‑10 (MQ=38)
                   tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa                                                                                 >  1:3602494/1‑100 (MQ=255)
                   tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa                                                                                 <  2:120318/100‑1 (MQ=255)
                   tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa                                                                                 >  1:2602209/1‑100 (MQ=255)
                      aaCGGCACCGGAAGTGTTTGCAGCAATCAGACTGACTAGTGCAAGACGCGACTTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                     <  1:1858188/93‑6 (MQ=25)
                      aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAtt                                                                              <  2:3062176/100‑1 (MQ=255)
                      aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAtt                                                                              >  1:592670/1‑100 (MQ=255)
                      aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAtt                                                                              <  2:2602223/100‑1 (MQ=255)
                      aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAtt                                                                              <  1:4240632/100‑1 (MQ=255)
                        cGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                            >  1:342895/1‑100 (MQ=255)
                         ggCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGCGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                            >  1:1111152/1‑99 (MQ=255)
                          gCACCGGCAGTGATTGCAGAAGTCATACAGATGATTGCAAGAAGAAAAAATAAAAAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                          <  1:1645167/100‑1 (MQ=255)
                          gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                          <  2:3517276/100‑1 (MQ=255)
                            aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                        <  2:2941165/100‑1 (MQ=255)
                             ccGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTTGCGAATAAAAAAATACCCCCTAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                       >  2:571326/1‑100 (MQ=255)
                              cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                      <  2:1125991/100‑1 (MQ=255)
                                  aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                  >  1:2299075/1‑100 (MQ=255)
                                   gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                 <  2:2891942/100‑1 (MQ=255)
                                     ggTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCaa                                                               >  2:2199614/1‑100 (MQ=255)
                                      gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                              >  2:1154311/1‑100 (MQ=255)
                                          cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                          >  2:2157168/1‑100 (MQ=255)
                                           agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                         >  1:4064048/1‑100 (MQ=255)
                                            gccttCACACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCCCGCGGGGGGATCCCACAACCATTAGTGGCACTCACCAATCGGCAg                                                        >  1:2649055/4‑100 (MQ=255)
                                            gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCCCCAATCgg                                                           >  1:1942185/1‑97 (MQ=255)
                                            gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                        <  2:4054161/100‑1 (MQ=255)
                                            gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                        >  2:2505639/1‑100 (MQ=255)
                                             caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGa                                                       >  2:1059334/1‑100 (MQ=255)
                                               gacAGACCGAAGATGGCAAGACACAACATTAGCGAATAAAAAAATCCCTCCGAACGGGGGTATCTAACAAAAATTAACAGAATTCAACAATCGGCAGaa                                                      >  1:1561013/3‑99 (MQ=255)
                                                 cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                   >  1:850290/1‑100 (MQ=255)
                                                   gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                 <  1:922606/100‑1 (MQ=255)
                                                    aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                                >  2:4240545/1‑100 (MQ=255)
                                                    aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                                <  2:3136096/100‑1 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                               >  1:224242/1‑100 (MQ=255)
                                                     ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGACAGAACGGTGCg                                               >  1:1340123/1‑100 (MQ=255)
                                                       gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                             <  2:187863/100‑1 (MQ=255)
                                                          gagTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAAGGTGc                                                >  1:346913/1‑94 (MQ=255)
                                                              gCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATATCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTg                                      >  1:1353078/1‑100 (MQ=255)
                                                                  gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCACAACAATTACTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                  <  2:874453/100‑1 (MQ=255)
                                                                    cgcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                >  2:3916388/1‑100 (MQ=255)
                                                                     gcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                               <  2:1450045/100‑1 (MQ=255)
                                                                      cgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                              <  2:1580309/100‑1 (MQ=255)
                                                                        aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                            >  2:3739805/1‑100 (MQ=255)
                                                                        aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                            >  2:3237360/1‑100 (MQ=255)
                                                                          gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCACTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                          <  1:1412160/100‑1 (MQ=255)
                                                                          gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                          >  2:1490187/1‑100 (MQ=255)
                                                                             aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                       <  1:3395624/100‑1 (MQ=255)
                                                                                 tgttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                   >  2:1352822/10‑100 (MQ=38)
                                                                                 aaTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                   <  1:1138138/100‑1 (MQ=255)
                                                                                  aTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                  >  1:4400943/1‑100 (MQ=255)
                                                                                   tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                 >  2:1596810/1‑100 (MQ=255)
                                                                                   tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                 >  2:1734028/1‑100 (MQ=255)
                                                                                   tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                 <  2:2564513/100‑1 (MQ=255)
                                                                                    aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                <  2:3329656/100‑1 (MQ=255)
                                                                                     ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg               >  1:4201436/6‑100 (MQ=255)
                                                                                      tgagaTCCCCCCGCTCGGGGGGATCTTAAAACAATCAGTGGGATTCACCAGTCGGCAGAACGGTGCGACCAAACTGCCCGTTCAGTACTTCACCCATCGc              >  1:381590/5‑100 (MQ=255)
                                                                                      ngagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc              >  2:3271477/5‑100 (MQ=255)
                                                                                       gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCCCCCATCg               >  1:2064649/4‑98 (MQ=255)
                                                                                       gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc             >  1:3158601/4‑100 (MQ=255)
                                                                                       aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc             >  1:354485/1‑100 (MQ=255)
                                                                                        agaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa            >  1:798649/3‑100 (MQ=255)
                                                                                        aaaTCCCCCCGAGCGGGGGGATCTCCAAACACTTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa            <  2:606764/100‑1 (MQ=255)
                                                                                        aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa            <  2:389/100‑1 (MQ=255)
                                                                                        aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa            <  2:1258878/100‑1 (MQ=255)
                                                                                         gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGCTTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag           >  1:3506994/2‑100 (MQ=255)
                                                                                           tctccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagat         <  2:201405/97‑1 (MQ=255)
                                                                                           tCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACCGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag           <  1:1012952/98‑1 (MQ=255)
                                                                                              ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagataga      >  1:2409610/1‑100 (MQ=255)
                                                                                                ccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCACACTGCTCGTTCAGTACTTCACCCACCGCCAGATAGAtt    <  1:2721540/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGC  >  NC_000913/9808‑10005

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: