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

Only 100 of 151 total aligned reads displayed.

GCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCG  >  NC_000913/9809‑10006
                                                                                                  |                                                                                                   
gCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAg                                                                                                    >  1:1718636/1‑100 (MQ=255)
    gAAGGGCCATACGTTGAACCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                                <  2:1210340/100‑1 (MQ=255)
     aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                               <  2:829231/100‑1 (MQ=255)
     aaGGGCCATACGTAGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                               <  2:2274834/100‑1 (MQ=255)
      aGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggggg                                                                                              >  1:2705112/1‑100 (MQ=255)
      aGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATGCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                               <  1:2546031/99‑1 (MQ=255)
        ggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACTTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGa                                                                                            <  2:257300/100‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAt                                                                                           <  2:2763641/100‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                           <  2:841215/100‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                           <  2:586458/100‑1 (MQ=255)
          ccATCCGTAGAAACGGCACCGGAAGTGGTTGCAGCATTAAGACCCAAGAGTGCCAGACGCGACGTTAGCGAATAACAAAATCCCCCCGAGCGGGGGGAtc                                                                                          <  1:1311831/100‑1 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                          <  1:1093718/100‑1 (MQ=255)
          ccATACGGTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                          <  1:1688611/100‑1 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                         <  2:1933633/100‑1 (MQ=255)
            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGCAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                         <  1:2771591/99‑1 (MQ=255)
             tACGTTGAAACGGCACCGGAAGTGGTTGTAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAACCCCCCCGCTCGGGGGGAtttt                                                                                        <  1:2907039/99‑4 (MQ=255)
             tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                       >  1:3045937/1‑100 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGCATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                      <  1:2237970/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                      <  2:216981/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                      <  2:2722608/100‑6 (MQ=255)
               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttttttt                                                                                     <  2:2621697/100‑7 (MQ=39)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                    <  1:795406/99‑1 (MQ=255)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                   <  2:1859771/100‑1 (MQ=255)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                   <  1:2461285/100‑1 (MQ=255)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATGCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                   <  1:1488721/100‑1 (MQ=255)
                   gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTACAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACaa                                                                                 >  1:3314302/1‑100 (MQ=255)
                    aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAt                                                                                <  2:1943805/100‑1 (MQ=255)
                      aCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGCCGTTAGCGAATACAAAAATCCCCCCGAGCGGGGGGATCCCAAAACAAtt                                                                               >  1:264284/1‑99 (MQ=255)
                      aCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTa                                                                              >  2:1065949/1‑100 (MQ=255)
                          cACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTACAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                          <  2:3314330/100‑1 (MQ=255)
                           aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                         <  1:435704/100‑1 (MQ=255)
                           aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                         <  1:2109996/100‑1 (MQ=255)
                              ggAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGCGAtt                                                                      >  1:1926000/1‑100 (MQ=255)
                               gAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAtt                                                                      <  1:3143694/99‑1 (MQ=255)
                                aaGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCa                                                                    <  2:2010581/100‑1 (MQ=255)
                                 aGTGGTTGTAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                   <  2:2707321/100‑1 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  >  1:197654/1‑100 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                  <  2:3374850/100‑1 (MQ=255)
                                   tGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCa                                                                 >  1:1632842/1‑100 (MQ=255)
                                    ggTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCaa                                                                >  1:1639956/1‑100 (MQ=255)
                                       tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                             >  2:1904523/1‑100 (MQ=255)
                                       tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                             >  2:1612061/1‑100 (MQ=255)
                                        gcagcaTTCAGACCGAAGAGGGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                            >  1:3356373/1‑100 (MQ=255)
                                         tccgatctaGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                           <  1:1007597/92‑1 (MQ=255)
                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                           >  1:67603/1‑100 (MQ=255)
                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                           >  2:2465932/1‑100 (MQ=255)
                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                           <  2:788785/100‑1 (MQ=255)
                                          agcaTTCAGCCCGAAGAGTGCAAGACGCGACGTTAGCCAATAAAAAAATCCCCCCGAGAGGGGGGAACTCAAAGCCATTAGTGAGATTCACCCATCgcc                                                           >  1:2955043/1‑97 (MQ=255)
                                          agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                          >  2:1757544/1‑100 (MQ=255)
                                              ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                      <  2:67603/100‑1 (MQ=255)
                                               tCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACg                                                     >  1:955397/1‑100 (MQ=255)
                                                cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                    >  1:2740843/1‑100 (MQ=255)
                                                cAGACCGAAGAGTGCAAGACGCGACGTTACCGAATAAAAAAATCCCCCCGAGCGGGGGGATTTCAAAACAATTAGTGGGCTGCACCAATCGACAGAACg                                                     >  1:2860993/1‑99 (MQ=255)
                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGATCgga                                                   >  2:1007600/1‑99 (MQ=255)
                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGt                                                   >  1:2105865/1‑100 (MQ=255)
                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                  <  1:2759381/100‑1 (MQ=255)
                                                   aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                                 <  2:401232/100‑1 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCGCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                >  1:1368519/1‑100 (MQ=255)
                                                     cGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCTCAGAGCGGGGGGATCTCACCACAATTAGTGGGATTCACCAATCGGCAGACCGCTGCg                                                >  1:630295/1‑99 (MQ=255)
                                                        agagTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGGGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCa                                            <  2:842750/100‑1 (MQ=255)
                                                          agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGGGGGATTCACCAATCGGCAGAACGGTGCGACCaaa                                          <  1:2204409/100‑1 (MQ=255)
                                                           gTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                         <  2:3262263/100‑1 (MQ=255)
                                                            tGCAAGACGCGACGTTAGCGAATAAAAGAATCCCCCCGAGCGGGGGGATCTCAAAACAATGAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACt                                        <  2:1414595/100‑1 (MQ=255)
                                                             gCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTg                                       >  2:3080374/1‑100 (MQ=255)
                                                                aGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                    <  1:524007/100‑1 (MQ=255)
                                                                 gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                   <  2:1106801/100‑1 (MQ=255)
                                                                   cgcgACGTTAGCGGATAAAAAAATCCACCCGATGGGTGGGAGGTCAAAACAAGTAGTTGGATGCACCAATCGTGAGAACGGTGCGACCAAACTGCTCg                                   >  1:1125670/1‑98 (MQ=255)
                                                                    gcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                <  1:2742559/100‑1 (MQ=255)
                                                                        cGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                            >  2:470931/1‑100 (MQ=255)
                                                                         gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                           >  1:2414155/1‑100 (MQ=255)
                                                                            aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                        <  2:2270302/100‑1 (MQ=255)
                                                                              cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                      <  1:1888176/100‑1 (MQ=255)
                                                                              cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                      <  2:513639/100‑1 (MQ=255)
                                                                               gAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                     >  2:3139132/1‑100 (MQ=255)
                                                                               gAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                     <  2:1346639/100‑1 (MQ=255)
                                                                               gAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                     <  2:508752/100‑1 (MQ=255)
                                                                                aaTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                     <  1:714270/99‑1 (MQ=255)
                                                                                 aTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                   >  2:642699/1‑100 (MQ=255)
                                                                                  ttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                  >  2:3006888/8‑100 (MQ=38)
                                                                                   tttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  2:755270/7‑100 (MQ=39)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAGCCATCg                >  2:1636119/1‑100 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  2:1481839/1‑99 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  2:1637569/1‑100 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                >  1:2361377/1‑100 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                <  2:897124/100‑1 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAAGTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                <  1:875915/100‑1 (MQ=255)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATTTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc               >  2:1632804/5‑100 (MQ=38)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACGCCTCGGCAGCACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                 >  1:2005934/5‑98 (MQ=255)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCCCCCATCGc               >  2:936098/5‑100 (MQ=255)
                                                                                     aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc               >  1:616777/1‑100 (MQ=255)
                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCGATGGGCAGAACGGTGCGACCAAACCGCTCGGCCCGTACTTCAccc                    >  1:3048217/4‑94 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  2:765293/2‑100 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag            >  1:3342727/2‑100 (MQ=255)
                                                                                         aTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCaga           >  1:1428493/1‑100 (MQ=255)
                                                                                         aTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTGCAGTACTTCACCCATCGCCaga           >  1:2801786/1‑100 (MQ=255)
                                                                                           ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata         <  2:279196/100‑1 (MQ=255)
                                                                                              cccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatagat      <  2:1293470/100‑1 (MQ=255)
                                                                                               ccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt     <  2:1687305/100‑1 (MQ=255)
                                                                                                  aGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcg  >  2:2473363/1‑100 (MQ=255)
                                                                                                  aGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcg  >  1:1133085/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: