Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,847,139 A→G G42G (GGT→GGC hycC ← membrane‑spanning protein of hydrogenase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,847,1390AG98.9% 872.7 / ‑4.3 277G42G (GGT→GGChycCmembrane‑spanning protein of hydrogenase 3
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (116/158);  total (118/159)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

TTGCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAACCA  >  NC_000913/2847049‑2847228
                                                                                          |                                                                                         
ttGCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGcc                                                                                         <  1:1732820/93‑1 (MQ=255)
  gCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCaa                                                                                       <  2:7054957/93‑1 (MQ=255)
  gCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCaa                                                                                       >  2:2174489/1‑93 (MQ=255)
   cacaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAc                                                                                      <  1:495726/93‑1 (MQ=255)
   cacaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAc                                                                                      <  2:5513583/93‑1 (MQ=255)
    acaTCGTAGCTTCCCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAcc                                                                                     <  2:678253/93‑1 (MQ=255)
    acaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAcc                                                                                     >  2:4268600/1‑93 (MQ=255)
     caTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcg                                                                                    <  2:2233283/93‑1 (MQ=255)
      aTCGTAGCTTACCAGCGACAGGGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgc                                                                                   >  1:2701865/1‑93 (MQ=255)
      aTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgc                                                                                   <  2:106336/93‑1 (MQ=255)
      aTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgc                                                                                   >  1:3585254/1‑93 (MQ=255)
       tCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgcg                                                                                  >  2:479487/1‑93 (MQ=255)
       tCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgcg                                                                                  <  2:2462813/93‑1 (MQ=255)
       tCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgcg                                                                                  >  2:517924/1‑93 (MQ=255)
       tCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgcc                                                                                  >  2:1844496/1‑92 (MQ=255)
        cGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCgcgc                                                                                 >  2:3776178/1‑93 (MQ=255)
         gTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGcgcc                                                                                <  2:6290003/93‑1 (MQ=255)
          tAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGcgccg                                                                               >  1:4669186/1‑93 (MQ=255)
           aGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGcgccgc                                                                              <  2:3354468/93‑1 (MQ=255)
           aGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGcgccgc                                                                              >  1:2865181/1‑93 (MQ=255)
            gCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCgccgcc                                                                             >  1:1691626/1‑93 (MQ=255)
            gCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCgccgcc                                                                             <  1:5766795/93‑1 (MQ=255)
              ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           <  1:7220786/93‑1 (MQ=255)
              ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           >  2:4108678/1‑93 (MQ=255)
              ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCCCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           <  2:5094566/93‑1 (MQ=255)
               tACCAGCGACAGCGCACCGCTCACGCCAACCGCGCGAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAt                                                                          <  2:6381875/93‑1 (MQ=255)
               tACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAt                                                                          >  1:7610615/1‑93 (MQ=255)
               tACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAt                                                                          <  1:1322206/93‑1 (MQ=255)
               tACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAt                                                                          <  2:541683/93‑1 (MQ=255)
                   aGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCa                                                                      <  1:4195289/93‑1 (MQ=255)
                   aGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCa                                                                      >  2:2642898/1‑93 (MQ=255)
                   aGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCa                                                                      >  1:2020692/1‑93 (MQ=255)
                    gCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAg                                                                     <  1:4290185/93‑1 (MQ=255)
                    gCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAg                                                                     >  1:298430/1‑93 (MQ=255)
                       aCAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTa                                                                  >  1:918784/1‑93 (MQ=255)
                       aCAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTa                                                                  >  2:7610527/1‑93 (MQ=255)
                         aGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATc                                                                <  2:3019056/93‑1 (MQ=255)
                          gcgcACCGCTCACGCAAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATcc                                                               <  2:1351551/93‑1 (MQ=255)
                             ccccGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGc                                                            <  1:7852379/91‑1 (MQ=255)
                             cACCGCTCACGCCAACCGCGCCAGTCAGTACAGAGAAGCCCGCGGCTGCCGTATACAGATTGCCAACCGCGCCGCCAATTCCAGCTATCCAGc                                                            <  2:2847018/93‑1 (MQ=255)
                              aCCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGcc                                                           >  2:4664679/1‑93 (MQ=255)
                               ccGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCa                                                          >  1:1193323/1‑93 (MQ=255)
                                cGCTCACCCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCAc                                                         <  2:6643339/93‑1 (MQ=255)
                                 ccgcACGCCACCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACt                                                        <  2:826369/90‑1 (MQ=255)
                                   tCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTa                                                      >  1:3251477/1‑93 (MQ=255)
                                   ccACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTa                                                      <  1:5141269/92‑1 (MQ=255)
                                    cACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGCCTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTaa                                                     <  2:7881710/93‑1 (MQ=255)
                                       cccACCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGc                                                  <  2:1500430/92‑1 (MQ=255)
                                        ccAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCt                                                 >  1:6494109/1‑93 (MQ=255)
                                         cAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtt                                                >  1:7528088/1‑93 (MQ=255)
                                         cAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtt                                                <  1:1022631/93‑1 (MQ=255)
                                         cAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtt                                                <  2:7163509/93‑1 (MQ=255)
                                           aCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtttt                                              >  1:5565024/1‑93 (MQ=255)
                                            ccGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCttttt                                             <  2:4451691/93‑1 (MQ=255)
                                              ccgcCCGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTga                                           <  1:7365975/92‑1 (MQ=255)
                                               cgcCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaa                                          <  2:6032557/93‑1 (MQ=255)
                                                 ccAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaa                                        <  2:2697239/93‑1 (MQ=255)
                                                  cagtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaag                                       >  2:2888055/1‑93 (MQ=255)
                                                  cagtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaag                                       >  1:1885610/1‑93 (MQ=255)
                                                    gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     >  1:3469370/1‑93 (MQ=255)
                                                     tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    <  2:3509773/93‑1 (MQ=255)
                                                     tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    >  1:41627/1‑93 (MQ=255)
                                                     tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    <  1:1566097/93‑1 (MQ=255)
                                                      cagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaaa                                   <  2:1564483/93‑1 (MQ=255)
                                                      cagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaaa                                   <  2:2019366/93‑1 (MQ=255)
                                                       agtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGaaaagaaaag                                  <  2:142414/93‑1 (MQ=255)
                                                       agtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGaaaagaaaag                                  <  1:2140684/93‑1 (MQ=255)
                                                        gtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAaaagaaaaga                                 <  2:4887186/93‑1 (MQ=255)
                                                            aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  1:5373541/93‑1 (MQ=255)
                                                            aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  2:5753722/93‑1 (MQ=255)
                                                            aGTGAAGCCCGCGGCTCCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  2:1884996/93‑1 (MQ=255)
                                                              tGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAAtc                           >  1:6714055/1‑92 (MQ=255)
                                                              tGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATg                           <  1:3971962/93‑1 (MQ=255)
                                                                aaGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGcc                         >  2:1261135/1‑93 (MQ=255)
                                                                 aGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCa                        <  2:1743768/93‑1 (MQ=255)
                                                                    ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     >  1:2447441/1‑93 (MQ=255)
                                                                    ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  2:5927482/93‑1 (MQ=255)
                                                                    ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     >  2:1702329/1‑93 (MQ=255)
                                                                    ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     >  1:4575888/1‑93 (MQ=255)
                                                                      gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   <  2:4713547/93‑1 (MQ=255)
                                                                      gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   <  1:12641/93‑1 (MQ=255)
                                                                      gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   <  1:2500705/93‑1 (MQ=255)
                                                                      gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   >  2:3672551/1‑93 (MQ=255)
                                                                       cgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAg                  <  1:7177328/93‑1 (MQ=255)
                                                                       cgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCAGAACAgc                 <  2:2334871/93‑1 (MQ=255)
                                                                        ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  1:223440/93‑1 (MQ=255)
                                                                        ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  2:2656749/93‑1 (MQ=255)
                                                                        ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  2:7363752/93‑1 (MQ=255)
                                                                         gCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcg                >  2:4857569/1‑93 (MQ=255)
                                                                              cGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcggc           <  1:4258635/93‑1 (MQ=255)
                                                                               gTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcggcg          >  1:597224/1‑93 (MQ=255)
                                                                                tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         <  1:3529425/93‑1 (MQ=255)
                                                                                tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         <  1:707451/93‑1 (MQ=255)
                                                                                 ataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGAc        >  1:4895091/1‑93 (MQ=255)
                                                                                  taCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACa       >  1:4426939/1‑93 (MQ=255)
                                                                                  taCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACa       >  1:1481753/1‑93 (MQ=255)
                                                                                    cAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     <  2:3925955/93‑1 (MQ=255)
                                                                                     aGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAc    >  2:5794436/1‑93 (MQ=255)
                                                                                     aGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAc    >  1:6188862/1‑93 (MQ=255)
                                                                                      gACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcc   <  1:5212350/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTGCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAACCA  >  NC_000913/2847049‑2847228

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